Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,045,258	T→C	8.3%	intergenic (+31/+75)	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,258	0	T	C	8.3%	165.0 / 4.9	60	intergenic (+31/+75)	PP_0905/PP_0906	hypothetical protein/RND family multidrug transporter
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base C (3/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

CCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGG  >  NC_002947/1045113‑1045406                                                                                                                                                     |                                                                                                                                                     ccAGTGCGTGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTcgc                                                                                                                                                   <  5:250701/149‑1 (MQ=255)         cgTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCgg                                                                                                                                           <  3:35425/149‑1 (MQ=255)                gcgGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCttt                                                                                                                                    >  4:362817/1‑149 (MQ=255)                   ggATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCCCGCAATCGCGAGGCCTCTggc                                                                                                                                 <  7:5102/149‑2 (MQ=255)                       gCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGttttt                                                                                                                             >  2:255280/1‑149 (MQ=255)                        cAAAACCGCAGCAGCATGCCCGGAGCGAG‑CAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTg                                                                                                                           <  3:57743/149‑1 (MQ=255)                        cAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                            <  3:244594/149‑1 (MQ=255)                        cAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTGTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGtttttt                                                                                                                            <  4:229776/149‑1 (MQ=255)                         aaaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTGCCGGACCGGGACGGGAGAACCCTGCTTTAGATGAAAAAACAACAGGACTCCCGATTGAGGGGCCTTTGGtttt                                                                                                                              >  1:122360/1‑149 (MQ=255)                         aaaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTg                                                                                                                           >  8:50296/1‑149 (MQ=255)                         aaaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTg                                                                                                                           >  7:14656/1‑149 (MQ=255)                          aaaCCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGc                                                                                                                          <  3:73996/149‑1 (MQ=255)                             ccGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTc                                                                                                                       >  2:52844/1‑149 (MQ=255)                              cGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTGCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCt                                                                                                                      >  3:33203/1‑149 (MQ=255)                                         gCCCGGAGCGAGCCAGAGCAGGACGTATTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGc                                                                                                           <  8:310791/149‑1 (MQ=255)                                         gCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGc                                                                                                           <  6:104325/149‑1 (MQ=255)                                                  gagcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTcgcg                                                                                                  <  8:212117/149‑1 (MQ=255)                                                    gcCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCggg                                                                                                >  5:300307/1‑149 (MQ=255)                                                     ccAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGAGATTGAGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCggg                                                                                                >  1:217812/1‑148 (MQ=255)                                                       aGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAg                                                                                             >  8:214226/1‑149 (MQ=255)                                                       aGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAg                                                                                             >  7:272800/1‑149 (MQ=255)                                                           cAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGc                                                                                                               <  6:253881/127‑1 (MQ=255)                                                           cAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccg                                                                                         >  6:288141/1‑149 (MQ=255)                                                             ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgcc                                                                                       <  7:265624/149‑1 (MQ=255)                                                             ggACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAgcccgcc                                                                                       <  3:9592/149‑1 (MQ=255)                                                                      gCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCg                                                                              >  4:259803/1‑149 (MQ=255)                                                                             cTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCaa                                                                        <  1:52844/148‑1 (MQ=255)                                                                                tttGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGa                                                                    <  1:255280/149‑1 (MQ=255)                                                                                         cGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGt                                                            <  7:214226/148‑1 (MQ=255)                                                                                           tCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGcc                                                         >  3:222598/1‑149 (MQ=255)                                                                                             aaGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGt                                                                            >  8:289820/1‑128 (MQ=255)                                                                                             aaGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTAAAGCGATTGCAGGTGCCtt                                                       >  2:260925/1‑149 (MQ=255)                                                                                                ggTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGgc                                                    <  6:300307/149‑1 (MQ=255)                                                                                                  ttG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGgcgc                                                  >  3:89053/1‑149 (MQ=255)                                                                                                      ctgTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTc                                                <  6:123242/146‑1 (MQ=255)                                                                                                             acaGGGATGAACACTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGGCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtg                                          >  6:182986/1‑149 (MQ=255)                                                                                                              caggtctGAAGCGTGCGTTAGGTGGAAAAGCCTGAGGCGTTGCGATTGGTTGTCCTTTGGTTTTTTGCGTCTGCCTTGCAGGCCTTTTCGCTGGCAGCGCGCCACATGGTATTTTAAAGCGATTTCAGGTGCATTTGCGCTCGCCTtgt                                         <  2:316326/143‑1 (MQ=255)                                                                                                              caGGGATGAACCCTGCGTTAGCGGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTtgt                                         >  7:237198/1‑149 (MQ=255)                                                                                                               aGGGATGAACACTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGACTTTGTTTTTTTGCGTCTGCATTGCCGGCCTCTTCGCGGGCAGCCCGACCCAAGGTCGGTTCAAGCGATTGCAGGTGCATTGGCGCTCTCCTtgtg                                        >  8:146284/1‑149 (MQ=255)                                                                                                                 ggATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGAGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTgcg                                      >  3:307625/1‑149 (MQ=255)                                                                                                                 ggATGAAACCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCAGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTgcg                                      >  3:332156/1‑149 (MQ=255)                                                                                                                     gAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCGTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGTTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcag                                  >  7:65134/1‑149 (MQ=255)                                                                                                                        cccTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCcagcag                               >  6:173516/1‑149 (MQ=255)                                                                                                                          cTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCg                                               >  2:348021/1‑131 (MQ=255)                                                                                                                          cTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCg                                               <  1:348021/131‑1 (MQ=255)                                                                                                                             cGTTAGCAGGAAAAGCCAGAGGCCTCGCGATGGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                            <  4:33203/147‑1 (MQ=255)                                                                                                                                  gCTGGAAAAGCCAGCGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGAt                      <  5:302956/148‑1 (MQ=255)                                                                                                                                  gCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGAt                      <  4:222598/148‑1 (MQ=255)                                                                                                                                   cTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATg                     >  1:528/1‑148 (MQ=255)                                                                                                                                      gAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACg                 >  5:319715/1‑149 (MQ=255)                                                                                                                                          aaccAAAGGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCaa             >  8:3843/3‑149 (MQ=255)                                                                                                                                           gCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAac            <  4:89053/149‑1 (MQ=255)                                                                                                                                           gCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAac            <  3:293076/149‑1 (MQ=255)                                                                                                                                            ccAAAGGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                            <  7:167053/132‑1 (MQ=255)                                                                                                                                            ccAAAGGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGgt                            >  8:167053/1‑132 (MQ=255)                                                                                                                                              aaagGCCCCGCAATCGCGAGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAacac          >  3:233414/3‑148 (MQ=255)                                                                                                                                                agGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaa       <  5:135598/149‑1 (MQ=255)                                                                                                                                                agGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaa       <  5:288141/149‑1 (MQ=255)                                                                                                                                                  gCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACacaaac      >  7:208083/1‑148 (MQ=255)                                                                                                                                                    cTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAg   <  2:122360/149‑1 (MQ=255)                                                                                                                                                     tCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAgg  >  2:52261/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCAGTGCGCGTTGCCGCGGGATGCAAAACCGCAGCAGCAGGCCCGGAGCGAGCCAGACCAGGACCTTTTGCTGTTCCTGTTTGGCCTGCGTCAAGGGTTG‑CGGTC‑‑ACAGGGATGAACCCTGCGTTAGCTGGAAAAGCCAGAGGCCTCGCGATTGCGGGGCCTTTGGTTTTTTGCGTCTGCCTTGCCGGCCTCTTCGCGGGCAGCCCGCCCCAAGGTCGGTTCAAGCGATTGCAGGTGCCTTGGCGCTCGCCTTGTGCGCCAGCAGGGTGTAGATGCACGGCAACACAAACAAGG  >  NC_002947/1045113‑1045406

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