Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,173 T→C 6.2% intergenic (+92/+27) ghoT → / ← lysU toxin of GhoTS toxin‑antitoxin pair; membrane‑lytic protein; stimulator of persister cell formation/lysine tRNA synthetase, inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,353,1730TC 6.2% 309.8 / 8.1 113intergenic (+92/+27)ghoT/lysUtoxin of GhoTS toxin‑antitoxin pair; membrane‑lytic protein; stimulator of persister cell formation/lysine tRNA synthetase, inducible
Reads supporting (aligned to +/- strand):  ref base T (54/50);  new base C (2/5);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 116 total aligned reads displayed.

GTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATT  >  NC_000913/4353025‑4353321
                                                                                                                                                    |                                                                                                                                                    
gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGTCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCt                                                                                                                                                      <  2:319652/149‑1 (MQ=255)
  cGGCATAACCTGGCCAATGCGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGa                                                                                                                                                    <  1:685986/149‑1 (MQ=255)
  cGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGa                                                                                                                                                    <  1:541041/149‑1 (MQ=255)
      aTAACCTCGCCAATGAGTCTGCCTGGGGCATTACTTTATTCTCTCATTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGGGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGc                                                                                                                                                <  1:253605/149‑1 (MQ=255)
       tAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCt                                                                                                                                               <  2:719508/149‑1 (MQ=255)
         aacTGGCCAATGAGTCTGCCTGTGGCATTACTTTTATCTCTCTTTTAGGAGCTGTAGGTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTCGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTc                                                                                                                                             <  1:788180/147‑1 (MQ=255)
             ggCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTc                                                                                                                                         >  2:692552/1‑149 (MQ=255)
             ggCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCGTGTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCGCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGtta                                                                                                                                         >  2:216900/1‑148 (MQ=255)
              gCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCGTTCa                                                                                                                                        <  2:396491/149‑1 (MQ=255)
               ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCGTTCAt                                                                                                                                       <  1:1012135/149‑1 (MQ=255)
                 aaTGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTa                                                                                                                                     >  2:457503/1‑149 (MQ=255)
                   tGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTaaa                                                                                                                                   <  1:965384/149‑1 (MQ=255)
                      gTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTTGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAACCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTg                                                                                                                                <  1:324901/149‑1 (MQ=255)
                       tCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGa                                                                                                                               >  1:997941/1‑149 (MQ=255)
                          gCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGAATGCTTCGTTCATTAAAGTGAAat                                                                                                                            >  1:988736/1‑149 (MQ=255)
                             tgtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGGTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAattta                                                                                                                         <  2:645600/149‑1 (MQ=255)
                             tgtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGATGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGGTCATTAAAGTGAAattta                                                                                                                         >  2:801604/1‑149 (MQ=255)
                              gtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGTATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAATGTGAAatttat                                                                                                                        >  1:752028/1‑149 (MQ=255)
                              gtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatttat                                                                                                                        >  1:252062/1‑149 (MQ=255)
                              gtgGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatttat                                                                                                                        >  2:866247/1‑149 (MQ=255)
                                ggCATTACTTTTTTCTCTCTTTTCGGAGCTGTAGTTGAGGGGCCCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCCCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatttattt                                                                                                                      <  1:731589/149‑1 (MQ=255)
                                 gCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTc                                                                                                                     >  2:200324/1‑149 (MQ=255)
                                   aTTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCACCTCCAAACTGGATTTGCGTAACGATGAACGAAGCAGTCGCCGGAATGCTTCGTTCATTAAAGTGAAATTTATTTCtg                                                                                                                   >  2:858976/1‑149 (MQ=255)
                                     tACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgtg                                                                                                                 <  2:5860/149‑1 (MQ=255)
                                     tACTTTTTTCTCACTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTTAAACGAAGCGTCTATCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGAGAAATTTATTTCtgag                                                                                                                 >  2:296962/1‑147 (MQ=255)
                                        ttttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                              >  2:134477/1‑149 (MQ=255)
                                        ttttgtCTCTCGTTTAGTATCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGATGCGTCTCTACTGTAAAGTGTATTTGCTGATCGTTGATAGAATCAGTCGGATGACTGCTTCTTGAATTTAAGTGTAATTTATTTATGTTGgc                                                                                                              >  2:671350/1‑149 (MQ=255)
                                          ttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgcgc                                                                                                            >  2:388851/1‑149 (MQ=255)
                                              ctctctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCg                                                                                                        >  2:147469/1‑149 (MQ=255)
                                               tctctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGc                                                                                                       >  2:206787/1‑149 (MQ=255)
                                               cctctTTTCTGATCTGTAGTTGAGTGGACCCGCTGGCTATCCGGAAACGAAGCGTCTCTCCTCAAACCTGGATTAGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGc                                                                                                       <  2:667836/148‑1 (MQ=255)
                                                 tctTTTAGGAGCTGTAGTTGAGTGGACGCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCg                                                                                                     <  1:948813/149‑1 (MQ=255)
                                                 tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGGAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCg                                                                                                     <  2:115323/149‑1 (MQ=255)
                                                 tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCg                                                                                                     >  2:363301/1‑149 (MQ=255)
                                                 tctTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCg                                                                                                     >  1:245794/1‑149 (MQ=255)
                                                   tttaAGGAGCTGAAGTTGATTGTACCCGCTCGTTTTCGTGAAACGAAGCGTCTATCCTCAATACTGGATTTGAGGAACGATGAACGTAGCAGTCGCCTGACTGCTTCGTTCATTAAATTGAAATTTATTTCTGTGGGCGCATCGCCggg                                                                                                   <  1:858123/149‑1 (MQ=255)
                                                     ttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGGTGATTTAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGaa                                                                                                 >  1:134941/1‑149 (MQ=255)
                                                       aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                               <  2:154527/149‑1 (MQ=255)
                                                       aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                               <  1:1013642/149‑1 (MQ=255)
                                                       aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                               <  1:451883/149‑1 (MQ=255)
                                                              gttgTTGAGTGGACCAGCTGGGTATCTTGCAACGCAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGCTGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAAc                                                                                        <  2:722747/146‑1 (MQ=255)
                                                              gTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTTATGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAAc                                                                                        >  2:945110/1‑149 (MQ=255)
                                                                   tGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTcgc                                                                                   <  1:895292/149‑1 (MQ=255)
                                                                         gACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAg                                                                             >  2:141226/1‑149 (MQ=255)
                                                                          aCCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATATGCGGAACGATGATCGAAGCAGTCGCCTGACTGCGTCGTTCATTAAAGGGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGATTAACGTCGCGAATAGt                                                                            <  2:414618/149‑1 (MQ=255)
                                                                            ccGCTGGTTATCGTGATACGAAGCGTCACTCCTCAAAACTGGATTTGCGGAACGATGCACGAAGCAGTCGCCGGACGGCTTCGTTAAGTAAAGTGAAATTTATGTCTGTGGTCGCAGCGCCGCGAACTGAATAACGTCGCGAAGAGTaa                                                                          >  2:248173/1‑148 (MQ=255)
                                                                            ccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGCAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTAt                                                                          >  1:189048/1‑149 (MQ=255)
                                                                            ccGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTAt                                                                          >  1:217171/1‑149 (MQ=255)
                                                                                   ttATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGtt                                                                   <  2:44320/149‑1 (MQ=255)
                                                                                   ttATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGtt                                                                   >  1:84955/1‑149 (MQ=255)
                                                                                     aTCGTGAAACGAAGCGTCTCTCCTCATAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGAGTGGTTCGTTGATTTATGTGATTTTTATTTGTGTGGGCGCATCGCCGGGAGGTGACTACCGTGGCGATTTGTATGGCGGtttg                                                                 >  1:244236/1‑147 (MQ=255)
                                                                                     aTCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAg                                                                 >  1:114792/1‑149 (MQ=255)
                                                                                        gTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGATCGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCCTCGCCGGGAAGAGAAAAACGTCGCGAATTGTATTGCTGATAGTaa                                                              >  2:420469/1‑149 (MQ=255)
                                                                                        gTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaa                                                              >  2:168865/1‑149 (MQ=255)
                                                                                        gTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGGTCATTAAAGTGAAAGTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGGaa                                                              >  1:210585/1‑149 (MQ=255)
                                                                                         tGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                             >  2:211015/1‑149 (MQ=255)
                                                                                         tGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                             >  1:48479/1‑149 (MQ=255)
                                                                                          gAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCTGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAAc                                                            <  1:857594/149‑1 (MQ=255)
                                                                                          gAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAAc                                                            <  2:63706/149‑1 (MQ=255)
                                                                                           ataCGAATCGTATCTACTCAAAACTGGATTTGCGGAAAGATGAAAGAAGCAGTCTCCTGAATGCTTCGTTAATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGAGAATAGTATGGCTGTTAGTAAACa                                                           <  2:116802/147‑1 (MQ=255)
                                                                                           aaaCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACa                                                           <  2:885190/149‑1 (MQ=255)
                                                                                            aaCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAg                                                          >  2:507193/1‑149 (MQ=255)
                                                                                            aaCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAg                                                          >  2:956272/1‑149 (MQ=255)
                                                                                             aCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGc                                                         <  1:206787/149‑1 (MQ=255)
                                                                                                aaGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATa                                                      >  2:525717/1‑149 (MQ=255)
                                                                                                aaGAGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATa                                                      <  2:953910/149‑1 (MQ=255)
                                                                                                 aGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATATCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATaa                                                     >  1:971504/1‑149 (MQ=255)
                                                                                                        ctcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCg                                              <  2:417829/149‑1 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGTGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTTTTAGTAAAGCGCATAATCCTTAgg                                             >  1:575526/1‑149 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAAACATTCgg                                             >  2:847809/1‑149 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCGTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGATTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCgg                                             >  1:578133/1‑149 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCATGACTGAATCGTTCATGAAAGAGAAATTTATTTCTGTGGGCGCAACGCCTGGAATAGAATAACGTCGCGAATTGTATGGCTGTTAGTAAACAGCAAAATCAATCgg                                             >  2:711874/1‑149 (MQ=255)
                                                                                                           ctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTc                                           <  2:561536/149‑1 (MQ=255)
                                                                                                            tcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                                          >  1:883274/1‑149 (MQ=255)
                                                                                                            tcAAAACTGGATATGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCg                                          >  1:328243/1‑149 (MQ=255)
                                                                                                                  cTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCa                                    >  2:359677/1‑149 (MQ=255)
                                                                                                                   tGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCAGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATAGCaa                                   >  1:616511/1‑149 (MQ=255)
                                                                                                                    ggATTTGCGGTACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGGGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAt                                  <  2:962503/149‑1 (MQ=255)
                                                                                                                      ttttGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAc                                <  1:330900/148‑1 (MQ=255)
                                                                                                                       tttGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                               <  1:507193/149‑1 (MQ=255)
                                                                                                                       tttGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                               <  1:956272/149‑1 (MQ=255)
                                                                                                                         tGCGGAACGATGAACGAAGCAGTCGCCTTACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCAGCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                             <  2:931663/149‑1 (MQ=255)
                                                                                                                          gCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAAt                                                    <  1:181487/125‑1 (MQ=255)
                                                                                                                          gCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAg                            >  1:966607/1‑149 (MQ=255)
                                                                                                                            actttaaTGAACGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTATTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAc                          >  1:668382/7‑149 (MQ=38)
                                                                                                                              tttaaTGAACGAAGCAGTCGCCTGCCTGCTTCGGTGATTAAAGTGAAATTTTTTTCTGTGGGCGCTTCGCGGGGTAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCCTTCGGTCGATCCCAATACCCTGAccc                        >  2:958619/5‑148 (MQ=255)
                                                                                                                              tttaaTGAACGAAGCAGTCAGGCGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATACCGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCa                        >  1:767959/5‑149 (MQ=39)
                                                                                                                              aaCGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACCGCATAATCATTCGGTCGATACCAATACCCAGACCa                        <  1:299263/149‑1 (MQ=255)
                                                                                                                                 gATGAACGAAGCAGTCCCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCg                     <  2:285528/149‑1 (MQ=255)
                                                                                                                                      aCGAAGCAGTCGCCTGACTGCGTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTgg                <  1:460464/149‑1 (MQ=255)
                                                                                                                                       agAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggc               <  2:641637/148‑1 (MQ=255)
                                                                                                                                        gAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcg              >  1:132232/1‑149 (MQ=255)
                                                                                                                                          aGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTggcggc            <  1:662155/149‑1 (MQ=255)
                                                                                                                                            cAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGGCGATACCAATACCCAGACCAGCGGTTGGCGGCAg          <  2:894996/149‑1 (MQ=255)
                                                                                                                                             aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGa         <  1:692552/149‑1 (MQ=255)
                                                                                                                                             aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGa         <  1:211015/149‑1 (MQ=255)
                                                                                                                                             aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGa         <  1:457503/149‑1 (MQ=255)
                                                                                                                                                   cTGAGTGCTTCGTTCATTAAAGTGAATTTTATTTCTGTGGGCGCATGGCCGGTAAGAGAATAACGGCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTGGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCatat   <  2:352633/149‑1 (MQ=255)
                                                                                                                                                    tGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATAtt  >  1:639305/1‑149 (MQ=255)
                                                                                                                                                    tGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATAtt  >  1:635657/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATATT  >  NC_000913/4353025‑4353321

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: