Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,145	0	C	G	22.1%	213.9 / 24.6	86	intergenic (+128/‑70)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (40/27);  new base G (0/19);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 7.85e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACC  >  NC_000913/3781047‑3781240                                                                                                     |                                                                                                  cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCg                                                                                                  >  2:1978694/1‑100 (MQ=255) cccTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCg                                                                                                  >  2:611330/1‑100 (MQ=255)   cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGCGGGCCggg                                                                                                >  1:741342/1‑100 (MQ=255)    tgaGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCgggg                                                                                               <  2:228382/98‑1 (MQ=255)    tCAGGGAGAGGGGACCTTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCgggg                                                                                               <  1:2067165/100‑1 (MQ=255)    tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCgggg                                                                                               <  2:2963885/100‑1 (MQ=255)    tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCgggg                                                                                               <  1:1155033/100‑1 (MQ=255)        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                                                                                  <  1:3356413/100‑1 (MQ=255)               ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑aaaa                                                                                  <  2:1049736/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:617376/1‑100 (MQ=255)                                     cTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAAtttt                                                            >  2:2637044/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              >  1:488150/1‑100 (MQ=255)                                                                                    ccTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt              >  2:2333824/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)                                                                                                     |                                                                                                  CCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCC‑GGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACC  >  NC_000913/3781047‑3781240

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