Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1520AG21.1% 176.7 / 16.0 95intergenic (+135/‑63)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (47/28);  new base G (2/18);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 2.58e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGA  >  NC_000913/3781053‑3781243
                                                                                                         |                                                                                            
gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                              >  2:1964002/1‑100 (MQ=255)
 ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                             <  1:93270/100‑1 (MQ=255)
  gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                            <  1:2612107/100‑1 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                          <  1:2778207/100‑1 (MQ=255)
     agGGGACCGTTCGGCGCTGTATGCACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGa                                                                                        <  1:856368/100‑1 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                       >  1:1203314/1‑100 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                                                                       >  1:849270/1‑100 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTCTCCTTGCCCCTCCCGGGTGAGGGCCGGGGGGGGGGGgaaa                                                                                     >  1:2255885/1‑100 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                     >  1:2017473/1‑100 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                     <  2:1886315/100‑1 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                    >  1:192169/1‑100 (MQ=255)
          aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                   >  2:321491/1‑100 (MQ=255)
          aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                   <  1:1964020/100‑1 (MQ=255)
          aCCGTTCGGCGCTGTATGTACTCCCCCCCCCTGAAACGGCACCGCACTCTTTTTGTCTCCCTAGACCGACCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                   <  2:2866474/100‑1 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                                  >  1:1663671/1‑100 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                                  <  1:2060839/100‑1 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCCCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGc                                                                                  <  1:2033296/100‑1 (MQ=255)
             gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                                                                >  2:2805295/1‑100 (MQ=255)
                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCa                                                                              >  2:2871980/1‑100 (MQ=255)
                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                             >  1:888123/1‑100 (MQ=255)
                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                             >  2:1251389/1‑100 (MQ=255)
                 ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACt                                                                            >  2:2233330/1‑100 (MQ=255)
                 ggCGCTGTATGTACCCCCTCACTCGGAAACGACACCGCACTCTTTTTTCCCCCCCCCCCCCCCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                             <  2:180428/99‑1 (MQ=255)
                 ggCGCGGAATGTACTCCCTCACGCTGAACCGACCCCGCACTTTTTTTTTCTCCCGCGCCCCTCCGGGGGGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                             <  1:1076511/99‑1 (MQ=37)
                   cgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGc                                                                          >  1:74421/1‑100 (MQ=255)
                       gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                      >  1:295664/1‑100 (MQ=255)
                         atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                                                    >  2:2150303/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                                                                  >  1:1790276/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTACTCACCAAt                                                                  >  2:501496/1‑100 (MQ=255)
                              cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                               >  1:1726623/1‑100 (MQ=255)
                               tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                              >  1:827612/1‑100 (MQ=255)
                               cccccccctctAAAACAACCCCGCACCTTTTTTTCCACCCCGCCCCCTCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                              <  1:30568/93‑1 (MQ=255)
                                ccctccctCTGAACCGACCCCGCCCTCTTTTTTTCCCCCCCGCCCCCCGGGGGGGGGGCCCGGGTGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                              <  2:1760606/99‑1 (MQ=255)
                                ccctccctCTAAACCGCCCCCGCCTTTTTTTTTCCCCCCCCCCCCCCCGGGGAAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                             <  1:1819183/100‑1 (MQ=255)
                                ccctccctCTAAAACGCCCCCCCCTTTTTTTTTTCCCCCCCCCCCCTCGGGGGAGGGGGCCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                              <  1:1953591/99‑1 (MQ=255)
                                 cccccctctGAAACGCCACGGCACTTTTTTTTGCCCCCTCCCCCCTCGGGGGAGAGGCCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                            <  2:919671/95‑1 (MQ=25)
                                  cTCCCTCGGAAACGACCCCGCCCTCTTTTTGTCTCCCTCCCCCCTCGGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           <  2:1605145/100‑1 (MQ=255)
                                  cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           >  1:1368371/1‑100 (MQ=255)
                                    ccctctGAAACGACCCCGCATTCTTTTTTTCCCCCTCCCCCCTCGGGGGGGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                                         <  2:1198879/98‑1 (MQ=255)
                                    cACTGTGAAACGCCCCCGCACTTTTTTTTTCCCCCTCGCCCCCCCGGGGAGAGGGCGGGGTGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           <  2:187608/98‑1 (MQ=255)
                                      ctgcaaaaCGCCCCCCCCCTTTTTTTTCCCCCCCCGCCCCTCCGGGGAGGGGGCCGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           <  1:1449570/91‑1 (MQ=255)
                                      ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCACCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                       >  2:1627567/1‑100 (MQ=255)
                                        ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                     >  2:1773972/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTCTTTTTTTTTTCCTTGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTACTCACCAATTTTCTGCTATCCt                                                     >  2:924481/1‑99 (MQ=38)
                                         tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    >  1:2455347/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    >  1:740311/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    >  2:666317/1‑100 (MQ=255)
                                          gAAACGACACCGCTCTCTTTTTTTCTCCCTCACCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                   >  1:2043032/1‑100 (MQ=255)
                                           aaaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                  >  2:1949050/1‑100 (MQ=255)
                                           aaaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                  >  2:2183996/1‑100 (MQ=255)
                                            aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                 >  2:1265944/1‑100 (MQ=255)
                                            aaCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                 >  2:1445173/1‑100 (MQ=255)
                                              cGACCCCCCTCTTTTTTTTTCCCCCCCCCCCCTCGGGGGAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                <  1:1949068/99‑1 (MQ=255)
                                                cccccGCCCTTTTTTTTTCCCCCCCCCCCCTCGGGGGGGGGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                               <  1:317129/95‑1 (MQ=255)
                                                 ccccGCTCTCTTTTTTTCCCCCCCCCCCCCCCCGGGGGGAGGCCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                <  1:2872006/95‑1 (MQ=255)
                                                 ccccGCACTCTTTTTGTCTCCCCCTCCCCCCCGGGGAGCGGCCTGGGTTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                             <  2:1304575/97‑1 (MQ=255)
                                                   ccGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                           >  1:1733623/1‑99 (MQ=255)
                                                    ccgcctctTTTTTTCCCCCCCCCCCCCCCGGGGGGGAGGCCCGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                           <  2:74420/95‑1 (MQ=255)
                                                      ccttctTTTTTCCCCCCCCCCCCCCCCGGGGGGAGGGCGGGGTGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                           <  2:1517357/93‑1 (MQ=255)
                                                      ccctctTTTTTCCGCCCTCCCCCCGCGGGGGGGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                          <  2:629626/95‑1 (MQ=255)
                                                        ctctTTTTTTCTTTCTTGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                     >  2:68267/1‑100 (MQ=255)
                                                         ctttttttttCTCCCTCCCCCCGCGGGGGAGAGGGCCGGGTTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                              <  2:192166/88‑1 (MQ=38)
                                                          ctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   >  2:2192909/1‑100 (MQ=255)
                                                          ctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   >  1:2903263/1‑100 (MQ=255)
                                                            tttttttCCCCCTCGCCCCCCCGGGGGGAGGCCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   <  2:695845/99‑1 (MQ=255)
                                                            tttttttCCCCCCCGCCCCTCGGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                  <  2:2115140/100‑1 (MQ=255)
                                                            tttttttCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                     <  1:2184015/98‑1 (MQ=255)
                                                              tttttCCCCCCCGCCCCCCGGGGGAGAGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                <  2:952736/100‑1 (MQ=255)
                                                                     tcCGTCGCCCCTCCGGGGGGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                          <  2:1603293/100‑1 (MQ=255)
                                                                      ccccCGCCCCCCGGGGGAGGGGCCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                          <  1:2740521/99‑1 (MQ=255)
                                                                      cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  1:1912668/1‑100 (MQ=255)
                                                                      cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  1:638415/1‑100 (MQ=255)
                                                                      cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  2:1408672/1‑100 (MQ=255)
                                                                      cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  2:757474/1‑100 (MQ=255)
                                                                       ccgcGCCCCTCGGGGGGGAGGGCCGGGGTGGTGGGAATAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        <  2:2043014/97‑1 (MQ=255)
                                                                       ccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        >  1:2182433/1‑100 (MQ=255)
                                                                       ccTCCCCCCCCGGGGGGGGGGGCGGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        <  1:1673293/100‑1 (MQ=255)
                                                                        cccccccccccCGGGGGGGGGGGCGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        <  1:405354/95‑1 (MQ=255)
                                                                        cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                       >  1:1718698/1‑100 (MQ=255)
                                                                        cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                       >  1:193005/1‑100 (MQ=255)
                                                                             ccccTCCGGGGAGAGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                      <  2:875511/97‑1 (MQ=255)
                                                                             ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                   >  1:541489/1‑100 (MQ=255)
                                                                             ccccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                   >  1:240891/1‑100 (MQ=255)
                                                                               ccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                 >  1:1279061/1‑100 (MQ=255)
                                                                                 cccgggggAGGGGGCGGGGGAGGGAGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                 <  1:1126574/95‑1 (MQ=255)
                                                                                  ccGGGGGGGGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc               <  2:453081/99‑1 (MQ=255)
                                                                                  ccGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt              >  1:2063344/1‑100 (MQ=255)
                                                                                  ccGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt              >  2:1545874/1‑100 (MQ=255)
                                                                                      gggAGGGGGCCGGGGTGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa           <  2:1974232/100‑1 (MQ=255)
                                                                                       ggAGAGGGCCGGGGCGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc          <  2:2063325/100‑1 (MQ=255)
                                                                                        gagaGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg         >  2:2634682/1‑100 (MQ=255)
                                                                                           gagGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa       <  2:125697/100‑1 (MQ=255)
                                                                                             gggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa        <  2:1203303/97‑1 (MQ=255)
                                                                                             gggCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc     >  1:2175468/1‑100 (MQ=255)
                                                                                             ggccgggggggAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc          <  2:629571/90‑1 (MQ=255)
                                                                                               gcgggggggggggAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg   <  1:2233350/97‑1 (MQ=255)
                                                                                                ccGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa  <  2:2087822/100‑1 (MQ=255)
                                                                                                         |                                                                                            
GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGA  >  NC_000913/3781053‑3781243

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: