Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	535,818	G→T	8.6%	E85* (GAA→TAA) ‡	rplA →	50S ribosomal protein L1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	535,818	0	G	T	8.6%	268.5 / 9.4	86	E85* (GAA→TAA) ‡	rplA	50S ribosomal protein L1
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base T (3/5);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CAAGTTCGTCGAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGG  >  NC_002947/535673‑535960                                                                                                                                                   |                                                                                                                                                 cAAGTTCGTCGAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGTCACTGGCAAGACCGTTCTCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAg                                                                                                                                               <  1:346754/149‑1 (MQ=255)  aaGTTCGTCGAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAgc                                                                                                                                              >  8:319513/1‑149 (MQ=255)   aGTTCGTCGAGTCCTACGACATCGCCGTTAAGCCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAgc                                                                                                                                              >  1:218620/1‑149 (MQ=255)    gTTCGTCGAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAgctg                                                                                                                                            >  2:37889/1‑149 (MQ=255)      tcgtagAGTACTATGACATCGACGTTAA‑CCTCGGTGTTGACCAGCGTAATTCCGACCAGGTCGTACGTAGCTCTACCGTGCTGCCACAAGGCACTGTCAATACCTTTCGCGTTGCTGTCTTCACCCAGTGTCCAGCAGCTGAAgctgc                                                                                                                                           <  7:288775/148‑1 (MQ=255)         tcgAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACGAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgctgga                                                                                                                                       <  5:48080/149‑4 (MQ=255)         tcgAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgctgga                                                                                                                                       <  3:89858/149‑4 (MQ=255)          cgAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑gg                                                                                                                                     >  2:79908/1‑149 (MQ=255)            aGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCt                                                                                                                                   <  4:250509/149‑1 (MQ=255)                cTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCg                                                                                                                               >  4:162307/1‑149 (MQ=255)                 tACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCgg                                                                                                                              >  5:222779/1‑149 (MQ=255)                       aTCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACGGTTCGCGTTGCTGTCTTCACCCAGGCTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGa                                                                                                                        >  8:117571/1‑149 (MQ=255)                       aTCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGa                                                                                                                        >  4:105476/1‑149 (MQ=255)                         cGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑c                                                                                                                     >  4:232349/1‑149 (MQ=255)                         cGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑c                                                                                                                     >  3:128246/1‑149 (MQ=255)                           ccGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑cgc                                                                                                                   >  2:5462/1‑149 (MQ=255)                           ccGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑cgc                                                                                                                   <  6:222779/149‑1 (MQ=255)                              ttAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑cgcg                                                                                                                  >  7:281209/1‑147 (MQ=255)                              ttAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTa                                                                                                                <  2:218620/149‑1 (MQ=255)                                 a‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGt                                                                                                             <  8:281209/149‑1 (MQ=255)                                 a‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGt                                                                                                             <  3:232349/149‑1 (MQ=255)                                    cTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTAt                                                                                                           >  6:139229/1‑149 (MQ=255)                                    cTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTAt                                                                                                           >  6:23564/1‑149 (MQ=255)                                    cTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTAt                                                                                                           >  2:239145/1‑149 (MQ=255)                                       ggTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGa                                                                                                        <  3:105476/149‑1 (MQ=255)                                          gtTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGa                                                                                                     >  6:82358/1‑149 (MQ=255)                                             gaccatcgTAAATCATACCAGGTCTTACGTAGAGCTACGTGGCTGCCACACGGCACTGTAAAGACGGTTCGCGTTGGTGTCTTCAGCCAGTGTCAAGCAGCTTAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATc                                                                                                   <  6:209117/142‑1 (MQ=255)                                                ccGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGc                                                                                               <  3:162307/149‑1 (MQ=255)                                                        aTCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAg                                                                                                               >  2:251047/1‑125 (MQ=255)                                                        aTCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGCGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAg                                                                                                               <  1:251047/125‑1 (MQ=255)                                                         tCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGaaa                                                                                       <  1:79908/149‑1 (MQ=255)                                                         tCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGATCCGCGTAGGTATGGACGATCTGGCTGCCGaaa                                                                                       <  1:37889/149‑1 (MQ=255)                                                            gACCAGGTCGTACGTATCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGaa                                                                                   <  5:219791/149‑1 (MQ=255)                                                            gACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGaa                                                                                   <  1:43562/149‑1 (MQ=255)                                                                 ggTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAggcg                                                                              >  2:170057/1‑149 (MQ=255)                                                                   tCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGa                                                                                                     >  8:215147/1‑124 (MQ=255)                                                                   tCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGa                                                                                                     <  7:215147/124‑1 (MQ=255)                                                                   tCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGa                                                                                                     <  7:209601/124‑1 (MQ=255)                                                                   tCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAggcggc                                                                            <  7:319513/149‑1 (MQ=255)                                                                   tCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTCTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGa                                                                                                     >  8:209601/1‑124 (MQ=255)                                                                    cgtacgtaGCGCTACCGTGATGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACACAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAggcggc                                                                            >  2:109233/1‑148 (MQ=255)                                                                      tacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGAc                                                                                                       >  5:82470/1‑119 (MQ=255)                                                                      tacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGAc                                                                                                       <  6:82470/119‑1 (MQ=255)                                                                      tacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGAc                                                                         <  5:159292/149‑1 (MQ=255)                                                                            gcgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAAc                                                                   <  5:139229/149‑1 (MQ=255)                                                                             cgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAggcg                                                                              <  6:3703/137‑1 (MQ=255)                                                                             cgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAggcg                                                                              >  5:3703/1‑137 (MQ=255)                                                                                  ccGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATg                                                               <  5:82358/147‑1 (MQ=255)                                                                                     tgctgcCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGAcgtc                                                          <  1:239145/149‑1 (MQ=255)                                                                                         gcCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGa                                                                     <  5:306032/134‑1 (MQ=255)                                                                                         gcCACACGGCACTGGCAAGACCGTTAGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGa                                                                     >  6:306032/1‑134 (MQ=255)                                                                                           cacaCGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatc                                                    <  4:128246/149‑1 (MQ=255)                                                                                                  cACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcgcatcgc                                             <  5:23564/149‑1 (MQ=255)                                                                                                   aCTGGCAAGACCGTTCGCGTTGCTGTCATCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGTCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGcc                                            <  6:113471/149‑1 (MQ=255)                                                                                                       gCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGAt                                        <  1:109233/149‑1 (MQ=255)                                                                                                         aaGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGc                                      >  8:307056/1‑149 (MQ=255)                                                                                                         aaGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGc                                      >  7:181627/1‑149 (MQ=255)                                                                                                          aGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGATCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGGTGc                                      >  5:69900/1‑148 (MQ=255)                                                                                                           gcccGGTCGCGGTGCTGGCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGACCCGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCa                                    <  8:181627/147‑1 (MQ=255)                                                                                                            aCCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCAt                                   >  5:274711/1‑149 (MQ=255)                                                                                                                ttCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTAt                                                                >  2:86165/1‑116 (MQ=255)                                                                                                                ttCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTAt                                                                <  1:86165/116‑1 (MQ=255)                                                                                                                ttCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATgcgc                               >  3:336964/1‑149 (MQ=255)                                                                                                                        gCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAg                        <  4:336964/148‑1 (MQ=255)                                                                                                                              ttCACCCAGGGTCCAGTAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCt                 >  2:177673/1‑149 (MQ=255)                                                                                                                               tCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGaa                                                                                   >  5:165806/1‑82 (MQ=255)                                                                                                                               tCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGaa                                                                                   <  6:165806/82‑1 (MQ=255)                                                                                                                                 aCCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTggg              >  6:154338/1‑149 (MQ=255)                                                                                                                                 aCCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTgg               <  1:5462/148‑1 (MQ=255)                                                                                                                                   ccAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGaaa                                                                                  <  1:131272/79‑1 (MQ=255)                                                                                                                                   ccAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGaaa                                                                                  >  2:131272/1‑79 (MQ=255)                                                                                                                                   ccAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGATAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt             >  1:160131/1‑148 (MQ=255)                                                                                                                                   ccAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGCCCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt             >  6:237203/1‑148 (MQ=255)                                                                                                                                    cAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCa           <  6:274711/149‑1 (MQ=255)                                                                                                                                        ccagagcagcTTCAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCAt                                   <  4:159669/117‑1 (MQ=255)                                                                                                                                        ccagagcagcTTCAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCAt                                   >  3:159669/5‑121 (MQ=255)                                                                                                                                        ccagagcagcTTCAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAg                        <  5:252902/128‑1 (MQ=255)                                                                                                                                        ccagagcagcTTCAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAg                        >  6:252902/5‑132 (MQ=255)                                                                                                                                        ccagagcagcTTCAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGACTGATGACATGCGCGTTGTAGGTCAGCTGGGTCAGGt        >  6:232353/5‑148 (MQ=255)                                                                                                                                           cagcagCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGATAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCt     <  2:160131/148‑1 (MQ=255)                                                                                                                                           cagcagCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCt     <  1:159450/148‑1 (MQ=255)                                                                                                                                            agcagcTGAAGCTGCTCT‑GGCTGGCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg   <  7:117571/149‑1 (MQ=255)                                                                                                                                             gcagcTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTggg  <  5:154338/149‑1 (MQ=255)                                                                                                                                             gcagcTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTggg  >  2:310094/1‑149 (MQ=255)                                                                                                                                             gcagcTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg   >  5:170364/1‑148 (MQ=255)                                                                                                                                             gcagcTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg   <  6:170364/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CAAGTTCGTCGAGTCCTACGACATCGCCGTTAA‑CCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCT‑GGCTGCCGGCGCTGAT‑CGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGG  >  NC_002947/535673‑535960

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