Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,353,172 C→G 5.3% intergenic (+91/+28) 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,1720CG 5.3% 289.5 / 9.6 113intergenic (+91/+28)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 C (54/51);  new base G (2/4);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-02

Only 100 of 115 total aligned reads displayed.

GTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATAT  >  NC_000913/4353025‑4353320
                                                                                                                                                   |                                                                                                                                                    
gTCGGAATAACCTGGCCAATGAGTCTGCCTGTGTCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCt                                                                                                                                                     <  2:319652/149‑1 (MQ=255)
 tCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGcccg                                                                                                                                                    >  1:561536/1‑147 (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)
           cTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGt                                                                                                                                          >  2:460464/1‑149 (MQ=255)
             ggCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTc                                                                                                                                        >  2:692552/1‑149 (MQ=255)
              gCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCGTTCa                                                                                                                                       <  2:396491/149‑1 (MQ=255)
               ccAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCAGGCGACTGCTTCGTTCAt                                                                                                                                      <  1:1012135/149‑1 (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)
                           ccTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAatt                                                                                                                          >  2:394783/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                                                                                                                       >  2:866247/1‑149 (MQ=255)
                                 gCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTc                                                                                                                    >  2:200324/1‑149 (MQ=255)
                                   aTTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCACCTCCAAACTGGATTTGCGTAACGATGAACGAAGCAGTCGCCGGAATGCTTCGTTCATTAAAGTGAAATTTATTTCtg                                                                                                                  >  2:858976/1‑149 (MQ=255)
                                    ttACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCtgt                                                                                                                 >  2:662155/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)
                                               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                                                                                                    >  1:245794/1‑149 (MQ=255)
                                                     ttAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGGTGATTTAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGaa                                                                                                >  1:134941/1‑149 (MQ=255)
                                                       aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                              <  2:154527/149‑1 (MQ=255)
                                                       aGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAaga                                                                                              <  1:1013642/149‑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                                                                            >  1:953744/1‑149 (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                                                            >  1:48479/1‑149 (MQ=255)
                                                                                         tGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                            >  2:211015/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)
                                                                                                        ctcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCg                                             <  2:417829/149‑1 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAAACATTCgg                                            >  2:847809/1‑149 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCGTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGATTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCgg                                            >  1:578133/1‑149 (MQ=255)
                                                                                                         tcctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCATGACTGAATCGTTCATGAAAGAGAAATTTATTTCTGTGGGCGCAACGCCTGGAATAGAATAACGTCGCGAATTGTATGGCTGTTAGTAAACAGCAAAATCAATCgg                                            >  2:711874/1‑149 (MQ=255)
                                                                                                          cctcAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTTACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGt                                           <  1:678737/149‑1 (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: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                                                   >  2:181487/1‑125 (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)
                                                                                                                                 gATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCg                    <  1:839170/149‑1 (MQ=255)
                                                                                                                                 gATGAACGAAGCAGTCCCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCg                    <  2:285528/149‑1 (MQ=255)
                                                                                                                                      aCGAAGCAGTCGCCTGACTGCGTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTgg               <  1:460464/149‑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)
                                                                                                                                              gTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCAGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAACCAGACCAGCGCTTGGCGGCAGAc       >  2:800468/1‑149 (MQ=255)
                                                                                                                                                cGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCa     <  2:113185/149‑1 (MQ=255)
                                                                                                                                                   cTGAGTGCTTCGTTCATTAAAGTGAATTTTATTTCTGTGGGCGCATGGCCGGTAAGAGAATAACGGCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTGGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCatat  <  2:352633/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCGGAATAACCTGGCCAATGAGTCTGCCTGTGGCATTACTTTTTTCTCTCTTTTAGGAGCTGTAGTTGAGTGGACCCGCTGGTTATCGTGAAACGAAGCGTCTCTCCTCAAAACTGGATTTGCGGAACGATGAACGAAGCAGTCGCCTGACTGCTTCGTTCATTAAAGTGAAATTTATTTCTGTGGGCGCATCGCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTTGGCGGCAGACCATAT  >  NC_000913/4353025‑4353320

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: