Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	854,663	G→A	7.1%	A88T (GCC→ACC)	moeB →	molybdopterin synthase adenylyltransferase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	854,663	0	G	A	7.1%	209.4 / 4.0	70	A88T (GCC→ACC)	moeB	molybdopterin synthase adenylyltransferase subunit
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

CTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCGG  >  NC_002947/854519‑854807                                                                                                                                                  |                                                                                                                                                  cTGGTCTGGCCAGTCGCCCTTTACCTGGCCGCCGCTGTCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACTCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGTTCGACTC‑GGCCtt                                                                                                                                               <  6:259620/149‑1 (MQ=255) cTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCtt                                                                                                                                               <  4:60295/149‑1 (MQ=255)    ggCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCa                                                                                                                                            <  2:221915/149‑1 (MQ=255)         gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCCCGACCGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCt                                                                                                                                       >  4:206610/1‑149 (MQ=255)         gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAgcg                                                                                                                                         <  8:31382/147‑1 (MQ=255)         gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCt                                                                                                                                       >  6:158693/1‑149 (MQ=255)         gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTC‑GACCTTGctcat                                                                                                                                         <  2:124303/147‑5 (MQ=255)         gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTC‑GACCTTGctcat                                                                                                                                         <  1:98340/147‑5 (MQ=255)          ccAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCtt                                                                                                                                      >  3:64311/1‑149 (MQ=255)             gTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCTAGGTCATCCTCGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCa                                                                                                                                   <  8:265308/149‑1 (MQ=255)                gCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAgg                                                                                                                                 <  5:158693/148‑1 (MQ=255)                 cccTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGc                                                                                                                                >  3:188562/1‑148 (MQ=255)                  ccTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCaa                                                                                                                              >  2:67965/1‑149 (MQ=255)                      tACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTC‑GACCTTGctc                                                                                                                                           >  1:283864/1‑130 (MQ=255)                      tACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTC‑GACCTTGctc                                                                                                                                           <  2:283864/132‑3 (MQ=255)                         cTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AAccc                                                                                                                      >  1:241130/1‑149 (MQ=255)                               gccgcTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCACCCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑ACCCCGGAAGt                                                                                                                >  4:243577/1‑149 (MQ=255)                               gccgcTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACTACAGTGCCAGCGTTGGCATGAGCAAGTTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGTAAGt                                                                                                                <  6:354663/149‑1 (MQ=255)                                      gCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGcc                                                                                                           <  1:67965/147‑1 (MQ=255)                                        gTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGt                                                                                                       <  7:258838/149‑1 (MQ=255)                                        gTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGt                                                                                                       >  5:175131/1‑149 (MQ=255)                                              gAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCa                                                                                                 >  1:343443/1‑149 (MQ=255)                                                gCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCTGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCAcc                                                                                               >  3:50171/1‑149 (MQ=255)                                                gCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGcc                                                                                                           <  1:227936/137‑1 (MQ=255)                                                gCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGcc                                                                                                           >  2:227936/1‑137 (MQ=255)                                                 cTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCg                                                                                              >  8:185843/1‑149 (MQ=255)                                                    cACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGcc                                                                                            <  2:241130/148‑1 (MQ=255)                                                            cGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTg                                                                                    >  3:352112/1‑148 (MQ=255)                                                             gACTTCGACACCTTCGACCTGACAAACCTGCAACGCCAGTTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCTGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGa                                                                                  <  7:214461/149‑1 (MQ=255)                                                              aCTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGAc                                                                                 <  1:349628/149‑1 (MQ=255)                                                              aCTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGAc                                                                                 <  3:243577/149‑1 (MQ=255)                                                                        cGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGCACTCGc                                                                       >  3:295390/1‑149 (MQ=255)                                                                         gTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCt                                                                      >  6:319641/1‑149 (MQ=255)                                                                            gACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATCTAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGa                                                                               <  1:76685/138‑1 (MQ=255)                                                                            gACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGa                                                                               >  2:76685/1‑137 (MQ=255)                                                                            gACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGATGACTCGCtggt                                                                   <  4:352112/149‑1 (MQ=255)                                                                             aCCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtggtg                                                                  <  3:268650/149‑1 (MQ=255)                                                                                  aCCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGc                                                             >  7:292416/1‑149 (MQ=255)                                                                                    cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtg                                                                     >  1:267926/1‑139 (MQ=255)                                                                                    cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtg                                                                     <  2:267926/139‑1 (MQ=255)                                                                                    cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCgg                                                           <  3:339436/149‑1 (MQ=255)                                                                                       ccTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGGTGCGGTgg                                                        >  4:303104/1‑149 (MQ=255)                                                                                           cAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTggcggc                                                    <  6:175131/149‑1 (MQ=255)                                                                                           cAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTggcgg                                                     <  3:206610/148‑1 (MQ=255)                                                                                            aaCGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAg                                                                              <  2:256802/122‑1 (MQ=255)                                                                                            aaCGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAg                                                                              >  1:256802/1‑122 (MQ=255)                                                                                            aaCGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGgcggcg                                                   <  5:319641/149‑1 (MQ=255)                                                                                                 cAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGa                                              >  3:20332/1‑149 (MQ=255)                                                                                                    gTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCt                                           >  7:63938/1‑149 (MQ=255)                                                                                                     tCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTg                                          >  5:220478/1‑149 (MQ=255)                                                                                                       aTCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGt                                        >  2:194319/1‑149 (MQ=255)                                                                                                          cACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCt                                     >  6:165072/1‑149 (MQ=255)                                                                                                           aCGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTg                                    <  4:965/149‑1 (MQ=255)                                                                                                             gACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGa                                  <  2:343443/149‑1 (MQ=255)                                                                                                              aCAGTGCCAGCGTGGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGAc                                 <  4:64311/149‑1 (MQ=255)                                                                                                                  tGCCAGCGTTGGCATGAGCAAGGTCGACTCGGGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGc                              >  1:34028/1‑149 (MQ=255)                                                                                                                      aGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGa                         >  2:342027/1‑149 (MQ=255)                                                                                                                      aGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGa                         <  1:342027/149‑1 (MQ=255)                                                                                                                      aGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGa                         >  2:300744/1‑149 (MQ=255)                                                                                                                          ttGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAAt                     <  8:63938/149‑1 (MQ=255)                                                                                                                             gCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTtc                  >  1:117334/1‑149 (MQ=255)                                                                                                                               aTGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTgg                                         >  4:350918/1‑124 (MQ=255)                                                                                                                               aTGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTgg                                         <  3:350918/124‑1 (MQ=255)                                                                                                                                  tgCAAGGCCGAGTC‑GACCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGTCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTAc             >  6:349624/2‑149 (MQ=255)                                                                                                                                   gCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTAcc            <  5:165072/149‑1 (MQ=255)                                                                                                                                       ggTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGCCCTGGTGCTGGACTGCTCCGACAATTTCTCTACCcgc         >  6:337674/1‑148 (MQ=255)                                                                                                                                          cGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGAc                                             >  4:139732/1‑109 (MQ=255)                                                                                                                                          cGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGAc                                             <  3:139732/109‑1 (MQ=255)                                                                                                                                           gACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAgg     <  4:20332/148‑1 (MQ=255)                                                                                                                                             cTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCgg  >  2:151550/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTC‑GGCCTTGCAGCGCTTGCAGGCAATC‑AACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGCGAGGCGG  >  NC_002947/854519‑854807

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