Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,828,954 T→C 44.2% S4P (TCC→CCC)  xanP → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,828,9540TC44.2% 50.7 / 221.2 179S4P (TCC→CCC) xanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base T (56/43);  new base C (37/42);  total (93/86)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 179 total aligned reads displayed.

TTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCC  >  NC_000913/3828856‑3829050
                                                                                                  |                                                                                                
ttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTc                                                                                                 >  1:3044112/1‑100 (MQ=255)
 tAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTcc                                                                                                >  2:1550410/1‑100 (MQ=255)
 tAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGACTGTTTTCGATGTCTGTTTcc                                                                                                >  1:1036253/1‑100 (MQ=255)
  aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTTCTACCACGGGATTGTTTTCGATGTCTGTTTCCa                                                                                               >  1:1707045/1‑100 (MQ=255)
  aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCa                                                                                               <  1:4197345/100‑1 (MQ=255)
   aaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCAc                                                                                              <  2:1102559/100‑1 (MQ=255)
     tatTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCAccc                                                                                            <  2:388071/100‑1 (MQ=255)
         tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGa                                                                                        <  1:1117725/100‑1 (MQ=255)
         tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGa                                                                                        >  2:754100/1‑100 (MQ=255)
         tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGACTGTTTTCGATGTCTGTTTCCACCCTCGa                                                                                        >  1:1350757/1‑100 (MQ=255)
          aGATGAGGATGTGGCAAACCTTTGCTTTTCTCTTTCTGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAg                                                                                       <  1:693863/100‑1 (MQ=255)
           gATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGt                                                                                      >  2:2478667/1‑100 (MQ=255)
            aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTc                                                                                     >  1:1424138/1‑100 (MQ=255)
            aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTc                                                                                     <  1:3494693/100‑1 (MQ=255)
            aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTc                                                                                     <  2:1513326/100‑1 (MQ=255)
             tGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCa                                                                                    <  1:3891997/100‑1 (MQ=255)
               aGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGa                                                                                  <  1:566373/100‑1 (MQ=255)
               aGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGa                                                                                  <  2:1007020/100‑1 (MQ=255)
                ggATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCTCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGaa                                                                                 <  1:2335466/100‑1 (MQ=255)
                  aTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGaaaa                                                                               <  1:3302814/100‑1 (MQ=255)
                   tgtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGTTGTCTGTTTCGACCCTCGAGCCAGaaa                                                                                >  1:788494/1‑98 (MQ=255)
                   tgtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAAt                                                                              >  1:3232718/1‑100 (MQ=255)
                    gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATg                                                                             <  2:2283996/100‑1 (MQ=255)
                    gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATg                                                                             <  2:1880759/100‑1 (MQ=255)
                     tgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGTCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgc                                                                            <  2:2561453/100‑1 (MQ=255)
                      ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgcg                                                                           >  1:422954/1‑100 (MQ=255)
                        cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCa                                                                         >  1:729975/1‑100 (MQ=255)
                        cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCa                                                                         >  2:581556/1‑100 (MQ=255)
                             gTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                    >  2:3691487/1‑100 (MQ=255)
                              tctGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                   <  2:282181/98‑1 (MQ=255)
                                tGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTg                                                                 <  1:2051339/100‑1 (MQ=255)
                                 gCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGtaga                                                                >  1:282179/1‑97 (MQ=255)
                                 gCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                <  2:2936/100‑1 (MQ=255)
                                    tttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCCACCGGTTGCGCa                                                             <  2:2808360/100‑1 (MQ=255)
                                    tttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCa                                                             >  1:1873848/1‑100 (MQ=255)
                                       cgctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTGTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattt                                                          <  1:3724512/98‑4 (MQ=255)
                                        tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACt                                                         >  2:3573231/1‑100 (MQ=255)
                                        tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGCGATTGTTGTCGATGTCTGTTCCCACCGTCGAGTCAGAAAATGCGAAACCGGTTGAGCAGACt                                                         <  1:628997/100‑1 (MQ=255)
                                         ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTc                                                        >  1:1569159/1‑100 (MQ=255)
                                         ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTc                                                        >  2:1339320/1‑100 (MQ=255)
                                         ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTc                                                        >  2:1433335/1‑100 (MQ=255)
                                              ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                   >  2:3487647/1‑100 (MQ=255)
                                              ggTCAGCCAAAATGCCCGCTTTGCTACCACGGGACTGTTTTCGTTGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                   <  2:710023/100‑1 (MQ=255)
                                               gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttctgac                                                    >  1:1714458/1‑89 (MQ=255)
                                                 cAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                <  1:3509499/100‑1 (MQ=255)
                                                 cAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                <  1:723606/100‑1 (MQ=255)
                                                   gCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGa                                              >  2:2668975/1‑100 (MQ=255)
                                                     tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAc                                            <  2:2467173/100‑1 (MQ=255)
                                                     tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAc                                            >  2:2750913/1‑100 (MQ=255)
                                                     tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAc                                            <  1:929433/100‑1 (MQ=255)
                                                        aaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGa                                         >  2:3709385/1‑100 (MQ=255)
                                                        aaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGa                                         >  2:866563/1‑100 (MQ=255)
                                                         aTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAt                                        <  2:4252267/100‑1 (MQ=255)
                                                           gcgcGCTTGGCTACCTCGTGATTGTCTCCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGTTTGCGCAGATTCAAAACAGCGAACtta                                         <  1:1949607/94‑3 (MQ=255)
                                                           gCCCGCTTTGCTCCCACGGGCTTGTTTTCGCTGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACATCGAACTGAtt                                       <  1:3998044/99‑1 (MQ=255)
                                                           gCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAttt                                      <  2:4016849/100‑1 (MQ=255)
                                                           gCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAttt                                      >  1:454568/1‑100 (MQ=255)
                                                           gCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAttt                                      >  2:2618639/1‑100 (MQ=255)
                                                            cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTa                                     >  1:3756753/1‑100 (MQ=255)
                                                             ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGCCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAc                                    >  1:3022404/1‑100 (MQ=255)
                                                              cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAcc                                   >  1:2790949/1‑100 (MQ=255)
                                                                  ttGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCt                               <  2:30188/100‑1 (MQ=255)
                                                                    gCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTg                             >  2:1723623/1‑100 (MQ=255)
                                                                    gCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCACAACCGGTTGCGCAGACTCAAAACAGCGGACTGATTTACCGTCTTg                             >  2:3561775/1‑100 (MQ=255)
                                                                    gCTACCACGGGACTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTg                             >  1:295541/1‑100 (MQ=255)
                                                                     cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCACACGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGa                            >  2:183608/1‑100 (MQ=255)
                                                                     cTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGa                            <  1:2607664/100‑1 (MQ=255)
                                                                      tACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                           <  1:2496307/100‑1 (MQ=255)
                                                                        ccACGGGNTTGTNTTCGATATCTGTTCCNACCCNCGAGTCAGNAAATGCGCAACCGGTTGCNCAGACTCAANCNAGNGAACTGATNTACCNTCNTGANNa                         <  1:2226191/100‑1 (MQ=255)
                                                                         cACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                        <  2:4085066/100‑1 (MQ=255)
                                                                         cACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                        >  2:957908/1‑100 (MQ=255)
                                                                         cACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                        <  1:3522973/100‑1 (MQ=255)
                                                                               aTTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGATAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                  >  2:3021560/1‑100 (MQ=255)
                                                                               aTTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccg                  >  2:1786946/1‑100 (MQ=255)
                                                                                ttGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgc                 >  2:330439/1‑100 (MQ=255)
                                                                                 tGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                <  1:2164362/100‑1 (MQ=255)
                                                                                 tGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                <  2:1553729/100‑1 (MQ=255)
                                                                                 tGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                <  1:4075763/100‑1 (MQ=255)
                                                                                  gTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg               >  1:4132611/1‑100 (MQ=255)
                                                                                  gTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg               >  1:1016803/1‑100 (MQ=255)
                                                                                  gTTTTCGATGTCTGTTCCCACCCTCGCGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg               >  1:663656/1‑100 (MQ=255)
                                                                                    tttCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCt             <  2:3267853/100‑1 (MQ=255)
                                                                                     ttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCtt            <  2:1391160/100‑1 (MQ=255)
                                                                                      tCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTc           >  2:1684131/1‑100 (MQ=255)
                                                                                       cGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc          <  1:2150434/100‑1 (MQ=255)
                                                                                       cGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc          >  1:2291422/1‑100 (MQ=255)
                                                                                        gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt         <  2:1350761/100‑1 (MQ=255)
                                                                                        gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt         >  2:3833918/1‑100 (MQ=255)
                                                                                        gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt         >  2:2243092/1‑100 (MQ=255)
                                                                                        gATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt         <  1:2243086/100‑1 (MQ=255)
                                                                                         aTGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTc        <  2:1755731/100‑1 (MQ=255)
                                                                                         aTGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTc        <  2:3788116/100‑1 (MQ=255)
                                                                                           gTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaa      >  1:2028171/1‑100 (MQ=255)
                                                                                            tCTGTTTCCGCCCTCTAGTCAAAAAATGAGCAACCGGTTGCGCAGATTTAAAACAGCTCACTGATTTACCGTCTGGAAGATCGTCCGCCGCATCCTc        <  1:1218458/97‑1 (MQ=255)
                                                                                            tCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa     >  1:1765762/1‑100 (MQ=255)
                                                                                             cTGTTGCCACCCTCGAGTCAGAAAATGCGCAACCGGTGGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc    <  1:1614178/100‑1 (MQ=255)
                                                                                             cTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc    <  1:4250254/100‑1 (MQ=255)
                                                                                              tGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAcc   <  1:1139165/100‑1 (MQ=255)
                                                                                              tGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAcc   >  2:1014222/1‑100 (MQ=255)
                                                                                               gTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAccc  >  1:2626131/1‑100 (MQ=255)
                                                                                                  |                                                                                                
TTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCC  >  NC_000913/3828856‑3829050

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