Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,170 G→A 5.3% intergenic (+89/+30) 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,1700GA 5.3% 317.6 / 6.4 115intergenic (+89/+30)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 G (58/49);  new base A (2/4);  total (60/55)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-01

Only 100 of 116 total aligned reads displayed.

CTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCA  >  NC_000913/4353022‑4353317
                                                                                                                                                    |                                                                                                                                                   
cTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCg                                                                                                                                                     >  1:22558/1‑149 (MQ=255)
 tGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATAGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGc                                                                                                                                                    >  2:240632/1‑149 (MQ=255)
   gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGTCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCt                                                                                                                                                  <  2:319652/149‑1 (MQ=255)
    tCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGcccg                                                                                                                                                 >  1:561536/1‑147 (MQ=255)
     cGGCATAACCTGGCCAATGCGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGa                                                                                                                                                <  1:685986/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)
                                           ttttttCTCTCTTTTCGGAGCTGTAGTTGAGTGGACCCGCGGGTTATCGTGAAACGAAGCGTCGCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                          <  1:249031/149‑1 (MQ=255)
                                           ttttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgc                                                                                                          >  2:134477/1‑149 (MQ=255)
                                           ttttgtCTCTCGTTTAGTATCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGATGCGTCTCTACTGTAAAGTGTATTTGCTGATCGTTGATAGAATCAGTCGGATGACTGCTTCTTGAATTTAAGTGTAATTTATTTATGTTGgc                                                                                                          >  2:671350/1‑149 (MQ=255)
                                             ttttCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGgcgc                                                                                                        >  2:388851/1‑149 (MQ=255)
                                                tctctcTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATc                                                                                                     >  2:29453/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)
                                                      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)
                                                                     ttGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGCCTGCTTCGTTCATTAAAGTGATATTTATTTCTGTGGGCGCATCGCCGGGAAGAGATTAACGTcg                                                                                >  1:444809/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                                                               >  1:84955/1‑149 (MQ=255)
                                                                                      ttATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGtt                                                               <  2:44320/149‑1 (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:956272/1‑149 (MQ=255)
                                                                                               aaCGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAg                                                      >  2:507193/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)
                                                                                                                          tttGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                           <  1:956272/149‑1 (MQ=255)
                                                                                                                          tttGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                           <  1:507193/149‑1 (MQ=255)
                                                                                                                            tGCGGAACGATGAACGAAGCAGTCGCCTTACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCAGCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                         <  2:931663/149‑1 (MQ=255)
                                                                                                                            tGCCGAACGATGAACGTCGCAGTCGCCTGACTGCTTAGTTAATTAAAGTGAAATTTATTTCTGTGGGCGCATAGCCGGGAAGAGAAAAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                         <  1:566457/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)
                                                                                                                                         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:457503/149‑1 (MQ=255)
                                                                                                                                                aGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGa     <  1:211015/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGGTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCA  >  NC_000913/4353022‑4353317

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: