Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,712,337 A→T 5.2% intergenic (+179/‑432) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,712,3370AT 5.2% 703.8 / 17.1 191intergenic (+179/‑432)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base A (95/86);  new base T (6/4);  total (101/90)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 195 total aligned reads displayed.

TTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTT  >  NC_000913/1712074‑1712611
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                  
             gCCCGCAAGCCCCGCTGTGGCTCTTGTATCATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                                                                                                                                                                                                                                            >  2:251527/1‑275 (MQ=255)
                             gtgGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCAGTGATAACCTGTCTTTATCCTCTGTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                                                                                                                                                                                                                                            <  2:126180/275‑1 (MQ=255)
                                                      ttGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCAAAATTTTCTATTCATTCGTGATCAAGATACCGATGCTTTCGCAATGTAAAAACATTGCAAGCATCAGGTTAAACTCATTTAAATTTaca                                                                                                                                                                                                                   >  2:604476/1‑275 (MQ=255)
                                                                  agagAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGTTTAAAATCATTTAAATTTACACACGCAACAAat                                                                                                                                                                                                       >  1:82422/1‑275 (MQ=255)
                                                                             aCATCTGAAGAAAATGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACa                                                                                                                                                                                            <  1:215901/275‑1 (MQ=255)
                                                                             aCATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACa                                                                                                                                                                                            <  2:390490/275‑1 (MQ=255)
                                                                              cATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTTCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTa                                                                                                                                                                                              <  4:591797/272‑1 (MQ=255)
                                                                                              gTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTa                                                                                                                                                                           >  2:510001/1‑275 (MQ=255)
                                                                                                       aTTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGa                                                                                                                                                                  >  2:400606/1‑275 (MQ=255)
                                                                                                           ccccATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGaaa                                                                                                                                                                <  4:509999/273‑1 (MQ=255)
                                                                                                                   aTAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTAAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTa                                                                                                                                                                                              <  4:655441/235‑1 (MQ=255)
                                                                                                                    tAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTAAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCaa                                                                                                                                                                                                                                                                <  2:75829/168‑1 (MQ=255)
                                                                                                                       ccTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa                                                                                                                                                                                                                                                                                 <  5:3009/148‑1 (MQ=255)
                                                                                                                         tttctttATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATc                                                                                                                                                                                                                                                                               >  5:135105/1‑148 (MQ=255)
                                                                                                                         tttctttATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATc                                                                                                                                                                                                                                                                               >  5:366090/1‑148 (MQ=255)
                                                                                                                           tctttATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTa                                                                                                                                                                                    >  4:93774/1‑237 (MQ=255)
                                                                                                                               tATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTacacac                                                                                                                                                                                                                <  4:50593/205‑1 (MQ=255)
                                                                                                                               tATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTacaca                                                                                                                                                                                                                 <  2:415299/204‑1 (MQ=255)
                                                                                                                               tATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACt                                                                                                                                                                                <  2:93775/237‑1 (MQ=255)
                                                                                                                                  ccTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTg                                                                                                                                                                                                                                                                      >  6:901753/1‑148 (MQ=255)
                                                                                                                                   ctctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTc                                                                                                                                                                                                                                                                    >  7:55465/1‑149 (MQ=255)
                                                                                                                                   ctctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTc                                                                                                                                                                                                                                                                    >  8:179329/1‑149 (MQ=255)
                                                                                                                                     ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                                                                                                                                                                                                                                                   >  5:653903/1‑148 (MQ=255)
                                                                                                                                     ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTacac                                                                                                                                                                                                                  >  2:770070/1‑197 (MQ=255)
                                                                                                                                     ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCACAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACa                                                                                                                                                                                            >  2:172281/1‑219 (MQ=255)
                                                                                                                                     ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCACAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACa                                                                                                                                                                                            <  4:172280/219‑1 (MQ=255)
                                                                                                                                     ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATAACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAAAAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGAACTATCTAACTAGTAAAAACGCTa                                                                                                                                    >  1:102547/1‑275 (MQ=255)
                                                                                                                                      ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTc                                                                                                                                                                                                                                                                    >  1:77669/1‑146 (MQ=255)
                                                                                                                                      ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGa                                                                                                                                                           >  4:600144/1‑251 (MQ=255)
                                                                                                                                       tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCa                                                                                                                                                                                                                                                                 <  6:653903/148‑1 (MQ=255)
                                                                                                                                          aaaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCaa                                                                                                                                                                                                                                                                <  3:77668/146‑1 (MQ=255)
                                                                                                                                          aaaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                                                                                                                                                                                                                                             >  7:259324/1‑149 (MQ=255)
                                                                                                                                           aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGAGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                                                                                                                                                                                                                                            <  7:195455/149‑1 (MQ=255)
                                                                                                                                               aTCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaaa                                                                                                                                                                                                                                                        >  7:306830/1‑149 (MQ=255)
                                                                                                                                                  tttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTc                                                                                                                                                                                                                                                                    <  8:881498/134‑1 (MQ=255)
                                                                                                                                                  tttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTc                                                                                                                                                                                                                                                                    >  7:881498/1‑134 (MQ=255)
                                                                                                                                                  tttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAAtt                                                                                                                                                                                                                       >  2:546758/1‑179 (MQ=255)
                                                                                                                                                   ttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                                                                                                                                                     <  5:681738/148‑1 (MQ=255)
                                                                                                                                                    tttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGc                                                                                                                                                                                                                                                                  <  2:636289/134‑1 (MQ=255)
                                                                                                                                                    tttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                                                                                                                                                                                                                                            <  2:277505/156‑1 (MQ=255)
                                                                                                                                                     ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                                                                                                                                                                                                                                                  <  7:135172/149‑1 (MQ=255)
                                                                                                                                                        aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAg                                                                                                                                                                                                                                               <  8:749423/149‑1 (MQ=255)
                                                                                                                                                         cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                                                                                                                                                                                                                                            >  5:9978/1‑119 (MQ=255)
                                                                                                                                                            cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                                                                                                                                                                                                                                                     >  4:674186/1‑139 (MQ=255)
                                                                                                                                                              atgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAg                                                                                                                                                                                                                                         <  7:355777/149‑1 (MQ=255)
                                                                                                                                                                gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCa                                                                                                                                                                                                                                                 <  2:674188/139‑1 (MQ=255)
                                                                                                                                                                gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                                                                                                                                                                                                                                            <  4:272516/144‑1 (MQ=255)
                                                                                                                                                                gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                                         <  4:619902/243‑1 (MQ=255)
                                                                                                                                                                     cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaaaa                                                                                                                                                                                                                                  <  8:55465/149‑1 (MQ=255)
                                                                                                                                                                     cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaaaa                                                                                                                                                                                                                                  >  7:11321/1‑149 (MQ=255)
                                                                                                                                                                     cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaaa                                                                                                                                                                                                                                   <  6:370092/148‑1 (MQ=255)
                                                                                                                                                                       ttGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTg                                                                                                                                                                                                                                                                      >  3:188670/1‑111 (MQ=255)
                                                                                                                                                                       ttGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATc                                                                                                                                                                                                                                >  7:654776/1‑149 (MQ=255)
                                                                                                                                                                          cTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCaa                                                                                                                                                                                                                                                >  2:503578/1‑130 (MQ=255)
                                                                                                                                                                           ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGc                                                                                                                                                                                                                                                                  <  1:188671/111‑1 (MQ=255)
                                                                                                                                                                           ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAttt                                                                                                                                                                                                                            >  8:228883/1‑149 (MQ=255)
                                                                                                                                                                              gCGTTATCATCATGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTaaa                                                                                                                                                                                                                         <  7:752977/149‑1 (MQ=255)
                                                                                                                                                                               cGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTaaa                                                                                                                                                                                                                         <  6:142408/148‑1 (MQ=255)
                                                                                                                                                                                   atcatcAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTac                                                                                                                                                                                                                    <  8:494632/149‑1 (MQ=255)
                                                                                                                                                                                   atcatcAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTa                                                                                                                                                                                                                     >  1:202055/1‑148 (MQ=255)
                                                                                                                                                                                    tcatcaGGTAATCCGTTTGCGGCTAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaa                                                                                                                                                                                                                                                          >  2:599557/1‑110 (MQ=255)
                                                                                                                                                                                    tcatcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACg                                                                                                                                                                                                               >  1:118777/1‑153 (MQ=255)
                                                                                                                                                                                    tcatcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATACCCCCCCTCCcagtgc                                                                                      >  4:422607/1‑274 (MQ=255)
                                                                                                                                                                                    tcatcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATNTAAATTTACACACGCAACAAATATTGACCTACAAAACATCACACTGGCTATTTTTCAGAAACTGGACTATCTCCNTAGTAAAAACGCTAAATAGCATAACATAATNNCTTTAATAAATCACCCCCCCCTcaatt                                                                                       >  3:94363/1‑270 (MQ=255)
                                                                                                                                                                                       tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                                                                                                                                                                                                                                                         <  5:514974/108‑1 (MQ=255)
                                                                                                                                                                                       tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                                                                                                                                                                                                                                                         >  6:514974/1‑108 (MQ=255)
                                                                                                                                                                                       tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTacacac                                                                                                                                                                                                                >  8:534649/1‑149 (MQ=255)
                                                                                                                                                                                       tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTacaca                                                                                                                                                                                                                 <  3:202054/148‑1 (MQ=255)
                                                                                                                                                                                        caGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAATATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGcaac                                                                                                                                                                                                           <  3:118776/153‑1 (MQ=255)
                                                                                                                                                                                         aGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAtt                                                                                                                                                                                                                                                    >  8:664106/1‑111 (MQ=255)
                                                                                                                                                                                           gTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGcaa                                                                                                                                                                                                            >  8:829800/1‑149 (MQ=255)
                                                                                                                                                                                           gTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGca                                                                                                                                                                                                             >  6:479486/1‑148 (MQ=255)
                                                                                                                                                                                           gTAATCCGTTGGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAACACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGAcc                                                                                                                                                                                                >  1:47684/1‑161 (MQ=255)
                                                                                                                                                                                              aTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAat                                                                                                                                                                                                       >  4:763733/1‑151 (MQ=255)
                                                                                                                                                                                               tCCGTTTGCGGATAACCAAAAANGCAGGTTANTTTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACACAACAAATATTGACCTACa                                                                                                                                                                                            <  3:47683/161‑1 (MQ=255)
                                                                                                                                                                                                   tttGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGcaac                                                                                                                                                                                                           <  4:746441/142‑1 (MQ=255)
                                                                                                                                                                                                     tGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAat                                                                                                                                                                                                       >  1:125496/1‑144 (MQ=255)
                                                                                                                                                                                                      gCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGAc                                                                                                                                                                                                 <  7:437928/149‑1 (MQ=255)
                                                                                                                                                                                                         gATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGca                                                                                                                                                                                                             >  5:820997/1‑134 (MQ=255)
                                                                                                                                                                                                         gATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTg                                                                                                                                                                                                   <  3:125495/144‑1 (MQ=255)
                                                                                                                                                                                                                       cAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGa                                                                                                                                                                                                  <  7:505438/131‑1 (MQ=255)
                                                                                                                                                                                                                        aGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTac                                                                                                                                                                                                                    >  4:455448/1‑112 (MQ=255)
                                                                                                                                                                                                                            tAATTGTTTTTTTAATAGCGAAATTTTCTATTCATGCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACAGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                             <  2:714195/275‑1 (MQ=255)
                                                                                                                                                                                                                                  tttttttAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTc                                                                                                                                                                     >  8:129211/1‑149 (MQ=255)
                                                                                                                                                                                                                                   ttttttAATAGCGAAATTTTCTATTCATTCGCGATCAAGATCACGATGCTGTCGCTATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTc                                                                                                                                                                     <  5:260289/148‑1 (MQ=255)
                                                                                                                                                                                                                                            aGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTgg                                                                                                                                                            <  6:770882/148‑1 (MQ=255)
                                                                                                                                                                                                                                             gCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGa                                                                                                                                                           <  6:825112/148‑1 (MQ=255)
                                                                                                                                                                                                                                                         attcattcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACa                                                                                                                                                                                       >  6:28268/1‑108 (MQ=255)
                                                                                                                                                                                                                                                              ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTAttttt                                                                                                                                                                      >  6:313432/1‑120 (MQ=255)
                                                                                                                                                                                                                                                              ttcGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                                         >  7:364519/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               tcGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGAc                                                                                                                                                          <  5:838459/131‑1 (MQ=255)
                                                                                                                                                                                                                                                               tcGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                                         >  6:229986/1‑148 (MQ=255)
                                                                                                                                                                                                                                                               tcGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTaa                                                                                                                                            <  5:403051/145‑1 (MQ=255)
                                                                                                                                                                                                                                                               tcGTGATCTTGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAAc                                                                                                                                        >  8:400638/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 gTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAg                                                                                                           <  4:302424/176‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaaat                                                                                                                                                                                                                                  >  5:769580/1‑54 (MQ=255)
                                                                                                                                                                                                                                                                    aTCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGAcc                                                                                                                                                                                                >  4:239879/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                    aTCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                                         >  6:439551/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                    aTCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTaaaaa                                                                                                                                         <  5:439551/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aTCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTaa                                                                                                                                   >  8:850380/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                  
TTATACCTGCATTGCCCGCAAGCCCCGCTGTGGCTCTTGTATTATTGAAGATCTTTGTGAATACAAAGAGAAAGTTGACATCTGAAGAAAAGGGGTAACACCGATTACCCCATTGATAACCTTTCTTTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTT  >  NC_000913/1712074‑1712611

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