Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,272,965	C→G	6.2%	L38V (CTA→GTA)	rlmG →	ribosomal RNA large subunit methyltransferase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,272,965	0	C	G	6.2%	203.2 / 4.0	65	L38V (CTA→GTA)	rlmG	ribosomal RNA large subunit methyltransferase G
Reads supporting (aligned to +/- strand):  ref base C (27/33);  new base G (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GCAACGGGCAGTTTCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGG  >  NC_002947/5272821‑5273112                                                                                                                                                  |                                                                                                                                                        gCAACGGGCAGTTTCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTa‑‑‑ca                                                                                                                                                  <  3:228647/149‑1 (MQ=255) gCAACGGGCAGTTTCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTa‑‑‑ca                                                                                                                                                  <  7:226583/149‑1 (MQ=255)     acggCAGTTTCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGAt                                                                                                                                              <  4:275428/147‑1 (MQ=255)       ggCAGTTTCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCa                                                                                                                                            <  7:256919/149‑1 (MQ=255)        gCAGTTTCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCTACGACCCGTTGGAGGCCTTCGATGCCG‑CCGACGTGTACCTGCTTGCTCAGCTA‑‑‑GAGGATCAg                                                                                                                                           <  8:42104/149‑1 (MQ=255)            tttCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGACGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCAc                                                                                                                                       >  3:117239/1‑149 (MQ=255)                   cccGAGCCAGTCGCCATGCCACTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCAATCCGTAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CAGACGAGTACCTGCTTGCGCAGATA‑‑‑CACGAGCATTCACCCGACGc                                                                                                                                <  6:26824/149‑1 (MQ=255)                   cccGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCTAAGTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGc                                                                                                                                <  8:169241/149‑1 (MQ=255)                   cccGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGc                                                                                                                                <  3:192837/149‑1 (MQ=255)                    ccGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCGTCCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGc                                                                                                                                <  1:238723/149‑1 (MQ=255)                       aGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAAt                                                                                                                            <  8:31606/149‑1 (MQ=255)                        gCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAAtt                                                                                                                           <  6:120731/149‑1 (MQ=255)                              cGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCggg                                                                                                                     <  6:218595/149‑1 (MQ=255)                                ccATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGcc                                                                                                                        >  1:321056/1‑144 (MQ=255)                                ccATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGcc                                                                                                                        <  2:321056/144‑1 (MQ=255)                                 cATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTAc                                                                                                                  >  2:17786/1‑149 (MQ=255)                                   tGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTg                                                                                                                <  5:21684/149‑1 (MQ=255)                                     ccccTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGt                                                                                                              >  7:225972/1‑149 (MQ=255)                                          gTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTTCTGGTGCTc                                                                                                          >  7:67578/1‑148 (MQ=255)                                           ttGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCaa                                                                                                        <  3:152647/149‑1 (MQ=255)                                                   ccccTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCt                                                                                                >  8:98247/1‑149 (MQ=255)                                                   ccccTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCt                                                                                                >  4:230306/1‑149 (MQ=255)                                                    cccTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTCCTTGCGCAGCTA‑‑‑CACGATCat                                                                                                                                           >  6:280319/1‑104 (MQ=255)                                                      cTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCg                                                                                             >  5:257429/1‑149 (MQ=255)                                                         cGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTTCTGGTGCTCAATGACAGCTTCGGt                                                                                           <  8:67578/148‑1 (MQ=255)                                                            agAACTCGACCTGATCCGCCAGCCGGAGCAGGCAAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCATGCACCCGACGCCAATTGACGGCTACTGGTGCTCAATGACAGCTAAGGTGacc                                                                                       >  6:9467/2‑149 (MQ=255)                                                                 tCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGcc                                                                                  >  3:265225/1‑149 (MQ=255)                                                                   gACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑g                                                                                >  3:159549/1‑148 (MQ=255)                                                                    aCCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑Gc                                                                               >  2:192177/1‑148 (MQ=255)                                                                     ccTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCaa                                                                             <  3:232216/149‑1 (MQ=255)                                                                        gATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTAGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGcc                                                                          >  1:135409/1‑149 (MQ=255)                                                                         aTCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGATTACCTGCTTGCGCAGCTACACCACAATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCaa                                                                             <  1:267330/148‑1 (MQ=255)                                                                            cGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTg                                                                       <  3:230306/148‑1 (MQ=255)                                                                             gCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCa                                                                     <  2:80122/149‑1 (MQ=255)                                                                               cAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAgg                                                                   >  1:184385/1‑149 (MQ=255)                                                                                aGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGcc                                                                                  <  8:15131/134‑1 (MQ=255)                                                                                aGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGcc                                                                                  >  7:15131/1‑134 (MQ=255)                                                                                     gAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCCTGCAAGGCTTGC‑GGCCAggg                                                             <  4:117239/149‑3 (MQ=255)                                                                                            ccAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTAGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAggtg                                                      <  2:135409/149‑1 (MQ=255)                                                                                             cAACGACCCGTTGTAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAggtgg                                                     <  2:141604/149‑1 (MQ=255)                                                                                                cGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGCTTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCGCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTca                                                  <  5:9467/149‑1 (MQ=255)                                                                                                cGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTca                                                  <  7:98247/149‑1 (MQ=255)                                                                                                    ccGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTcagcag                                              >  8:71265/1‑149 (MQ=255)                                                                                                    ccGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGAACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGCCAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGCTGCTcagcag                                              >  3:118833/1‑149 (MQ=255)                                                                                                      gTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAg                                                                    >  6:12638/1‑125 (MQ=255)                                                                                                      gTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAg                                                                    <  5:12638/125‑1 (MQ=255)                                                                                                         gCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCTGGCCAGCTGCAGGTGGTCAGCAGCGGt                                          <  6:257429/148‑1 (MQ=255)                                                                                                         gCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGt                                          <  6:150674/148‑1 (MQ=255)                                                                                                           aGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGAc                                       <  8:225972/149‑1 (MQ=255)                                                                                                           aGGCCTTCGATGACG‑ACGACGAGTACCTGCTTGCGCAGCTA‑‑‑AAAGATCAGGCAACCGACTCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCAGTGAc                                       <  4:265225/149‑1 (MQ=255)                                                                                                             gCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTc                                     >  8:163401/1‑149 (MQ=255)                                                                                                               cTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGc                                   <  3:273658/149‑1 (MQ=255)                                                                                                                ttCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGGTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCa                                                          >  3:224961/1‑125 (MQ=255)                                                                                                                ttCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGGTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCa                                                          <  4:224961/125‑1 (MQ=255)                                                                                                                       ccg‑ccgACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCa                                                                     <  7:168556/107‑1 (MQ=255)                                                                                                                       ccg‑ccgACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCa                                                                     >  8:168556/1‑107 (MQ=255)                                                                                                                           ccgACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCAc                        <  7:163401/149‑1 (MQ=255)                                                                                                                                gaGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGc                   >  2:12949/1‑149 (MQ=255)                                                                                                                                     ccTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGCCAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGGAGGTGGTTAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTg               >  4:176417/1‑148 (MQ=255)                                                                                                                                      cTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTgg              >  2:57559/1‑148 (MQ=255)                                                                                                                                       tGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGa             >  2:28117/1‑148 (MQ=255)                                                                                                                                        gcgcaagcAGGTA‑‑‑CTCGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTgggg                           >  6:160834/7‑133 (MQ=255)                                                                                                                                        gcgcaagcAGGTA‑‑‑CTCGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTgggg                           <  5:160834/127‑1 (MQ=255)                                                                                                                                          ttGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGaa            >  6:153114/1‑146 (MQ=255)                                                                                                                                                 gCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTgg  >  6:145200/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                        GCAACGGGCAGTTTCGTACCCGAGCCAGACGCCATGCCCCTGTTGACCACCCCCTACGCCGAACTCGACCTGATCCGCCAGCCGGAGCAGGCCAACGACCCGTTGCAGGCCTTCGATGCCG‑CCGACGAGTACCTGCTTGCGCAGCTA‑‑‑CACGATCAGGCACCCGACGCCAATTGCCGGGTACTGGTGCTCAATGACAGCTTCGGTGCCCTGGCC‑GCAAGCCTTGCAGGCCAGCTGCAGGTGGTCAGCAGCGGTGACTCGCACCTGGGGCACCTGGCCCTGGAAAAGAACCTGG  >  NC_002947/5272821‑5273112

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