Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1400AG24.1% 158.8 / 16.1 86intergenic (+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 (40/23);  new base G (0/20);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.34e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACG  >  NC_000913/3781041‑3781236
                                                                                                     |                                                                                                 
cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                                   >  2:798305/1‑100 (MQ=255)
 ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                  <  1:3043894/100‑1 (MQ=255)
 ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                  <  1:1127620/100‑1 (MQ=255)
 ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                  <  2:600219/100‑1 (MQ=255)
 ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                  >  2:1794515/1‑100 (MQ=255)
 ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                  >  2:2336675/1‑100 (MQ=255)
  ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                 >  1:969519/1‑100 (MQ=255)
  ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                 >  2:118298/1‑100 (MQ=255)
  ctctCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                 >  2:2805960/1‑100 (MQ=255)
   tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                                                <  2:2320669/100‑1 (MQ=255)
   tctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                                                <  2:1871277/100‑1 (MQ=255)
      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                             >  2:1978694/1‑100 (MQ=255)
      cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCg                                                                                             >  2:611330/1‑100 (MQ=255)
        cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCggg                                                                                           >  1:741342/1‑100 (MQ=255)
         tgaGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                                          <  2:228382/98‑1 (MQ=255)
         tCAGGGAGAGGGGACCTTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                                          <  1:2067165/100‑1 (MQ=255)
         tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                                          <  2:2963885/100‑1 (MQ=255)
         tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                                                          <  1:1155033/100‑1 (MQ=255)
             ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                                                                      >  1:2229574/1‑99 (MQ=255)
              gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                    <  2:664225/100‑1 (MQ=255)
              gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                    >  1:3407573/1‑100 (MQ=255)
               agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                                                                   <  1:2087250/100‑1 (MQ=255)
                gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgggg                                                                                  <  2:2799371/100‑1 (MQ=255)
                    ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                              <  2:1049736/100‑1 (MQ=255)
                    ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                              <  1:3356413/100‑1 (MQ=255)
                     gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                             <  2:898882/100‑1 (MQ=255)
                     gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                             <  1:1314495/100‑1 (MQ=255)
                      aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                            <  1:2687343/100‑1 (MQ=255)
                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                                                         >  2:873958/1‑100 (MQ=255)
                         gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                                                         >  2:37306/1‑100 (MQ=255)
                            cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCAc                                                                      >  1:1463864/1‑100 (MQ=255)
                                gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                  >  1:3085973/1‑100 (MQ=255)
                                gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCt                                                                  >  2:2945583/1‑100 (MQ=255)
                                 ntgtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                 >  1:258363/2‑100 (MQ=255)
                                 cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                 >  1:2561186/1‑100 (MQ=255)
                                   gtatgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                               >  2:2162693/1‑100 (MQ=255)
                                    tttgtcCTCCCTCCCCCTGAACCGCCCCCGCCCTTTTTTTTGCCCCCCCGCCCCCCGGGGGAGGGGGCCGGGGGGGGGGGAAAAGGCCGCACTGCTCAcc                                                              <  1:798301/98‑1 (MQ=255)
                                     atgtaCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa                                                             >  2:1646512/1‑100 (MQ=255)
                                         aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt                                                         <  1:2885342/100‑1 (MQ=255)
                                          ccccctccctCTGAACCGCCACCGCCCTGTTTTTTTCCCCCCCGCCCCCCCGGGGGGGGGCCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCaa                                                            <  1:3062809/94‑1 (MQ=255)
                                          ccccccCACTGGGAAAGGCCCCCGCTCTTTTTTTTTCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCAc                                                               <  2:52617/91‑1 (MQ=255)
                                          cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                        >  1:3275173/1‑100 (MQ=255)
                                          cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                        >  2:2637044/1‑100 (MQ=255)
                                          cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                        >  2:617376/1‑100 (MQ=255)
                                            ccctccctCTGAAACGCCCCCGCCCTCTTTTTTTCCCCCTCGCCCCCCGGGGGAGAGGGCGGGGTTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                      <  1:2336658/100‑1 (MQ=255)
                                            cccTCCTTCTGAAACGACACCGCACTCTTTTTTGCCCCCTCGCCCCCCGGGGGAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                       <  1:2364840/99‑1 (MQ=255)
                                            cccTCACCCCGAACCGCCCCCGCACTCTTTTTTTCCCCCGCGCCCCCCGGGGGGGGGGGCCGGGGGGGGGGGAAAAGGCCGCACTGCTCAc                                                               <  1:2460721/91‑1 (MQ=255)
                                             cccccctctAAAACGACCCCCCACTCTTTTTTCCCCCCCCGCCCCCCCGGGGGGAGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                     <  1:201967/95‑1 (MQ=255)
                                              cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                    >  1:2028801/1‑100 (MQ=255)
                                              cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                    >  1:975399/1‑100 (MQ=255)
                                                   ttagAAAAGCCACCGCACCATTTTTTTCTCCTTAGCCCCAAAGGGGAAAGGGCCGGGTAGAGGGAAAAAGGCAGCACTGCTCACCAATTTTCTGCTAt                                                 <  1:2950959/95‑1 (MQ=255)
                                                   tctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                               >  2:1614579/1‑100 (MQ=255)
                                                    ctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                              >  2:2335781/1‑100 (MQ=255)
                                                    cggAAACGCCCCCCCCCTTTTTTTTTCTCCCCCGCCCCGCGGGGGGGGGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                               <  2:1162239/97‑1 (MQ=255)
                                                     tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                             >  2:49730/1‑100 (MQ=255)
                                                       aaccGACCCCGCACTCTTTTTTTCCCCCTCCCCCCTCGGGGGAGAGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                             <  1:2022987/95‑1 (MQ=37)
                                                           gACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                       >  1:3246/1‑100 (MQ=255)
                                                             ccccGCATTCTTTTTTTCTCCCTCGCCCCCCGGGGGAGGGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                       <  1:971729/96‑1 (MQ=255)
                                                              ccccacctttttttttCCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                        <  2:671936/88‑1 (MQ=255)
                                                              acCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGcgc                                    >  1:2151895/1‑98 (MQ=255)
                                                               cccccctttttttttCCCCCCCCCCCCCCCGGGGAGGGGGCCGGTGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                    <  2:2541793/91‑1 (MQ=255)
                                                                ccctctctTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGCCGGGGGTTGGGGGAAAAGGCCGCACTGCTCACCAAGTTTCTGCTATCCTGccccc                                        <  1:2825985/90‑1 (MQ=255)
                                                                cGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTa                                  >  1:2837773/1‑100 (MQ=255)
                                                                 gCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTaa                                 >  2:2252348/1‑100 (MQ=255)
                                                                    ctctTTTTTTCTCTCCCGGGCCCTCCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAg                                <  2:975403/100‑1 (MQ=255)
                                                                    ctctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                              >  2:1989083/1‑100 (MQ=255)
                                                                    ctctTTTTTTCTCCCTCGCCCCCCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                              >  2:508962/1‑100 (MQ=255)
                                                                      ttttttttCCCCTCCGCCCCCCGGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                             <  2:3073719/98‑1 (MQ=255)
                                                                      ctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                             >  1:1915363/1‑99 (MQ=255)
                                                                      ctTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                            >  1:1386193/1‑100 (MQ=255)
                                                                       ttttttcCCCCCCCCCCCCCCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                           <  2:969523/100‑1 (MQ=255)
                                                                         tttttCCCCGTCGCCCCTGGGGGATGGGGCCCGGGATGGGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCTTAGTGCGGGGGCAg                         <  2:1762945/100‑1 (MQ=255)
                                                                         tttttCCCCCTCCCCCCCCGGGGGAGGGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                         <  2:297941/100‑1 (MQ=255)
                                                                         ttttcCTTCCTCGGCCCTGCCGGGGGGGGGCCGGGGCTGTGGCGCAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                           <  2:84898/99‑1 (MQ=255)
                                                                         ttttcCCCCCCCGCCCCGCCGGGGGGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                         <  1:2636135/100‑1 (MQ=255)
                                                                            ttCTCCCCCCCCCCTCGGGGGGGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                      <  2:1463870/100‑1 (MQ=255)
                                                                              ctcCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                    >  2:28520/1‑100 (MQ=255)
                                                                              cgcccGCGCACCTCCGGGGAGTGGGCGGAGGAAAGGAGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCACACTAAGGGGGCAGTATg                     <  1:873954/97‑1 (MQ=255)
                                                                                cccTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                  >  1:2324033/1‑100 (MQ=255)
                                                                                cccTCGCCCCCCGGGGGGGAGGGCGGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                   <  1:1767331/99‑1 (MQ=255)
                                                                                   cTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                >  2:507219/1‑100 (MQ=255)
                                                                                     cGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACa              >  2:3492935/1‑100 (MQ=255)
                                                                                       ccccccGGGGGGGAGGGCCGGGTGGAGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                <  1:129342/96‑1 (MQ=255)
                                                                                        cccTCCGGGGGGGGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg           <  1:1970135/100‑1 (MQ=255)
                                                                                         ccTCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGAGGCAGTATGCTAAACATCGt          <  1:2162679/100‑1 (MQ=255)
                                                                                         ccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt          >  2:2333824/1‑100 (MQ=255)
                                                                                         ccTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt          >  1:488150/1‑100 (MQ=255)
                                                                                          cccggggggggggggCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt          <  2:3407599/95‑1 (MQ=255)
                                                                                          cTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa         >  1:93931/1‑100 (MQ=255)
                                                                                           tCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc        >  1:1597829/1‑100 (MQ=255)
                                                                                            ccGGGGGGAGGGCCGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa         <  1:1646506/98‑1 (MQ=255)
                                                                                             cgggggAGGGGCCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt       <  2:1711906/98‑1 (MQ=255)
                                                                                               ggggggAGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt     <  1:399557/100‑1 (MQ=255)
                                                                                                  gagaGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg  >  2:801571/1‑100 (MQ=255)
                                                                                                     |                                                                                                 
CCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACG  >  NC_000913/3781041‑3781236

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: