Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,657,893	C→A	9.4%	intergenic (‑199/‑53)	tktA ← / → sahR	transketolase A/methionine metabolism transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,657,893	0	C	A	9.4%	187.8 / 7.2	61	intergenic (‑199/‑53)	tktA/sahR	transketolase A/methionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (22/33);  new base A (4/2);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

ATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTG  >  NC_002947/5657745‑5658038                                                                                                                                                       |                                                                                                                                                  ccaagcGCCCCCAGTGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAAc                                                                                                                                                   >  7:77776/6‑149 (MQ=255)      cGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                               <  4:318721/148‑1 (MQ=255)      cGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                               <  8:246039/148‑1 (MQ=255)       gCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTTg                                                                                                                                             <  3:349723/149‑1 (MQ=255)        cccACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGa                                                                                                                                            >  2:78273/1‑149 (MQ=255)           aCTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGatat                                                                                                                                         >  1:125074/1‑149 (MQ=255)               ggggCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTg                                                                                                                                       <  6:28257/147‑1 (MQ=255)               ggggCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGc                                                                                                                                     >  2:110213/1‑149 (MQ=255)               ggggCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAAATTAAGAAGCACAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTg                                                                                                                                       <  1:77058/148‑1 (MQ=255)                     cTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTCGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTgtt                                                                                                                               <  5:243576/149‑3 (MQ=255)                         gggCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAggg                                                                                                                           <  5:138643/149‑1 (MQ=255)                         gggCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAgg                                                                                                                            <  1:144663/148‑1 (MQ=255)                           gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCAGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAggg                                                                                                                           >  1:101426/1‑148 (MQ=255)                            cAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCAt                                                                                                                        >  8:24107/1‑149 (MQ=255)                                 aTTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGata                                                                                                                                          >  2:88739/1‑126 (MQ=255)                                 aTTTAAGCGCCAAGCGCCGGCTTTAAGACGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGata                                                                                                                                          <  1:88739/126‑1 (MQ=255)                                  tttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGc                                                                                                                  <  2:246161/149‑1 (MQ=255)                                   ttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTTGATATTg                                                                                                                                       >  5:163650/1‑127 (MQ=255)                                   ttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTTGATATTg                                                                                                                                       <  6:163650/127‑1 (MQ=255)                                        cgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCt                                                                                                            <  7:24107/149‑1 (MQ=255)                                        cgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCt                                                                                                            <  1:78273/149‑1 (MQ=255)                                         gccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTa                                                                                                           <  5:161808/149‑1 (MQ=255)                                         accaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                >  2:5505/2‑144 (MQ=255)                                         accaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                <  1:5505/143‑1 (MQ=255)                                           caagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGa                                                                                                         <  3:141603/149‑1 (MQ=255)                                           caagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGa                                                                                                         <  6:244074/149‑1 (MQ=255)                                            aagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCtt                                                                                                                                  <  5:234246/123‑1 (MQ=255)                                            aagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCtt                                                                                                                                  >  6:234246/1‑123 (MQ=255)                                            aagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGAc                                                                                                        <  6:38705/149‑1 (MQ=255)                                                gccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTgcc                                                                                                    >  4:84334/1‑149 (MQ=255)                                                           gCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGaa                                                                                          >  7:265115/1‑148 (MQ=255)                                                             cGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAACACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCt                                                                                       <  2:309296/149‑1 (MQ=255)                                                                  cTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGAt                                                                                                                    <  4:137226/115‑1 (MQ=255)                                                                  cTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGAt                                                                                                                    >  3:137226/1‑115 (MQ=255)                                                                     aagaagCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCGAAACGTTTTGATATTGGCGTTGCTAGGGCATCGATGCCTGTCTAGACTGGCGCCCCGTGCATCtgcgtggg                                                                                >  4:350251/1‑146 (MQ=255)                                                                       gaagCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCaa                                                                             >  3:328554/1‑149 (MQ=255)                                                                             gagCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATc                                                                       <  8:77776/149‑1 (MQ=255)                                                                             gagCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCAGCCCCATGAATCTGCGTGCGCAATCGATc                                                                       >  8:83908/1‑149 (MQ=255)                                                                              agCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcc                                                                      >  7:218577/1‑149 (MQ=255)                                                                                cgc‑gACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcccc                                                                    <  2:101426/149‑1 (MQ=255)                                                                                cgc‑gACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcccc                                                                    <  8:265115/149‑1 (MQ=255)                                                                                    gACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                 >  7:248553/1‑149 (MQ=255)                                                                                          gggCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcgc                                                           <  1:110213/149‑1 (MQ=255)                                                                                           ggCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcgcg                                                          <  4:175652/149‑1 (MQ=255)                                                                                              gcCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTg                                                                                                                                       <  4:227869/69‑1 (MQ=255)                                                                                              gcCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTg                                                                                                                                       >  3:227869/1‑69 (MQ=255)                                                                                                 ccAAAGCGTGTAACCTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACAt                                                    >  4:74247/1‑149 (MQ=255)                                                                                                  cAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcccc                                                                    <  4:116855/132‑1 (MQ=255)                                                                                                  cAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcccc                                                                    >  3:116855/1‑132 (MQ=255)                                                                                                   aaaGCGTGTAACCTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTg                                                  <  3:74247/149‑1 (MQ=255)                                                                                                      gCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCa                                               <  3:84334/149‑1 (MQ=255)                                                                                                             aCTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcgcg                                                          >  7:277056/1‑131 (MQ=255)                                                                                                             aCTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcgcg                                                          <  8:277056/131‑1 (MQ=255)                                                                                                                    g‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGc                                 <  2:125074/149‑1 (MQ=255)                                                                                                                                         ggCCAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGcc             >  4:206736/2‑149 (MQ=255)                                                                                                                                           ccAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt            >  1:254414/1‑148 (MQ=255)                                                                                                                                           ccAATATCAAAAAGTTTTGAAATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt            >  8:246813/1‑148 (MQ=255)                                                                                                                                              aaatCAAACCTTTTTGATATTGGCCTTGCTAGGGCCTCGTTGGCTGTCTAGACTGCCGGCCCTTGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAAc        <  2:332922/147‑1 (MQ=255)                                                                                                                                                  cAAAACTTTTTGATAATGCCCTTGCTAGGGCAACGATGCCTGTCTAGACAGCCGCCCCATGAATCTGCGTGCGCAATAGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAt    >  1:72742/1‑149 (MQ=255)                                                                                                                                                   aaaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAtt   <  8:218577/149‑1 (MQ=255)                                                                                                                                                    aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg  <  3:350251/149‑1 (MQ=255)                                                                                                                                                    aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg  >  2:311434/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  ATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGC‑GACTCGGGGCGCCCCAAAGCGTGTAACTTGAGG‑CTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTG  >  NC_002947/5657745‑5658038

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