Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1520AG26.5% 129.2 / 21.5 84intergenic (+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 (37/24);  new base G (0/22);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 1.23e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGC  >  NC_000913/3781054‑3781251
                                                                                                     |                                                                                                    
ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGa                                                                                                      >  1:2229574/1‑99 (MQ=255)
 gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAgg                                                                                                    <  2:664225/100‑1 (MQ=255)
 gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAgg                                                                                                    >  1:3407573/1‑100 (MQ=255)
  agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAggg                                                                                                   <  1:2087250/100‑1 (MQ=255)
   gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAgggg                                                                                                  <  2:2799371/100‑1 (MQ=255)
       ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑aaaa                                                                                             <  2:1049736/100‑1 (MQ=255)
       ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑aaaa                                                                                             <  1:3356413/100‑1 (MQ=255)
        gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAg                                                                                            <  2:898882/100‑1 (MQ=255)
        gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAg                                                                                            <  1:1314495/100‑1 (MQ=255)
         aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAgg                                                                                           <  1:2687343/100‑1 (MQ=255)
            gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCg                                                                                        >  2:37306/1‑100 (MQ=255)
            gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCg                                                                                        >  2:873958/1‑100 (MQ=255)
               cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCAc                                                                                     >  1:1463864/1‑100 (MQ=255)
                   gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCt                                                                                 >  1:3085973/1‑100 (MQ=255)
                   gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCt                                                                                 >  2:2945583/1‑100 (MQ=255)
                    ntgtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTc                                                                                >  1:258363/2‑100 (MQ=255)
                    cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTc                                                                                >  1:2561186/1‑100 (MQ=255)
                      gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCAc                                                                              >  2:2162693/1‑100 (MQ=255)
                       tttgtcCTCCCTCCCCCTGAACCGCCCCCGCCCTTTTTTTTGCCC‑CCCCGCCCCCCG‑GGGGAGGGGGCCGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCAcc                                                                             <  1:798301/98‑1 (MQ=255)
                        atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCa                                                                            >  2:1646512/1‑100 (MQ=255)
                            aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGACCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAAttt                                                                        <  1:2885342/100‑1 (MQ=255)
                             ccccctccctCTGAACCGCCACCGCCCTGTTTTTTTCCC‑CCCCGCCCCCCC‑GGGGGGGGGCCCGGGG‑GGAGGGG‑AAAAGGCCGCACTGCTCACCaa                                                                           <  1:3062809/94‑1 (MQ=255)
                             ccccccCACTGGGAAAGGCCCCCGCTCTTTTTTTTTCCC‑CCCCCCCCCCCG‑GGGGGGGGGGCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCAc                                                                              <  2:52617/91‑1 (MQ=255)
                             cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAAtttt                                                                       >  2:2637044/1‑100 (MQ=255)
                             cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAAtttt                                                                       >  2:617376/1‑100 (MQ=255)
                             cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAAtttt                                                                       >  1:3275173/1‑100 (MQ=255)
                               ccctccctCTGAAACGCCCCCGCCCTCTTTTTTTCCC‑CCTCGCCCCCCG‑GGGGAGAGGGCGGGGT‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCt                                                                     <  1:2336658/100‑1 (MQ=255)
                               cccTCCTTCTGAAACGACACCGCACTCTTTTTTGCCC‑CCTCGCCCCCCG‑GGGGAGGGGGCGGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTc                                                                      <  1:2364840/99‑1 (MQ=255)
                               cccTCACCCCGAACCGCCCCCGCACTCTTTTTTTCCC‑CCGCGCCCCCCG‑GGGGGGGGGGCCGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCAc                                                                              <  1:2460721/91‑1 (MQ=255)
                                cccccctctAAAACGACCCCCCACTCTTTTTTCCCC‑CCCCGCCCCCCC‑GGGGGGAGGGCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTg                                                                    <  1:201967/95‑1 (MQ=255)
                                 cTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                   >  1:2028801/1‑100 (MQ=255)
                                 cTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGc                                                                   >  1:975399/1‑100 (MQ=255)
                                      ttagAAAAGCCACCGCACCATTTTTTTCTC‑CTTAGCCCCAAA‑GGGGAAAGGGCCGGGT‑AGAGGGA‑AAAAGGCAGCACTGCTCACCAATTTTCTGCTAt                                                                <  1:2950959/95‑1 (MQ=255)
                                      tctGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                              >  2:1614579/1‑100 (MQ=255)
                                       ctGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                             >  2:2335781/1‑100 (MQ=255)
                                       cggAAACGCCCCCCCCCTTTTTTTTTCTC‑CCCCGCCCCGCG‑GGGGGGGGGGCGGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                              <  2:1162239/97‑1 (MQ=255)
                                        tGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                            >  2:49730/1‑100 (MQ=255)
                                          aaccGACCCCGCACTCTTTTTTTCCC‑CCTCCCCCCTCG‑GGGGAGAGGGCCGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                            <  1:2022987/95‑1 (MQ=37)
                                              gACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                                      >  1:3246/1‑100 (MQ=255)
                                                ccccGCATTCTTTTTTTCTC‑CCTCGCCCCCCG‑GGGGAGGGGGCGGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                                      <  1:971729/96‑1 (MQ=255)
                                                 ccccacctttttttttCCC‑CCCCCCCCCCCCGGGGGGGGGGGCGGGGG‑GGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                       <  2:671936/88‑1 (MQ=255)
                                                 acCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGcgc                                                   >  1:2151895/1‑98 (MQ=255)
                                                  cccccctttttttttCCC‑CCCCCCCCCCCC‑GGGGAGGGGGCCGGTG‑GGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                                   <  2:2541793/91‑1 (MQ=255)
                                                   ccctctctTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGTGGCCGGGGG‑TTGGGGG‑AAAAGGCCGCACTGCTCACCAAGTTTCTGCTATCCTGccccc                                                       <  1:2825985/90‑1 (MQ=255)
                                                   cGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                                 >  1:2837773/1‑100 (MQ=255)
                                                    gCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                                >  2:2252348/1‑100 (MQ=255)
                                                       ctctTTTTTTCTCTCCCGGGCCCTCCGGGGGAGGGGGCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                               <  2:975403/100‑1 (MQ=255)
                                                       ctctTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                             >  2:1989083/1‑100 (MQ=255)
                                                       ctctTTTTTTCTC‑CCTCGCCCCCCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                             >  2:508962/1‑100 (MQ=255)
                                                         ttttttttCCC‑CTCCGCCCCCCG‑GGGGGGGGGGCGGGGG‑GGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                            <  2:3073719/98‑1 (MQ=255)
                                                         ctTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCCCTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                            >  1:1915363/1‑99 (MQ=255)
                                                         ctTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                           >  1:1386193/1‑100 (MQ=255)
                                                          ttttttcCCC‑CCCCCCCCCCCG‑GGGGGGGGGGCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                          <  2:969523/100‑1 (MQ=255)
                                                            tttttCCC‑CGTCGCCCCTGG‑GGGATGGGGCCCGGGA‑TGGGGGT‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCTTAGTGCGGGGGCAg                                        <  2:1762945/100‑1 (MQ=255)
                                                            tttttCCC‑CCTCCCCCCCCG‑GGGGAGGGGGCCGGGG‑GGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                        <  2:297941/100‑1 (MQ=255)
                                                            ttttcCTT‑CCTCGGCCCTGC‑CGGGGGGGGGCCGGGGCTGTGGCG‑CAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                          <  2:84898/99‑1 (MQ=255)
                                                            ttttcCCC‑CCCCGCCCCGCC‑GGGGGGGGGGCGGGGG‑GGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                        <  1:2636135/100‑1 (MQ=255)
                                                               ttCTC‑CCCCCCCCCTCG‑GGGGGGGGGGCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                     <  2:1463870/100‑1 (MQ=255)
                                                                 ctc‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                                   >  2:28520/1‑100 (MQ=255)
                                                                 cgc‑ccGCGCACCTCC‑GGGGAGTGGGCGGAGG‑AAAGGAG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCACACTAAGGGGGCAGTATg                                    <  1:873954/97‑1 (MQ=255)
                                                                   c‑ccTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                 >  1:2324033/1‑100 (MQ=255)
                                                                   c‑ccTCGCCCCCCG‑GGGGGGAGGGCGGGGG‑GGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                                  <  1:1767331/99‑1 (MQ=255)
                                                                      cTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                               >  2:507219/1‑100 (MQ=255)
                                                                        cGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa                             >  2:3492935/1‑100 (MQ=255)
                                                                          ccccccG‑GGGGGGAGGGCCGGGT‑GGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                               <  1:129342/96‑1 (MQ=255)
                                                                           cccTCC‑GGGGGGGGGGCCGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                          <  1:1970135/100‑1 (MQ=255)
                                                                            ccTCG‑GGGGAGGGGGCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGAGGCAGTATGCTAAACATCGt                         <  1:2162679/100‑1 (MQ=255)
                                                                            ccTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                         >  2:2333824/1‑100 (MQ=255)
                                                                            ccTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                         >  1:488150/1‑100 (MQ=255)
                                                                             cccggggggggggggCGGGGG‑GGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                         <  2:3407599/95‑1 (MQ=255)
                                                                             cTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                        >  1:93931/1‑100 (MQ=255)
                                                                              tCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                       >  1:1597829/1‑100 (MQ=255)
                                                                               cc‑GGGGGGAGGGCCGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                        <  1:1646506/98‑1 (MQ=255)
                                                                                cgggggAGGGGCCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                      <  2:1711906/98‑1 (MQ=255)
                                                                                  ggggggAGGGCGGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                    <  1:399557/100‑1 (MQ=255)
                                                                                     gagaGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg                 >  2:801571/1‑100 (MQ=255)
                                                                                     aagggggCCGGTG‑TAAGGGA‑AAAAGGCCGCACTGCTCAACAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                          <  1:3001616/87‑1 (MQ=255)
                                                                                         gggCCGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc             <  1:2945563/100‑1 (MQ=255)
                                                                                         gggCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc             >  1:118008/1‑100 (MQ=255)
                                                                                         ggcccGGGG‑GGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc             <  1:260398/97‑1 (MQ=255)
                                                                                          gggcGGGG‑TGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc             <  2:1111040/97‑1 (MQ=255)
                                                                                               ggg‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt       <  1:617373/100‑1 (MQ=255)
                                                                                               ggg‑TGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt       <  2:1260398/100‑1 (MQ=255)
                                                                                                    gAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTcc    <  2:2028814/99‑1 (MQ=255)
                                                                                                    gAGGGG‑AAAAGGCAGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCg   >  1:2193031/1‑100 (MQ=255)
                                                                                                     aGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGc  >  1:1211623/1‑100 (MQ=255)
                                                                                                     |                                                                                                    
GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGC  >  NC_000913/3781054‑3781251

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: ‑