Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1400AG21.6% 186.4 / 11.0 89intergenic (+123/‑75)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (49/20);  new base G (0/19);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 7.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAC  >  NC_000913/3781042‑3781239
                                                                                                  |                                                                                                     
ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                      <  1:539326/100‑1 (MQ=255)
ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                      <  1:1325857/100‑1 (MQ=255)
ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                      >  1:1380303/1‑100 (MQ=255)
ccTCGCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                      <  2:1003784/100‑1 (MQ=255)
 ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                      <  2:1696428/99‑1 (MQ=255)
    tcCCTCAGCTAGAGGGTCCCCTTCGGCCGATTATTTCCTCCCTCGCCCCGAAACGCCACCGCACTCCTCTTCTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                                                                  <  2:2544205/100‑1 (MQ=255)
     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTGCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                                 >  2:1665816/1‑100 (MQ=255)
     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                                 >  1:2252491/1‑100 (MQ=255)
     cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                                 <  1:1737375/100‑1 (MQ=255)
      ccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                >  1:1385676/1‑100 (MQ=255)
       cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                               >  1:295987/1‑100 (MQ=255)
       cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                               <  1:1131285/100‑1 (MQ=255)
       cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                               >  2:2086946/1‑100 (MQ=255)
       cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                               >  2:1951945/1‑100 (MQ=255)
        tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                                              <  1:610957/100‑1 (MQ=255)
         cAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                                                             >  1:2523194/1‑100 (MQ=255)
           gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCAGGGGAGAGGGCCggg                                                                                               >  2:2155768/1‑96 (MQ=255)
           gggAGAGGAGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                          >  2:238409/1‑100 (MQ=255)
             gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                        >  1:81515/1‑100 (MQ=255)
               gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                      <  1:1138662/100‑1 (MQ=255)
               gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                      <  1:2482769/100‑1 (MQ=255)
               gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                      <  2:2673560/100‑1 (MQ=255)
                 ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                   >  1:200350/1‑100 (MQ=255)
                 ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                   <  2:2045527/100‑1 (MQ=255)
                 ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                   >  1:257849/1‑100 (MQ=255)
                  gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa                                                                                  <  1:2088184/100‑1 (MQ=255)
                   ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                 >  1:2133575/1‑100 (MQ=255)
                    gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                <  1:1741938/100‑1 (MQ=255)
                    gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                <  1:1741815/100‑1 (MQ=255)
                      ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                              <  1:192232/100‑1 (MQ=255)
                       cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                <  1:2198222/97‑1 (MQ=255)
                       cGTTCGGCCCTGTATGTACTCCCTCGCTCTGAAACGACACCGCACGCTTTTTTTCTCCCGCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGcc                                                                             <  2:607765/100‑1 (MQ=38)
                         ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGc                                                                           >  2:423217/1‑100 (MQ=255)
                           ccgcgcTGTATGTACTCCCTCACTCTGAAACCACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                         >  1:1813878/3‑100 (MQ=255)
                            ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACt                                                                        >  2:777541/1‑100 (MQ=255)
                            ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACt                                                                        >  1:306824/1‑100 (MQ=255)
                              cgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGc                                                                      >  2:87275/1‑100 (MQ=255)
                               ncTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCTCTCTTCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                     >  2:1973384/2‑100 (MQ=255)
                                  gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTTTTTTCTCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                  >  2:1875493/1‑100 (MQ=255)
                                      gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                              >  2:931777/1‑100 (MQ=255)
                                      gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                              >  2:2687268/1‑100 (MQ=255)
                                        ccgcccTCACTCGAAAACAACCCCCCTCTTTTTTTTTCCCCCCCCCCCGGGGGGGGGGGGGGCGGGGGTGGGGGAAAAAGGCCGCACTGCTCAcc                                                                 <  1:1768350/92‑1 (MQ=255)
                                        aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                            >  2:1592138/1‑100 (MQ=255)
                                        aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                            >  1:289573/1‑100 (MQ=255)
                                           ccctccctCGAAAGCGACCCCTCTCTTTTTTTTTCCCCCGCGCCCCTGCGGGGGGGGGCCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                          <  1:669988/99‑1 (MQ=255)
                                           ccccccccctGAAGCGACCCCGCCTTTTTTTTTTCCCCCCCGCCCCCCGGGGGGGGGGGCGGGGTGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                         <  1:1114753/92‑1 (MQ=255)
                                           cccTCCCCCTGAAACGACACCGCTCTTTTTTTTTCCCCCTCGCCCCCCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                            <  2:1554204/97‑1 (MQ=37)
                                           cccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                         >  1:853067/1‑100 (MQ=255)
                                            cccccctctGAACCGCCCCGGCCCTCTTTTTTTCCCCCTCCCCCCTCCGGGGAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                        <  2:1016746/95‑1 (MQ=255)
                                            ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                        >  1:2104139/1‑100 (MQ=255)
                                            ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                        >  2:337479/1‑100 (MQ=255)
                                            ccTCACTCTGAAACGACACCACACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCGGGGTGAGGATAAAAGGCCGCACTGCTCACCAATTTTc                                                          >  2:1194264/1‑98 (MQ=255)
                                            ccTCACTCGAAAACGACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                        <  1:1932069/100‑1 (MQ=255)
                                              tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt                                                      >  2:1711085/1‑100 (MQ=255)
                                                cctcgGAACCGACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                       <  2:2018216/96‑1 (MQ=255)
                                                  tctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                  >  1:518909/1‑100 (MQ=255)
                                                  tataaaaCTACCCCGCCCTTTTTTTTTCCCCCCCCCCCCCCGGGGGAGGGGGCCGGGGGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                  <  1:1105424/96‑1 (MQ=255)
                                                   ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                 >  2:1968622/1‑100 (MQ=255)
                                                     gAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                               >  2:395785/1‑100 (MQ=255)
                                                      aaaCCACCCCGCACTTTTTTTTTCCCCCCCGCCCCCCCGGGGAGAGGGCCGTGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                <  2:1495741/98‑1 (MQ=255)
                                                           acacCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                 >  2:1704475/1‑92 (MQ=255)
                                                              ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                      >  1:980923/1‑100 (MQ=255)
                                                              ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                      >  1:2487016/1‑100 (MQ=255)
                                                              ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                      >  1:508464/1‑100 (MQ=255)
                                                                 cACTCTTTTTTCCCCCCTCGCCCCCCGGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                       <  2:1406017/96‑1 (MQ=255)
                                                                   ctctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                 >  1:1175043/1‑100 (MQ=255)
                                                                   ctctTTTTTTCTCCCCCCCCCCCCCGGGGGGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgg                                  <  2:1433862/99‑1 (MQ=255)
                                                                      tttttttCCCCCCCCCCCCCCCGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                 <  1:2157465/97‑1 (MQ=255)
                                                                      tttttttCCCCCCCCCCCCCCCGGGGGGGGGCCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                       <  2:2064707/91‑1 (MQ=255)
                                                                      ttttttcCGCCCCCCCCCCCCCGGGGGGAGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                               <  2:1871177/99‑1 (MQ=255)
                                                                       ttttttCTCCTCCGCCCCGCGGGGGGGGGGGCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                 <  2:124733/96‑1 (MQ=255)
                                                                        tttttCTCCCCCGCCCCCCCGGGGGGAGGGCCGGGGTGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                            <  2:1380291/100‑1 (MQ=255)
                                                                          tttCCCCCCCGCCCCTCCGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAgga                          <  2:306821/100‑3 (MQ=255)
                                                                          tcccccccTCCCCCCGCGGGGGAGGGGGCGGGGTGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTa                          <  1:2130933/95‑1 (MQ=255)
                                                                           tcccccccccccccGCCGGGGGGAGGGCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                         <  1:605424/90‑1 (MQ=255)
                                                                             ctcCCTCGCCCCCCCGGGGGGAGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                       <  2:2069749/100‑1 (MQ=255)
                                                                              tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                      >  1:2458212/1‑100 (MQ=255)
                                                                              tcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                      >  1:2413242/1‑100 (MQ=255)
                                                                               cccTCTCTCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACGAATTTTCTGCTATCCTGCCCCCGCACTAAGGTGGCAGTAttc                       >  2:2282255/1‑96 (MQ=255)
                                                                               cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                     >  2:944564/1‑100 (MQ=255)
                                                                                ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                    >  1:610666/1‑100 (MQ=255)
                                                                                ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                   >  1:934117/1‑100 (MQ=255)
                                                                                   cGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                 >  2:1939190/1‑100 (MQ=255)
                                                                                    gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                >  1:90773/1‑100 (MQ=255)
                                                                                    gCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                >  1:625332/1‑100 (MQ=255)
                                                                                     cccgactggggAGGGGGCGGGTCGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                      <  1:2583137/86‑1 (MQ=255)
                                                                                     ccccggaggggCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACCTc               >  1:1223776/8‑100 (MQ=255)
                                                                                     ccccGCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                   <  1:2445145/96‑1 (MQ=255)
                                                                                      cccccgggggggggggCGGGGGGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg              <  1:1875515/94‑1 (MQ=255)
                                                                                       ccccgcgggAGGGGGCCGGGGTCGGGGTGAAAAGGCCGCACTTCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg              <  1:2078513/94‑1 (MQ=255)
                                                                                       ccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt             >  1:402264/1‑100 (MQ=255)
                                                                                        ccccGGGGAGGGTGCGGGGGGCGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt             <  1:139552/97‑1 (MQ=255)
                                                                                         tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc           >  1:2638336/1‑100 (MQ=255)
                                                                                              ggAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc      <  1:568275/100‑1 (MQ=255)
                                                                                              ggAGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc      <  1:238411/100‑1 (MQ=255)
                                                                                                  agGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc  >  1:426630/1‑100 (MQ=255)
                                                                                                  |                                                                                                     
CCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAC  >  NC_000913/3781042‑3781239

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