Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	542,698	G→C	R312P (CGT→CCT)	rpoC →	DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	542,698	0	G	C	98.5%	228.3 / NA	67	R312P (CGT→CCT)	rpoC	DNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (33/33);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCG  >  NC_002947/542560‑542832                                                                                                                                           |                                                                                                                                       aCGACCTGTATCGCCGGGTGATCCACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCa                                                                                                                             <  1:1475043/150‑1 (MQ=255)   gACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATc                                                                                                                           <  2:78575/150‑1 (MQ=255)      cTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTc                                                                                                                                   <  2:1532642/139‑1 (MQ=255)      cTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTc                                                                                                                                   >  1:1532642/1‑139 (MQ=255)       tGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCGGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTg                                                                                                                       <  1:710832/150‑1 (MQ=255)              ccGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGaa                                                                                                                <  2:1612803/150‑1 (MQ=255)                 ggTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACaa                                                                                                             <  2:1326775/150‑1 (MQ=255)                     aTCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGCACAACGGCCGTCCTGGTCGCGCGATCACTGGCTCGAACAAGCGt                                                                                                         <  2:677055/150‑1 (MQ=255)                      tCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTc                                                                                                        >  1:401517/1‑150 (MQ=255)                          ccGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCt                                                                                                    >  1:411998/1‑150 (MQ=255)                           cGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTg                                                                                                   <  2:1478027/150‑1 (MQ=255)                             tAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGaa                                                                                                 <  2:401517/150‑1 (MQ=255)                                caacCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGACGTGCTGGTCGCGCCATCACTGGGTAGAACAAGCGTCAGCTGACGTc                                                                                              >  1:617911/1‑150 (MQ=255)                                caacCGTCTGAAGCGCCTGCTGGATCTGTCGGCACCGGACATCATCGGGCGCAACGAAAAGCGCATACATTTTGTAGCGGTTCAATCAGTGCGATGCAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAACATCGGCTGAAGTc                                                                                              >  1:1338633/1‑150 (MQ=255)                                 aacCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGaa                                                                                                 <  2:269574/146‑1 (MQ=255)                                 aacCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGaa                                                                                                 >  1:269574/1‑146 (MQ=255)                                 aacCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTcc                                                                                             <  2:411998/150‑1 (MQ=255)                                 aacCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTcc                                                                                             >  1:1328790/1‑150 (MQ=255)                                 aacCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTcc                                                                                             >  2:582824/1‑150 (MQ=255)                                   ccGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCt                                                                                           <  2:1328790/150‑1 (MQ=255)                                   cagTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGATGCAGGAGGCGGTTCCAGCTCTGCTGGCCAGCTGCCGTCGTGGCCGCGCTATCATTTGCTCGAACAAGCGTCCGCTGAAGTCCCt                                                                                           <  1:1473920/148‑1 (MQ=255)                                    cGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTg                                                                                          >  1:1479908/1‑150 (MQ=255)                                           aGCGCCTGCTGGATCTGTCGGCGCCGGACATCCTTGTGATCTGCTGAAAGCGCGTTTTGCAGGAAGCGGTTGCCGCTCTGCTGGACAACGTCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACa                                                                                   <  1:1585420/150‑1 (MQ=255)                                               ccTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGAt                                                                               >  1:110484/1‑150 (MQ=255)                                                ctgctgGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACaa                                                                                                             <  2:1077922/119‑1 (MQ=255)                                                ctgctgGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACaa                                                                                                             >  1:1077922/1‑119 (MQ=255)                                                    tgGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAg                                                                          <  1:973536/150‑1 (MQ=255)                                                    tgGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAg                                                                          <  1:582824/150‑1 (MQ=255)                                                     ggATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAgg                                                                         >  1:1223116/1‑150 (MQ=255)                                                       aTCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTa                                                                       >  2:24660/1‑150 (MQ=255)                                                        tCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTaa                                                                      <  2:1479908/150‑1 (MQ=255)                                                        tCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTaa                                                                      <  1:24660/150‑1 (MQ=255)                                                            tCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCAt                                                                                                                            >  2:1475095/1‑92 (MQ=255)                                                            tCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCAt                                                                                                                            <  1:1475095/92‑1 (MQ=255)                                                            tCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCaa                                                                  >  1:904787/1‑150 (MQ=255)                                                             cGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAg                                                                 <  1:221080/150‑1 (MQ=255)                                                              ggCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATc                                                                              <  1:1459879/136‑1 (MQ=255)                                                              ggCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATc                                                                              >  2:1459879/1‑136 (MQ=255)                                                                cgcCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTc                                                              <  2:1403539/150‑1 (MQ=255)                                                                 gcCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCg                                                             >  1:829238/1‑150 (MQ=255)                                                                 gcCGGACATCATCGTGCGCAACGAAAAACGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTGCGCTGAAGTCCCTCGCCGACATGATGAAAGGTAAGCAAGGTCg                                                             >  2:621055/1‑150 (MQ=255)                                                                          catcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCa                                                    >  1:229186/1‑150 (MQ=255)                                                                          catcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCa                                                    <  2:260691/150‑1 (MQ=255)                                                                           atcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAg                                                   >  2:1459061/1‑150 (MQ=255)                                                                            tcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGc                                                                                                           <  2:245232/93‑1 (MQ=255)                                                                            tcGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGc                                                                                                           >  1:245232/1‑93 (MQ=255)                                                                             cGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGaa                                                 <  1:808597/150‑1 (MQ=255)                                                                                gcgcAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACtt                                              <  2:1338633/150‑1 (MQ=255)                                                                                  gcAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGc                                            >  1:659992/1‑150 (MQ=255)                                                                                    aaCGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTc                                          >  2:839449/1‑150 (MQ=255)                                                                                    aaCGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTc                                          >  1:681598/1‑150 (MQ=255)                                                                                    aaCGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTc                                          >  2:1118691/1‑150 (MQ=255)                                                                                          aaGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAg                                    >  2:1148010/1‑150 (MQ=255)                                                                                            gcgcATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACACGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCg                                  >  1:588498/1‑150 (MQ=255)                                                                                                aTGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGtt                              <  2:1223116/150‑1 (MQ=255)                                                                                                    tgcAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGtt                              >  1:340689/1‑146 (MQ=255)                                                                                                    tgcAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGtt                              <  2:340689/146‑1 (MQ=255)                                                                                                           aGCGGTTGACGCTCTGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCgg                   <  2:659992/150‑1 (MQ=255)                                                                                                            gCGGTTGACGCTTTGCTGGTCAACGGCCGCCCTGGCCACGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGc                  >  2:1249197/1‑150 (MQ=255)                                                                                                                      ctctGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGc                                                                                        >  1:1391774/1‑70 (MQ=255)                                                                                                                      ctctGCTGGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGc                                                                                        <  2:1391774/70‑1 (MQ=255)                                                                                                                      ctctGCTGGAAAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTaa        <  1:1148010/150‑1 (MQ=255)                                                                                                                         tgctgGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGtt                              <  2:500682/125‑1 (MQ=255)                                                                                                                         tgctgGACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGtt                              >  1:500682/1‑125 (MQ=255)                                                                                                                            tgGACAACGGCCTTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTTAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCg  >  2:877170/1‑150 (MQ=255)                                                                                                                              gACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGc                  >  1:1153364/1‑132 (MQ=255)                                                                                                                              gACAACGGCCGTCCTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGc                  <  2:1153364/132‑1 (MQ=255)                                                                                                                                        gtcctGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGaa                                                                                                 >  1:936394/1‑43 (MQ=255)                                                                                                                                           |                                                                                                                                       ACGACCTGTATCGCCGGGTGATCAACCGTAACAACCGTCTGAAGCGCCTGCTGGATCTGTCGGCGCCGGACATCATCGTGCGCAACGAAAAGCGCATGCTGCAGGAAGCGGTTGACGCTCTGCTGGACAACGGCCGTCGTGGTCGCGCCATCACTGGCTCGAACAAGCGTCCGCTGAAGTCCCTGGCCGACATGATCAAAGGTAAGCAAGGTCGTTTCCGTCAGAACTTGCTCGGTAAGCGTGTTGACTACTCCGGCCGTTCGGTAATTACCG  >  NC_002947/542560‑542832

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