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

GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGA  >  NC_000913/3781053‑3781243
                                                                                                         |                                                                                            
gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑Ga                                                                                              >  2:1964002/1‑100 (MQ=255)
 ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAg                                                                                             <  1:93270/100‑1 (MQ=255)
  gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAgg                                                                                            <  1:2612107/100‑1 (MQ=255)
    gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAgggg                                                                                          <  1:2778207/100‑1 (MQ=255)
     agGGGACCGTTCGGCGCTGTATGCACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑a                                                                                        <  1:856368/100‑1 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑aa                                                                                       >  1:1203314/1‑100 (MQ=255)
      ggggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑aa                                                                                       >  1:849270/1‑100 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑TCTCCTTG‑CCCCTCC‑CGGG‑TGAGGGCCGGGGG‑GGGGGG‑gaaa                                                                                     >  1:2255885/1‑100 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑aaaa                                                                                     >  1:2017473/1‑100 (MQ=255)
        ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑aaaa                                                                                     <  2:1886315/100‑1 (MQ=255)
         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAg                                                                                    >  1:192169/1‑100 (MQ=255)
          aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAgg                                                                                   >  2:321491/1‑100 (MQ=255)
          aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAgg                                                                                   <  1:1964020/100‑1 (MQ=255)
          aCCGTTCGGCGCTGTATGTACTCCCCCCCCCTGAAACGGCACCGCACTC‑TTTTTGT‑CTCCCTAG‑ACCGACC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAgg                                                                                   <  2:2866474/100‑1 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGc                                                                                  >  1:1663671/1‑100 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGc                                                                                  <  1:2060839/100‑1 (MQ=255)
           ccGTTCGGCGCTGTATGTACTCCCCCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGc                                                                                  <  1:2033296/100‑1 (MQ=255)
             gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCg                                                                                >  2:2805295/1‑100 (MQ=255)
                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCa                                                                              >  2:2871980/1‑100 (MQ=255)
                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCAc                                                                             >  1:888123/1‑100 (MQ=255)
                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCAc                                                                             >  2:1251389/1‑100 (MQ=255)
                 ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACt                                                                            >  2:2233330/1‑100 (MQ=255)
                 ggCGCTGTATGTACCCCCTCACTCGGAAACGACACCGCACTC‑TTTTTTC‑CCCCCCCC‑CCCCCCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCAc                                                                             <  2:180428/99‑1 (MQ=255)
                 ggCGCGGAATGTACTCCCTCACGCTGAACCGACCCCGCACTT‑TTTTTTT‑CTCCCGCG‑CCCCTCC‑GGGG‑GGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCAc                                                                             <  1:1076511/99‑1 (MQ=37)
                   cgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGc                                                                          >  1:74421/1‑100 (MQ=255)
                       gtatgtaCTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCAc                                                                      >  1:295664/1‑100 (MQ=255)
                         atgtaCTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCa                                                                    >  2:2150303/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAAt                                                                  >  1:1790276/1‑100 (MQ=255)
                           gtaCTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTACTCACCAAt                                                                  >  2:501496/1‑100 (MQ=255)
                              cTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAAtttt                                                               >  1:1726623/1‑100 (MQ=255)
                               tCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTc                                                              >  1:827612/1‑100 (MQ=255)
                               cccccccctctAAAACAACCCCGCAC‑C‑TTTTTTT‑CCACCCCGCCCCCTCG‑GGGG‑AGGGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTc                                                              <  1:30568/93‑1 (MQ=255)
                                ccctccctCTGAACCGACCCCGCCCTC‑TTTTTTT‑CCCCCCCG‑CCCC‑CC‑GGGG‑GGGGGGCCCGGGTGGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTc                                                              <  2:1760606/99‑1 (MQ=255)
                                ccctccctCTAAACCGCCCCCGC‑CTT‑TTTTTTT‑CCCCCCCC‑CCCCCCC‑GGGGAAGGGGGCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                             <  1:1819183/100‑1 (MQ=255)
                                ccctccctCTAAAACGCCCCCCCCTTT‑TTTTTTT‑CCCCCCCC‑CCCCTCG‑GGGG‑AGGGGGCCGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTc                                                              <  1:1953591/99‑1 (MQ=255)
                                 cccccctctGAAACGCCACGGCACTT‑TTTTTTG‑CCCCCTCC‑CCCCTCG‑GGGG‑AGAGGCCGGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                            <  2:919671/95‑1 (MQ=25)
                                  cTCCCTCGGAAACGACCCCGCCCTC‑TTTTTGT‑CTCCCTCC‑CCCCTCG‑GGGG‑GGGGGGCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           <  2:1605145/100‑1 (MQ=255)
                                  cTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGCGGGCCGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           >  1:1368371/1‑100 (MQ=255)
                                    ccctctGAAACGACCCCGCATTC‑TTTTTTT‑CCCCCTCC‑CCCCTCG‑GGGG‑GGGGGGCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTa                                                         <  2:1198879/98‑1 (MQ=255)
                                    cACTGTGAAACGCCCCCGCACTT‑TTTTTTT‑CCCCCTCG‑CCCCCCC‑GGGG‑AGAGGG‑CGGGGT‑GGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           <  2:187608/98‑1 (MQ=255)
                                      ctgcaaaaCGCCCCCCCCCTT‑TTTTTTC‑CCCCCCCG‑CCCCTCC‑GGGG‑AGGGGGCCGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                           <  1:1449570/91‑1 (MQ=255)
                                      ctctGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CACCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                       >  2:1627567/1‑100 (MQ=255)
                                        ctGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                     >  2:1773972/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTC‑TTTTTTT‑TTTCCTTG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTACTCACCAATTTTCTGCTATCCt                                                     >  2:924481/1‑99 (MQ=38)
                                         tGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    >  1:2455347/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    >  1:740311/1‑100 (MQ=255)
                                         tGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                    >  2:666317/1‑100 (MQ=255)
                                          gAAACGACACCGCTCTC‑TTTTTTT‑CTCCCTCA‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                   >  1:2043032/1‑100 (MQ=255)
                                           aaaCGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                  >  2:1949050/1‑100 (MQ=255)
                                           aaaCGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                  >  2:2183996/1‑100 (MQ=255)
                                            aaCGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                 >  2:1265944/1‑100 (MQ=255)
                                            aaCGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                 >  2:1445173/1‑100 (MQ=255)
                                              cGACCCCCCTCTT‑TTTTTTT‑CCCCCCCC‑CCCCTCG‑GGGG‑AGGGGGCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                <  1:1949068/99‑1 (MQ=255)
                                                cccccGCCCTT‑TTTTTTT‑CCCCCCCC‑CCCCTCG‑GGGG‑GGGGGGCGGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                               <  1:317129/95‑1 (MQ=255)
                                                 ccccGCTCTC‑TTTTTTTCCCCCCCCC‑CCCCCCC‑GGGG‑GGAGGCCGGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                <  1:2872006/95‑1 (MQ=255)
                                                 ccccGCACTC‑TTTTTGT‑CTCCCCCT‑CCCCCCC‑GGGG‑AGCGGCCTGGGTT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                             <  2:1304575/97‑1 (MQ=255)
                                                   ccGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                           >  1:1733623/1‑99 (MQ=255)
                                                    ccgcctc‑tTTTTTT‑CCCCCCCC‑CCCCCCCGGGGG‑GGAGGCCCGGGGG‑GGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                           <  2:74420/95‑1 (MQ=255)
                                                      ccttc‑tTTTTTC‑CCCCCCCC‑CCCCCCC‑GGGG‑GGAGGG‑CGGGGT‑GGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                           <  2:1517357/93‑1 (MQ=255)
                                                      ccctc‑tTTTTTC‑CGCCCTCC‑CCCCGCG‑GGGG‑GGGGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                          <  2:629626/95‑1 (MQ=255)
                                                        ctc‑tTTTTTT‑CTTTCTTG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                     >  2:68267/1‑100 (MQ=255)
                                                         cttttttttt‑CTCCCTCC‑CCCCGCG‑GGGG‑AGAGGGCCGGGTT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                              <  2:192166/88‑1 (MQ=38)
                                                          c‑tTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   >  2:2192909/1‑100 (MQ=255)
                                                          c‑tTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   >  1:2903263/1‑100 (MQ=255)
                                                            ttttttt‑CCCCCTCG‑CCCCCCC‑GGGG‑GGAGGCCCGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                   <  2:695845/99‑1 (MQ=255)
                                                            ttttttt‑CCCCCCCG‑CCCCTCG‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                  <  2:2115140/100‑1 (MQ=255)
                                                            ttttttt‑CCCCCCCC‑CCCCCCCGGGGG‑GGGGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                     <  1:2184015/98‑1 (MQ=255)
                                                              ttttt‑CCCCCCCG‑CCCCCCG‑GGGG‑AGAGGGCGGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                <  2:952736/100‑1 (MQ=255)
                                                                     tcCGTCG‑CCCCTCC‑GGGG‑GGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                          <  2:1603293/100‑1 (MQ=255)
                                                                      ccccCG‑CCCCCCG‑GGGG‑AGGGGCCCGGGGG‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                          <  1:2740521/99‑1 (MQ=255)
                                                                      cccTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  1:1912668/1‑100 (MQ=255)
                                                                      cccTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  1:638415/1‑100 (MQ=255)
                                                                      cccTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  2:1408672/1‑100 (MQ=255)
                                                                      cccTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                         >  2:757474/1‑100 (MQ=255)
                                                                       ccgcG‑CCCCTCG‑GGGG‑GGAGGGCCGGGGT‑GGTGGG‑AATAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        <  2:2043014/97‑1 (MQ=255)
                                                                       ccTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        >  1:2182433/1‑100 (MQ=255)
                                                                       ccTCC‑CCCCCCG‑GGGG‑GGGGGGCGGGGGG‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        <  1:1673293/100‑1 (MQ=255)
                                                                        cccccccccccC‑GGGG‑GGGGGGGCGGGGG‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                        <  1:405354/95‑1 (MQ=255)
                                                                        cTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                       >  1:1718698/1‑100 (MQ=255)
                                                                        cTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                       >  1:193005/1‑100 (MQ=255)
                                                                             ccccTCC‑GGGG‑AGAGGGCCGGGGT‑GGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                      <  2:875511/97‑1 (MQ=255)
                                                                             ccccTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                   >  1:541489/1‑100 (MQ=255)
                                                                             ccccTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                   >  1:240891/1‑100 (MQ=255)
                                                                               ccTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                 >  1:1279061/1‑100 (MQ=255)
                                                                                 cccggggg‑AGGGGGCGGGGGA‑GGGAGT‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                 <  1:1126574/95‑1 (MQ=255)
                                                                                  cc‑GGGG‑GGGGGGCGGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc               <  2:453081/99‑1 (MQ=255)
                                                                                  cc‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt              >  1:2063344/1‑100 (MQ=255)
                                                                                  cc‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt              >  2:1545874/1‑100 (MQ=255)
                                                                                      ggg‑AGGGGGCCGGGGT‑GAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa           <  2:1974232/100‑1 (MQ=255)
                                                                                       gg‑AGAGGGCCGGGGC‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc          <  2:2063325/100‑1 (MQ=255)
                                                                                        g‑agaGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg         >  2:2634682/1‑100 (MQ=255)
                                                                                           gagGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa       <  2:125697/100‑1 (MQ=255)
                                                                                             gggCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa        <  2:1203303/97‑1 (MQ=255)
                                                                                             gggCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc     >  1:2175468/1‑100 (MQ=255)
                                                                                             ggccgggggg‑gAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc          <  2:629571/90‑1 (MQ=255)
                                                                                               gcgggggg‑gggggA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg   <  1:2233350/97‑1 (MQ=255)
                                                                                                ccGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa  <  2:2087822/100‑1 (MQ=255)
                                                                                                         |                                                                                            
GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTT‑CTCCCTCG‑CCCCTCC‑GGGG‑AGAGGGCCGGGGT‑GAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGA  >  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: ‑