Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,353,340	C→G	11.0%	P141R (CCG→CGG)	yrpB →	2‑nitropropane dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,353,340	0	C	G	11.0%	187.7 / 12.7	64	P141R (CCG→CGG)	yrpB	2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (24/32);  new base G (4/3);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATG  >  NC_002947/4353195‑4353488                                                                                                                                                   |                                                                                                                                                     cAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAg                                                                                                                                                  <  1:1635415/150‑1 (MQ=255)  gggtccTTGAGCGTTATTAAGGGGAGGTGGGCGACGTTTTCCCGCCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAgg                                                                                                                                                 <  2:1633339/145‑1 (MQ=255)  aGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCGGCAgg                                                                                                                                                 >  2:237516/1‑150 (MQ=255)      aCTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGa                                                                                                                                             >  1:1100529/1‑150 (MQ=255)       cTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGaa                                                                                                                                            <  1:1340384/150‑1 (MQ=255)       cTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGaa                                                                                                                                            >  1:1564132/1‑150 (MQ=255)       cTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCAAT‑TTCGGCCTGCCGCAGGCCGaa                                                                                                                                            <  2:180151/150‑1 (MQ=255)       cTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGAAGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGaa                                                                                                                                            >  1:1691736/1‑150 (MQ=255)         tGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCGGCAGGCCGaaat                                                                                                                                          <  1:237516/150‑3 (MQ=255)         tGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCGGCAGGCCGaaat                                                                                                                                          <  1:458852/150‑3 (MQ=255)             ccTTATTACAGCGAAGTGGGCGCCGTTTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAActgctg                                                                                                                                      <  1:1672957/150‑1 (MQ=255)               ttattaCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCa                                                                                                                                    <  1:1613274/150‑1 (MQ=255)                   taCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                <  1:205861/150‑1 (MQ=255)                    aCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCggg                                                                                                                               <  2:1468282/150‑1 (MQ=255)                      aGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTa                                                                                                                             >  1:1988572/1‑150 (MQ=255)                      aGCGAAGTGGGCGCCGAGTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTa                                                                                                                             >  1:1841148/1‑150 (MQ=255)                       gCGAAGTGGGCGCAGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCAGCAGGCCGAACTGCTGCAGCGGGTaa                                                                                                                            >  1:562546/1‑150 (MQ=255)                         gAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAg                                                                                                                          <  1:670248/150‑1 (MQ=255)                         gAAGTGGGCCCTGATTTCACGCCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAg                                                                                                                          <  1:2147804/150‑1 (MQ=255)                                gcgcCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCg                                                                                                                   >  2:786781/1‑150 (MQ=255)                                      aTTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCa                                                                                                             >  2:1805436/1‑150 (MQ=255)                                           aCGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                        >  2:935772/1‑150 (MQ=255)                                            cGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTAc                                                                                                       >  2:1440227/1‑150 (MQ=255)                                             ggCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGCGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACt                                                                                                      >  1:2288090/1‑150 (MQ=255)                                                  ccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCa                                                                                                 >  1:2149097/1‑150 (MQ=255)                                                  ccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCa                                                                                                 >  1:638638/1‑150 (MQ=255)                                                    aCGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAG‑TCCAGTGTTGGAGGTGTCGCAGTGCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAgc                                                                                               >  1:1113002/1‑150 (MQ=255)                                                           tgtCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGa                                                                                        <  1:935772/150‑1 (MQ=255)                                                             tCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCCCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGGGGTTCCAAGGTACTGTCCAGAGCgacgaca                                                                                      >  1:2276985/1‑149 (MQ=255)                                                                aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                   >  1:1568298/1‑150 (MQ=255)                                                                aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                   <  2:1466880/150‑1 (MQ=255)                                                                aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGGGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCCGATAAAGGACAGGGGTGTCGAGGTACGGTCCAGCGCCACGACCGTg                                                                                   >  2:1423761/1‑150 (MQ=255)                                                                   cGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaa                                                                                <  1:2142381/150‑1 (MQ=255)                                                                   cGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaa                                                                                >  1:1554818/1‑150 (MQ=255)                                                                           tttCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGc                                                                        <  2:2149097/150‑1 (MQ=255)                                                                             tCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCTCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCct                                                                      <  2:1179515/150‑1 (MQ=255)                                                                             tCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCct                                                                      <  2:1232618/150‑1 (MQ=255)                                                                                acgaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggc                                                                   <  2:1849624/150‑1 (MQ=255)                                                                                       gagCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCg                                                                                                                   >  2:828962/1‑95 (MQ=255)                                                                                       gagCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCg                                                                                                                   <  1:828962/95‑1 (MQ=255)                                                                                            gcctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAGCTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTCCTGTCCAGCGCCACGACCGTCGAAGAAGCGGCCTGGCTGGAGCGCactc                                                       >  2:1002713/1‑147 (MQ=255)                                                                                            gcctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATg                                                       <  1:1440227/150‑1 (MQ=255)                                                                                                gctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtg                                                   >  1:1942779/1‑150 (MQ=255)                                                                                                 ctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgc                                                  <  2:1841148/150‑1 (MQ=255)                                                                                                 ctgGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgc                                                  <  2:1100529/150‑1 (MQ=255)                                                                                                      cGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgc                                             <  1:1805436/150‑1 (MQ=255)                                                                                                       gAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgcg                                            <  2:1799615/150‑1 (MQ=255)                                                                                                                   ccGTAGGTGGTGAGCGTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATg                                <  2:1113002/150‑1 (MQ=255)                                                                                                                   ccGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTc                                                                                                   >  1:380157/1‑83 (MQ=255)                                                                                                                   ccGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTc                                                                                                   <  2:380157/83‑1 (MQ=255)                                                                                                                        ggtggtGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTa                           <  2:1564132/150‑1 (MQ=255)                                                                                                                                  ttCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                 <  1:786781/150‑1 (MQ=255)                                                                                                                                  ttCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                 <  2:1988572/150‑1 (MQ=255)                                                                                                                                  ttCCAT‑TTCGGCCTGCCGCAGACAAAAGTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                 >  1:901724/1‑150 (MQ=255)                                                                                                                                   tCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                <  1:115435/150‑1 (MQ=255)                                                                                                                                   tCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGACGGTg                <  2:638638/150‑1 (MQ=255)                                                                                                                                   tCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGCCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGCAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                <  2:1682031/150‑1 (MQ=255)                                                                                                                                   tCCAT‑TGCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGGCCGTGGAAGAAGCGGCCGGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                <  2:450288/150‑1 (MQ=255)                                                                                                                                       t‑ttCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCa            >  1:341674/1‑150 (MQ=255)                                                                                                                                         ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATggg                              >  2:1119947/1‑131 (MQ=255)                                                                                                                                         ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATggg                              <  1:1119947/131‑1 (MQ=255)                                                                                                                                         ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCAt           >  1:1729844/1‑150 (MQ=255)                                                                                                                                         ttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCAt           >  1:2055397/1‑150 (MQ=255)                                                                                                                                                 gCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCAt   >  1:842517/1‑150 (MQ=255)                                                                                                                                                  ccGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATg  >  1:250916/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAT‑TTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATG  >  NC_002947/4353195‑4353488

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