Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,828,980 A→C 5.3% Q12H (CAA→CAC xanP → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,828,9800AC 5.3% 950.2 / 2.5 357Q12H (CAA→CACxanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base A (181/147);  new base C (5/24);  total (186/171)
Fisher's exact test for biased strand distribution p-value = 7.41e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

TATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828832‑3829129
                                                                                                                                                    |                                                                                                                                                     
  taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    >  1:3352926/1‑150 (MQ=255)
  taAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCg                                                                                                                                                    >  2:5458766/1‑150 (MQ=255)
   aaaaaaGTGCATTTGTCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                                                                                                   >  1:1667697/1‑150 (MQ=255)
   aaaaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                                                                                                   >  2:5367917/1‑150 (MQ=255)
     aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                 <  1:4739649/150‑1 (MQ=255)
     aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                 >  1:2110580/1‑150 (MQ=255)
     aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                 >  1:4405161/1‑150 (MQ=255)
     aaaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGACAGCTAAAAAGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCATCGAGTCAGAAAATGCGCAACCGGtt                                                                                                                                                 >  2:1507742/1‑150 (MQ=255)
      aaaGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTg                                                                                                                                                <  2:1128767/150‑1 (MQ=255)
          tGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCa                                                                                                                                            <  2:5664316/150‑1 (MQ=255)
            ccgttctctcATTTAAATATTTAGATGAGGAGGTGGCAAACGTTTGCTTTTCTCGGTCGGTCAGCTAAAATGCCCGCTTTGCTACAACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGACAATGCGCCCCCGGGTGCGCAGa                                                                                                                                          <  2:1742408/144‑1 (MQ=255)
             aTTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAc                                                                                                                                         <  2:200161/150‑1 (MQ=255)
                  tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaa                                                                                                                                    >  2:2262448/1‑150 (MQ=255)
                  tctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaa                                                                                                                                    >  2:2941093/1‑150 (MQ=255)
                   ctcATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaaa                                                                                                                                   <  2:2058439/150‑1 (MQ=255)
                      aTTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAg                                                                                                                                >  1:1005043/1‑150 (MQ=255)
                      aTTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAAAAg                                                                                                                                >  1:2382722/1‑150 (MQ=255)
                       tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                                                                                               <  2:6229411/150‑1 (MQ=255)
                       tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCACACTCAAAACAGc                                                                                                                               >  1:2026168/1‑150 (MQ=255)
                                ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                                                                                                      <  2:397774/150‑1 (MQ=255)
                                   gATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAc                                                                                                                   >  2:513875/1‑150 (MQ=255)
                                     tGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAACg                                                                                                                 >  1:3556908/1‑150 (MQ=255)
                                        ggATGTGGCAAACGTGTGCTTTTCTCTTTCGGTCCGCTACAATGCCCGCTTTGCTACCACGGGAGTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCt                                                                                                              <  2:6327683/150‑1 (MQ=255)
                                         gttgtgGCAAACGTTTGCTTTTGTCTGTCGGACCGATAAAAGGCCCGATTGGCTACCCGGTTTTTGTTTTCGTTGTCTGTTTCCACCCTCGAGCAACAAACTGCGCAACGGGTCGCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                                                                                                             <  1:4298662/148‑1 (MQ=255)
                                         gATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCTCAGACTCAAAACAGCGAACTGATTTACCGTCtt                                                                                                             >  2:5581602/1‑150 (MQ=255)
                                            gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAAGGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                                                                                                          <  2:3277563/150‑1 (MQ=255)
                                               gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTTCGCAGACTCAAAAAAGCCAACAGATTTAACGTCTTGaataa                                                                                                       >  1:1525180/1‑147 (MQ=255)
                                                cAATCGTTTGTTTTTCTCTTTCGGTCCGATAAATGGCCCGCTTTGATACCTCGGGTTTGTTTTCGATGTCTGTTTCCACCCTCTAGTCATAAATTGCGCAACGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                      <  1:1259291/150‑1 (MQ=255)
                                                cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                                                                                                      <  2:2805987/150‑1 (MQ=255)
                                                 aacggTTTGCTTTTCTGTGTGGGTCAGCTAAATTGCCCGCTTTGCTACCTGGTGTTTTTTTTTGTTGTCTTTTTCAACCCTCGAGTCACAAATTGCTCACGGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCg                                                                                                     <  1:5185035/146‑1 (MQ=255)
                                                 aaaCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCg                                                                                                     >  2:759939/1‑150 (MQ=255)
                                                 aaaCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGGTGCGCAGACTCAAAAAAGCGAACTGATTTACCGGCTTGAAGATCg                                                                                                     >  1:4034789/1‑150 (MQ=255)
                                                     gTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                                                                                                 >  2:5303509/1‑150 (MQ=255)
                                                      tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                                                                                                >  1:1447739/1‑150 (MQ=255)
                                                           ttttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           <  2:4848880/150‑1 (MQ=255)
                                                           ttttCTCTTGGGGTCAGCTAAAAGGCCCGCTTTGCTCCCCGGGGTTTTTTTGCGTTGTCTGTCCCCCCCCTCGAGTCAGAAACGGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt                                                                                           <  1:6198155/150‑1 (MQ=255)
                                                            tttCTCTTTCGGACAGATAAACTGCCCGCTTTGCTACCACGGTATTTTTTTCGTTGTCTGTTTCCACCCTCGTGTCATAAATTGCGCAACGGTTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc                                                                                          <  1:2693925/150‑1 (MQ=255)
                                                             ttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc                                                                                         <  2:2652178/150‑1 (MQ=255)
                                                             ttCTCTGTCGGGCAGATAAACGGCCCGCTTTGCTACCAGGGTTTTTTTTTCTATGTCTTTACCCCCCCTCGAGTCAGAAAAGGCCCAACCGGTTGCCCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc                                                                                         <  1:6427796/150‑1 (MQ=255)
                                                                 ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaa                                                                                     <  2:3797336/150‑1 (MQ=255)
                                                                  tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCTTCCTCaaa                                                                                    >  1:1099274/1‑150 (MQ=255)
                                                                  tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGGCTTGAAGATCGGCCGCCGCTGccgccaa                                                                                    >  2:4528372/1‑145 (MQ=255)
                                                                  tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAAGCGAACTGATTTACCGGCTTGAAGATCGTCCGCCGCTTCCTCaaa                                                                                    >  1:1881981/1‑150 (MQ=255)
                                                                     cGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAccc                                                                                 <  2:2110580/150‑1 (MQ=255)
                                                                      ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCt                                                                                <  2:6641372/150‑1 (MQ=255)
                                                                         cAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAACATGCGCAACCGGTTGCTCAGACTCAAAACAGCGAACTGATTTAGCGTCTTGAAGATCGTCCGCCGCTTCCTCATACCCTGtt                                                                             <  2:4249181/150‑1 (MQ=255)
                                                                          aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGttt                                                                            >  2:5939384/1‑150 (MQ=255)
                                                                            cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGCACAGATTTACCGTCTTGAAGATCGGCCGCCGCGGCCTCAAACCCGGTTTgc                                                                          >  1:2649423/1‑150 (MQ=255)
                                                                             tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCACAAACCCTGTTTgcc                                                                         >  1:976706/1‑150 (MQ=255)
                                                                             tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGGCTTGAAGATCGTCCGCCGCGGCCACAAACCCTGTTTgcc                                                                         >  1:1808985/1‑150 (MQ=255)
                                                                               aattGCCCGGTTTGCTACCAGGGGTTTGTTTTCGTTGTCTGTTTCCACCCTCGTGTCAGAAATTGCGCACCGGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgc                                                                       <  1:844802/147‑1 (MQ=255)
                                                                               aaaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTgccgc                                                                       >  2:650511/1‑150 (MQ=255)
                                                                                 aTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCt                                                                     <  2:3051907/150‑1 (MQ=255)
                                                                                    cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACACTGTTTGCCGCCTGTc                                                                  >  1:5410785/1‑150 (MQ=255)
                                                                                    cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACACAAAACAGCGAACTGATTTACCGACTTGTAGATCGTCCGCCGCGTCCTCAAAACCAGTTTGCCGCCAGTc                                                                  >  1:4497844/1‑150 (MQ=255)
                                                                                     ccGCTTTGCTACCTCGGGATTGTTTTCGTTGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACGGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCa                                                                 <  1:5681677/150‑1 (MQ=255)
                                                                                     ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACAGGTTGCGCAGACACAAAACAGCGAACTGATTTACCGTCATGAATAACGTCCGCCGCTTCCACAAAACCTGTTTGCCGCCTGtaa                                                                 >  1:4979252/1‑148 (MQ=255)
                                                                                      cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTAAAAACACCGAACAGATTTACCGTCTTGAATATAGTCCGCCGCTTCCTCAAAACCAGTTTTGCGCCAGTCAg                                                                >  1:4611696/1‑150 (MQ=255)
                                                                                              tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgct                                                        >  2:3364081/1‑150 (MQ=255)
                                                                                              tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTGCGCAGACTCAAAACAGCGAACAGATTTACCGTCTTGAAGATCGTGCGCCGCTTCCTCAAACCCTGTTTGGCGCCTGTCAGCATctgat                                                        >  1:6461446/1‑148 (MQ=255)
                                                                                                 cACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGTTTCGCAGAAACAAAACAAAGAAAAGATATAACGTCTTTAATATAGTGCGCCGCGAACACAAAAACAGTTTGCCGCCAGTCAACAAATACTGGc                                                     >  1:710644/1‑150 (MQ=255)
                                                                                                  aCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCg                                                    <  1:3725387/150‑1 (MQ=255)
                                                                                                  aCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCg                                                    >  2:373081/1‑150 (MQ=255)
                                                                                                        ttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATctgc                                                         <  1:2417713/139‑1 (MQ=255)
                                                                                                        ttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTc                                              <  2:3556908/150‑1 (MQ=255)
                                                                                                        ttGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAAACCGAACTGATTTACCGTCTTGAAGATAGTCCGCCGCGTCCACAAACCCAGTTTGCCGCCGGGCAGCAACAGCTGGCGATGTTc                                              >  1:5463941/1‑150 (MQ=255)
                                                                                                            tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACACCGAACTGATTTACCGACTTGAAGATCGTCCGCCGCTTCCTCAAACCCAGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                          >  2:5091370/1‑150 (MQ=255)
                                                                                                            tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTCGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTg                                          <  2:5345345/150‑1 (MQ=255)
                                                                                                               cGATGTCTGTTTCCACCCTCGAGTCATAAATTGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  1:6027013/150‑1 (MQ=255)
                                                                                                               cGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCgg                                       <  2:1831826/150‑1 (MQ=255)
                                                                                                                 aTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTg                                     >  2:1413412/1‑150 (MQ=255)
                                                                                                                      tGTTTCCACCCTCGAGTCAGAAAATGCGCAACCCGGTGCGCAGACTCAAAAAAGCGAACAGATTTAACGTCTTTAATATCGTCCCCCGCGTCCTCAAACCCCGTTTTCCGCCAGTCAGCAACTGCTGGCGATGTTCGTTGCGGTGATCAc                                >  1:4409222/1‑150 (MQ=255)
                                                                                                                      tGTTTCCACCCTCGAGTCAGAAAATGCGAAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCAc                                <  2:5248080/150‑1 (MQ=255)
                                                                                                                        tttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATTTTCGTTGCGGTGATCACGc                              <  2:402343/150‑1 (MQ=255)
                                                                                                                        tttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGc                              >  2:2454698/1‑150 (MQ=255)
                                                                                                                          tCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCa                            <  1:2214210/150‑1 (MQ=255)
                                                                                                                           ccACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAg                           >  2:4609303/1‑150 (MQ=255)
                                                                                                                            cccccTCGTGACACAAATGGCGCACCCGGGCGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgc                          <  1:6338591/148‑1 (MQ=255)
                                                                                                                            cACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgc                          >  2:1070629/1‑150 (MQ=255)
                                                                                                                            cACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCAGCTGGCGATGTTCGTTGCGGTGATCACGCCAgc                          >  2:1287202/1‑150 (MQ=255)
                                                                                                                              cccTCGAGTCAGAAAATGCGCACCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcgc                        <  1:6567169/150‑1 (MQ=255)
                                                                                                                              cccTCGAGACATAAAATGGGCAACGGTGTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAgcgc                        <  1:3979035/150‑1 (MQ=255)
                                                                                                                               ccTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAGTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCt                       >  2:1635871/1‑150 (MQ=255)
                                                                                                                                cTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTa                      <  1:4163736/150‑1 (MQ=255)
                                                                                                                                 tCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTAt                     >  2:4426356/1‑150 (MQ=255)
                                                                                                                                     gTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAAt                 <  2:6571880/150‑1 (MQ=255)
                                                                                                                                      tccgAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGTACTGATTTACCGTCTTGAAGATCGTCCGCCGCATCCTCAAACCCGGGATGCCGCCTGTCAGAATCTGCAGGCGATGTTCGTTGCGGTGATCACGCCCGCGCTATTAATc                >  2:6394439/4‑150 (MQ=255)
                                                                                                                                      tCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc                <  2:4824923/150‑1 (MQ=255)
                                                                                                                                      tCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATc                <  2:6199054/150‑1 (MQ=255)
                                                                                                                                      tCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCGTCCTCAAACCCTGTTTGCCGCCGGTCAGCATCTGCTGGCGATGTTCGGTGCGGGGATCACGCCAGCGCTATTAATc                >  2:2269492/1‑150 (MQ=255)
                                                                                                                                         gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             >  2:4688029/1‑150 (MQ=255)
                                                                                                                                         gAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAATTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGc             >  2:1676117/1‑150 (MQ=255)
                                                                                                                                           aaaGGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCa           <  2:2693218/150‑1 (MQ=255)
                                                                                                                                                cgcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGgcgc      >  2:4975770/1‑150 (MQ=255)
                                                                                                                                                 gcAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCt     <  1:2941093/150‑1 (MQ=255)
                                                                                                                                                  cAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTg    >  1:3857545/1‑150 (MQ=255)
                                                                                                                                                   aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  1:3035423/150‑1 (MQ=255)
                                                                                                                                                   aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   >  2:3890204/1‑150 (MQ=255)
                                                                                                                                                   aaCCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  1:5836587/150‑1 (MQ=255)
                                                                                                                                                   aaCCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTgg   <  2:3141083/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TATAAAAAAGTGCATTTATCTCATTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGTTTGCCGCCTGTCAGCATCTGCTGGCGATGTTCGTTGCGGTGATCACGCCAGCGCTATTAATCTGCCAGGCGCTGGG  >  NC_000913/3828832‑3829129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: