Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,045,262	G→A	7.3%	intergenic (+35/+71)	PP_0905 → / ← PP_0906	hypothetical protein/RND family multidrug transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,045,262	0	G	A	7.3%	205.2 / 4.0	69	intergenic (+35/+71)	PP_0905/PP_0906	hypothetical protein/RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base A (2/3);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTA  >  NC_002947/1045122‑1045408                                                                                                                                              |                                                                                                                                                    gTTGCCGCGGGATGCAAAACCGCAGCAGCAGTCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAAGCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCggg                                                                                                                                             <  4:261226/149‑1 (MQ=255) gTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCggg                                                                                                                                             <  5:44714/149‑1 (MQ=255)      cgcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCtt                                                                                                                                        >  2:30477/1‑149 (MQ=255)      cgcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCtt                                                                                                                                        <  7:139678/149‑1 (MQ=255)       gcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCttt                                                                                                                                       <  7:195030/149‑1 (MQ=255)        cgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTg                                                                                                                                      <  6:226453/149‑1 (MQ=255)           gATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtt                                                                                                                                   <  7:218198/149‑1 (MQ=255)               cAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAACGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGttttt                                                                                                                                <  1:60225/149‑1 (MQ=255)               cAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGGCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                               <  1:209770/149‑1 (MQ=255)               cAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                               <  4:118957/149‑1 (MQ=255)               cAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                               <  1:154584/149‑1 (MQ=255)                 aaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGc                                                                                                                             >  3:290920/1‑149 (MQ=255)                   aCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGt                                                                                                                           >  2:195257/1‑149 (MQ=255)                      gcagcagcaGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTg                                                                                                                        >  3:307186/1‑149 (MQ=255)                       cagcagcagGCCCGGAGCGAGCCAGACCCGGACCTTATGCTGGTCCTGTTTGGCCTGCGGCAAGGGGTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGc                                                                                                                       <  4:239198/149‑1 (MQ=255)                        agcagcagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGcc                                                                                                                      <  1:195257/149‑1 (MQ=255)                             cagGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACACTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCGTTTGGTTTTTTGCGTCTGCCTTGcc                                                                                                                 >  6:59454/1‑149 (MQ=255)                                gCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGc                                                                                                              <  4:192639/149‑1 (MQ=255)                                  gcggttcgagcCTGACCTGGTCCTTTTGCTGTTCCTGTTTGTCGTGCGTCAAGGGTTGCGGTGACAGGGTTGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGGCTCGGGATTGGGGGTCCTTTGGTTTTTTGCGTCTGGGTTGGCGGCct                                                                                                            <  4:79980/143‑1 (MQ=255)                                         gagcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGAGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgcg                                                                                                     >  2:12742/1‑149 (MQ=255)                                           gcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCggg                                                                                                   >  1:16429/1‑149 (MQ=255)                                             cAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCAATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCa                                                                                                 <  2:108608/149‑1 (MQ=255)                                                 ccAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCgggtaaacc                                                                                             <  7:68897/149‑7 (MQ=255)                                                 ccAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgccc                                                                                             <  4:307186/149‑1 (MQ=255)                                                    ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGcc                                                                                                                      >  3:69945/1‑121 (MQ=255)                                                    ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGcc                                                                                                                      <  4:69945/121‑1 (MQ=255)                                                    ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgcc                                                                                          >  3:203501/1‑149 (MQ=255)                                                             gCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGc                                                                                                              <  6:12222/120‑1 (MQ=255)                                                             gCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGc                                                                                                              >  5:12222/1‑120 (MQ=255)                                                                gtaccTCTTTGGCCAGCGTCAAGGGTTGCCGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGACGCCACGCCATAGCGGGGCCTTTGGTTTTTTGCGTCTTCCTTGCAGGCCTCTTCCCGGGCAGACCGCCCCAAGGACGGtt                                                                              >  8:250589/4‑149 (MQ=255)                                                                 ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTc                                                                             <  2:196585/149‑1 (MQ=255)                                                                 ttCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTc                                                                             <  3:74833/149‑1 (MQ=255)                                                                    cTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAg                                                                          <  2:49112/149‑1 (MQ=255)                                                                    cTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGAGGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCTCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAg                                                                          <  1:111832/149‑1 (MQ=255)                                                                      gTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGc                                                                         >  5:338819/1‑148 (MQ=255)                                                                           gCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgc                                                                                           >  4:102661/1‑125 (MQ=255)                                                                           gCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgc                                                                                           <  3:102661/125‑1 (MQ=255)                                                                             cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑g                                                                >  4:1310/1‑149 (MQ=255)                                                                             cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑g                                                                >  6:308530/1‑149 (MQ=255)                                                                             cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑g                                                                <  1:30477/149‑1 (MQ=255)                                                                             cTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGAGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGTCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑g                                                                <  1:12742/149‑1 (MQ=255)                                                                                cGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGt                                                              <  6:5508/148‑1 (MQ=255)                                                                                 gTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGc                                                            >  3:277986/1‑149 (MQ=255)                                                                                  tCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGc                                                            >  2:80972/1‑148 (MQ=255)                                                                                    aaGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAg                                                                          <  4:248718/133‑1 (MQ=255)                                                                                    aaGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAg                                                                          >  3:248718/1‑133 (MQ=255)                                                                                       ggTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGgc                                                      <  2:16429/149‑1 (MQ=255)                                                                                           gCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTCCA‑GGTGCCTTGCCGCta                                                  >  8:322586/1‑148 (MQ=255)                                                                                            cGGTCACACGCATGACCCCGGCGTTTGCTGGAAAA‑CCCAGAGGCCCCGCGATTGCGGGGCCTTGGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCg                                                 <  4:277986/149‑1 (MQ=255)                                                                                              gTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGcc                                               >  7:53575/1‑149 (MQ=255)                                                                                              gTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGcc                                               >  4:351341/1‑149 (MQ=255)                                                                                                 acaGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTtg                                            >  8:27267/1‑149 (MQ=255)                                                                                                      gATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTgcgc                                       <  5:59454/149‑1 (MQ=255)                                                                                                       aTGAACCCTGCGGTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGcc                                      <  7:250589/149‑1 (MQ=255)                                                                                                         gAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCcag                                    >  6:255761/1‑149 (MQ=255)                                                                                                               tGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                              >  2:232242/1‑149 (MQ=255)                                                                                                                 cGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgtgt                            <  1:232242/149‑1 (MQ=255)                                                                                                                  gTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTa                           <  4:290920/149‑1 (MQ=255)                                                                                                                          gAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCAc                    <  1:19604/148‑1 (MQ=255)                                                                                                                          gAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACg                   >  3:332724/1‑149 (MQ=255)                                                                                                                            aaa‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGc                 <  3:1310/149‑1 (MQ=255)                                                                                                                            aaa‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGGCTGCCTTGCCGGCCTCTTCGCGGGCATCCCGCCCCAAGGTCGGTTCAAGCGATTGCATTGTGCCATGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACgg                  <  1:80972/149‑1 (MQ=255)                                                                                                                                gCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCaa               >  5:304262/1‑148 (MQ=255)                                                                                                                                   aaagGCCCCGCAATCGCGTGGACTTTGGTTTGTTGTGTGTGCCTTGGCGGCCTCTTCGGGTGAAGCGCGACCCAAGGTCTTTTTAAGAGATTGCT‑GGTTCCTTGGCGGTCGCCTTGTGCGCCAGCAGTGTGTa                           <  6:359918/131‑1 (MQ=255)                                                                                                                                   aaagGCCCCGCAATCGCGACGCCTTTGGTTTTTTGTGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGAATGCA‑GGTGCCTTGGCGCTCGGCTTGTGCGCCAGCAGGGTGTa                           >  5:359918/3‑133 (MQ=255)                                                                                                                                     agGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                <  5:251772/126‑1 (MQ=255)                                                                                                                                     agGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAgg                                >  6:251772/1‑126 (MQ=255)                                                                                                                                       gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaac        >  5:202140/1‑148 (MQ=255)                                                                                                                                         cTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGCAGATGCACGGCAACACAAACAAg     >  5:317691/1‑149 (MQ=255)                                                                                                                                            gcgATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTa  <  5:255761/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    GTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCGGTCACAGGGATGAACCCTGCGTTAGCTGGAAAA‑GCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCA‑GGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGGTA  >  NC_002947/1045122‑1045408

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