Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	869,126	0	T	G	42.2%	107.1 / 24.7	66	V536G (GTG→GGG)	gsiA	glutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (5/32);  new base G (27/0);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGAT  >  NC_000913/868979‑869274                                                                                                                                                    |                                                                                                                                                                 ccGACGAAGCCGTTTCGGCGCTGGATGTTTCTTTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGAt                                                                                                                                                               <  2:427543/151‑1 (MQ=255)      gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑At‑c                                                                                                                                                        <  2:314008/151‑1 (MQ=255)      gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑At‑c                                                                                                                                                        <  2:378099/151‑1 (MQ=255)       aaGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑At‑ct                                                                                                                                                       <  1:1297340/151‑1 (MQ=255)        aGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑At‑ct‑c                                                                                                                                                     <  1:1369821/151‑1 (MQ=255)              ttCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGT‑AT‑TT‑C‑GGGcg                                                                                                                                               >  2:540693/1‑149 (MQ=255)              ttCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCa                                                                                                                                               <  1:746704/150‑1 (MQ=255)                cGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTTGCGGTGATTT‑TT‑CT‑C‑CGGCCAAt                                                                                                                                            >  1:316139/1‑151 (MQ=255)                cGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGCGGGGGGGT‑AT‑TT‑C‑gggggaa                                                                                                                                             >  2:242516/1‑146 (MQ=255)                 ggCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAAtt                                                                                                                                           <  2:812084/151‑1 (MQ=255)                    gcTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGGGGTGT‑AT‑TT‑C‑GGGGaaaatttt                                                                                                                                        >  1:107971/1‑146 (MQ=255)                       ggATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑Ga                                                                                                                                    <  1:700569/150‑1 (MQ=255)                         aTGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAAt                                                                                                                                 <  2:480845/151‑1 (MQ=255)                           gTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTTGTCATCGTGGGGCGGGGGTGG‑AT‑TTCC‑GGGGAAATTT‑TT‑GAAAtt                                                                                                                                >  1:212260/1‑151 (MQ=255)                             ttCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCCTCGTGTGGCGGGGGTGT‑AT‑TT‑C‑GGGGGAATTT‑TT‑TTAAATGGt                                                                                                                             >  2:205070/1‑151 (MQ=255)                               cTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGGGGGGGTGT‑AT‑CT‑CGGGGCAAATTT‑TT‑TAATTTgtgc                                                                                                                            >  1:766288/1‑148 (MQ=255)                                 aTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTcc                                                                                                                           <  1:1186386/149‑1 (MQ=255)                                 aTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTcc                                                                                                                           <  1:201102/149‑1 (MQ=255)                                 aTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTcc                                                                                                                           <  1:309802/149‑1 (MQ=255)                                 aTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGCGGGAAAGCGGGTTAGTTCTCTGGTTGCGGCGCTTT‑AT‑CT‑CTGGGCAAA‑TG‑TT‑GAAAATGGTCCAc                                                                                                                         >  1:1280488/1‑151 (MQ=255)                                    cgcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGgc                                                                                                                      <  2:1340664/151‑1 (MQ=255)                                    cgcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGGGGGGGTGT‑GT‑CT‑CGGGGGAAATTT‑TT‑AAATTTGGTCCCCgg                                                                                                                       >  1:334970/1‑151 (MQ=255)                                    cgcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGCGGGGGTGTGTT‑TT‑C‑GGGGAAATTT‑TT‑AAAATTGGTCCCCgc                                                                                                                       >  1:293406/1‑150 (MQ=255)                                       gggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTTGTCATCGTGTGGGGGGGGTGT‑TT‑TTCC‑GGGGAAATTT‑TG‑TAAATTgggcaccgcgcc                                                                                                                    >  2:1433244/1‑141 (MQ=255)                                       gggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGCGGGGGTGT‑AT‑AT‑C‑GGGGGAAATG‑TT‑TAAATTTGTCCACcgcgcg                                                                                                                   >  2:453107/1‑151 (MQ=255)                                         gCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGcgcgcg                                                                                                                 <  2:1356973/151‑1 (MQ=255)                                           aGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTTGTCATCGTGTGGCGGTGGTGT‑AT‑TTCC‑GGGCAAATT‑‑TT‑GTAAATTGTCCCCGGCGCGCggg                                                                                                               >  2:289142/1‑150 (MQ=255)                                                aTCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCg                                                                                                        <  1:762683/151‑1 (MQ=255)                                                aTCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGGGATTT‑CTCCC‑C‑GGGCAAATT‑‑TT‑GTAAATTGTGCCCGGCGCGCGG‑GC‑TTct                                                                                                        >  1:600093/1‑150 (MQ=255)                                                  cAACTTGTTTCTCCGTCTCCAGCGTTATTTCGGCATTGCGTATCTTTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGaa                                                                                                      <  1:1160603/151‑1 (MQ=255)                                                              cGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGAGTTGTCGTCGGGTGGGGGGGGTGG‑AT‑AT‑T‑GGGGAAATTG‑TTGGAAGTTGGGTAAC‑GCGCGCGG‑GG‑TTTGAAAAACCGCAGCAt                                                                                          >  2:107737/1‑151 (MQ=255)                                                              cGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAACGGGATAGTTAACGTCTGGGGGGGGTGG‑TT‑AA‑CGGGGCAAATCG‑TT‑‑AAATTCGTCCCCGGCCCGCGGCCC‑TTCCAAAAccgccaca                                                                                            >  1:1286941/1‑144 (MQ=255)                                                               gATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGGTAGTCATCGTGTGGCGGGGGTGT‑AT‑TT‑C‑CGGGGAAATG‑TT‑GTAAATTGTGCACCGCGCGCGG‑GG‑CTCGAAAAACCCCCGCATc                                                                                         >  2:536232/1‑151 (MQ=255)                                                                    ccAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCt                                                                                       <  2:1261288/148‑1 (MQ=255)                                                                      aGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGGCATCCTGTGGCGGGGGGCT‑AT‑CTCC‑CGGCAAATTT‑TT‑TAAATTGG‑CCACCGCCCGCGG‑GC‑TTCGAAAAACCCCCCCATCCCTTTAc                                                                                  >  2:1137100/1‑151 (MQ=255)                                                                           gATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTa                                                                            <  1:452801/151‑1 (MQ=255)                                                                                cGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTa                                                                       <  1:897195/151‑1 (MQ=255)                                                                                  gCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑Ct                                                                    <  1:1414889/151‑1 (MQ=255)                                                                                  gCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGT‑AT‑CT‑C‑GGGGGAAATT‑TT‑GAAATTTGGTCACCGCGCGCGG‑GG‑TTCGGAAAACCCCAGCAACCTTATACGC‑GGaaaattat                                                                     >  2:1112406/1‑146 (MQ=255)                                                                                        cGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGT‑AT‑CTCC‑GGGGAAATTT‑TT‑GAAATTGTTCC‑CGGCGCGCGG‑GC‑TTTGGAAAACCCCCGCATCCTTATACGC‑GGAAAATA‑CTtgc‑ggc                                                             >  2:79920/1‑151 (MQ=255)                                                                                        cGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGGGGGGGGGT‑AT‑TT‑C‑GGGCAAAATG‑TT‑GGAAATTGTGCACCGCGCGCGG‑GC‑TTCGGAAAACCCCAGCACTCTTATAAGC‑GGAAAATA‑CTtgc‑ggg                                                             >  1:151299/1‑150 (MQ=255)                                                                                            tCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGCGGGGGTGT‑AT‑CT‑C‑GGGCGAAATG‑TT‑TTAAATTGTCCCCGGCGCGCGG‑GC‑CTCTAAAAACCCCAGCACCCCTTTAAGC‑GCAAAATA‑TTGGG‑GGGA‑ggt                                                        >  2:360264/1‑149 (MQ=255)                                                                                             cTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGCGGGGATGT‑AT‑TT‑CGGGGAAAATTT‑TT‑AAAATTGG‑CCCCCGCGCGCGG‑GC‑TTCGAAAAACCCCAGCATCCCTTTAAGC‑GGAAAATT‑CTTGC‑GGCA‑ggtt                                                       >  1:34801/1‑148 (MQ=255)                                                                                               gTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATCCGC‑GTAAATTA‑CTGGC‑GGCA‑GTTcc                                                      <  2:1244524/150‑1 (MQ=255)                                                                                               gTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTcc                                                      <  1:723871/150‑1 (MQ=255)                                                                                               gTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGCGGGGATGG‑AT‑CTCC‑GGGCAATTGT‑TT‑GAAATTGG‑CCCCCGCGCGCGG‑GCTTTCGAAAACCC‑CCGCATCCCTTTACGC‑GGAAAATA‑ATTGG‑GGCA‑GGTccc                                                     >  1:1070272/1‑150 (MQ=255)                                                                                                 ttATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCgg                                                    <  1:485763/150‑1 (MQ=255)                                                                                                 ttATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGT‑TT‑TTCC‑GGGCAAATTT‑TT‑AAAATTGG‑CCCCCGCGCGCGG‑GC‑TTCGAAAAACCCCAGCAACCCTTTACGC‑GGAAAATA‑Ctttggggaggttccgg                                                     >  1:1340973/1‑136 (MQ=255)                                                                                                 ttATCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGCGGGGGTGT‑TT‑CT‑C‑GGGGAAAATGTTT‑TAAATTGG‑CCCCCGCGCGCGG‑GC‑TTTGAAAAACCCCCGCAACCCTTTACGC‑GGAAAATA‑CTTGC‑GGGA‑GTTcccgg                                                   >  1:566065/1‑148 (MQ=255)                                                                                                   aTCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGGGTGGCGGGGGAGT‑TT‑TTCC‑GGGGAAATTT‑TT‑AAAATT‑GTCCCCCGCGCGCGG‑GC‑TTTGAAAAACCCCAGGATCCTTTTAACC‑GGTAATTATTTGGC‑GGCG‑GTTCCggcg                                                  >  1:912908/1‑149 (MQ=255)                                                                                                   aTCTCCCACGATATGGCGGTGGGAGAGCGGATTAGTCATCGTGTGGCGGGGGAGT‑TT‑TT‑T‑GGGGGAAATG‑TT‑TTAAATTGTTCACCGGGCGCGG‑GC‑TTTCAAAAACCCCAACCTCCCTTTATCC‑GGAAAATA‑CTTGCGGGCG‑GTTCCggcg                                                  >  2:983252/1‑149 (MQ=255)                                                                                                    tctcCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGGGGTGT‑AT‑TT‑CGGGGCAAATT‑‑TT‑GTAAATTGTCCCCAGCGCGCGG‑GC‑TTTCAAAAACCCCCGCATCCTTTTAACC‑GGAAAATA‑CTTGC‑GGGAGGTTCCggcgc                                                 >  1:1277313/1‑148 (MQ=255)                                                                                                       cccACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCt                                               <  2:316143/149‑1 (MQ=255)                                                                                                       cccACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTg                                              <  1:1321967/150‑1 (MQ=255)                                                                                                            gATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGG‑AT‑TT‑C‑GGGGGAAATG‑TT‑TTAAATTGTTCACCGCGCGCGG‑GCTTTCAAAAACCCCCGCCACCCTT‑TACGCGGAAAATTC‑TTGGG‑GGCG‑GTTCCggggctggaccgc                                         >  2:86954/1‑140 (MQ=255)                                                                                                                 ggCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGATTTGT‑AT‑CG‑C‑AGGCAAATTG‑GT‑TAAATTGGTCCACGGCGGGCGG‑TC‑TTCCAAAACCCCCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑TTTCCTGTCCCTTAACCGTccctc                                   >  2:612585/1‑149 (MQ=255)                                                                                                                 ggCGGTGGTAGAGCGGATGAGTCGTCGTGTGGCAGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCTTCCCg                                    <  2:1305561/150‑1 (MQ=255)                                                                                                                      tggtAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACg                              >  2:354538/1‑151 (MQ=255)                                                                                                                          agagCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGAcc                           <  1:540694/150‑1 (MQ=255)                                                                                                                          agagCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCg                          >  1:904955/1‑151 (MQ=255)                                                                                                                          agagCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCg                          <  1:1037222/151‑1 (MQ=255)                                                                                                                                    aGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGGCCGCAGCGTGTa                 <  2:1050918/150‑1 (MQ=255)                                                                                                                                           gtgtGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTc         >  2:168489/1‑151 (MQ=255)                                                                                                                                                gCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACg    <  1:421402/151‑1 (MQ=255)                                                                                                                                                  ggTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGAt  <  1:156242/151‑1 (MQ=255)                                                                                                                                                  ggTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGAt  <  2:902749/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                                 CCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGT‑AT‑CT‑C‑GGGCAAATTG‑TT‑GAAATTGGTCCACGGCGCGCGG‑TC‑TTCGAAAACCCGCAGCATCCTTATACGC‑GTAAATTA‑CTGGC‑GGCA‑GTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAGCGTGTACTGCTGTCGGACGAT  >  NC_000913/868979‑869274

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑