Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,506,923	A→G	9.5%	intergenic (+44/‑73)	petC → / → sspA	ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,506,923	0	A	G	9.5%	168.4 / 7.9	63	intergenic (+44/‑73)	petC/sspA	ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (2/4);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

AACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGG  >  NC_002947/1506776‑1507071                                                                                                                                                    |                                                                                                                                                      aaCTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAg                                                                                                                                                      >  6:150702/1‑149 (MQ=255)  aCTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAg                                                                                                                                                      >  4:285418/1‑148 (MQ=255)   cTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCa                                                                                                                                                    <  5:42936/149‑1 (MQ=255)     ggAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCaa                                                                                                                                                   >  5:142095/1‑147 (MQ=255)        aaGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt                                                                                                                                               <  3:285418/149‑1 (MQ=255)        aaGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt                                                                                                                                               <  4:296405/149‑1 (MQ=255)          gCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGc                                                                                                                                             <  7:146277/149‑1 (MQ=255)          gCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGc                                                                                                                                             <  6:152775/149‑1 (MQ=255)             aGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAg                                                                                                                                          <  7:58112/149‑1 (MQ=255)             aGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg                                                                                                                                          <  2:202582/149‑1 (MQ=255)             aGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAg                                                                                                                                          <  7:53982/149‑1 (MQ=255)              gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAgg                                                                                                                                         <  7:221692/149‑1 (MQ=255)              gcgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTCGCTCAAGCGTGAATACCGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAgg                                                                                                                                         >  2:87539/1‑149 (MQ=255)               cgcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc                                                                                                                                        <  6:142095/149‑1 (MQ=255)                gcATCGGCAACTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGTTGGCCTGGCAATACTGCTAGGcc                                                                                                                                       <  5:304563/149‑1 (MQ=255)                                    ctttttCTGGCTTACTTCTTCGTATGCGCATACTTGCTCAAGCGTCAATACTGTAAGGACGGGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGGGGCCTAGCAGTATTGCCAGGCCACCCGGGGCGCGTTTTCAt                                                                                                                    <  4:5747/143‑1 (MQ=255)                                         tCTGGCTTTCTTCTTCGTATTCGCCTACTCGCTCAAGCGTGAATACCGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTaaa                                                                                                               <  1:87539/148‑1 (MQ=255)                                             gCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTg                                                                                                           <  5:150702/148‑1 (MQ=255)                                              atttCTTCTTCGTATTCGCCTACTTGCTCAAGCTTGAATACTGGAAGGACGTGCACTTATCACGCAGTATTTTCTGGTTTTTGAACGCGCCCGGGGTGGGATAGCa                                                                                                                                                    <  6:133858/105‑1 (MQ=255)                                                   ttcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTAAATTAAACTTTCACAAg                                                                                                     >  3:190944/1‑148 (MQ=255)                                                     cttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCa                                                                                                   <  6:323070/148‑1 (MQ=255)                                                        cGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATcc                                                                                                <  6:371152/148‑1 (MQ=255)                                                           aTTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGc                                                                                            >  3:160576/1‑149 (MQ=255)                                                               gCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCCGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCga‑g                                                                                        >  1:122554/1‑148 (MQ=255)                                                                 cTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACACAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCga‑ggag                                                                                     >  1:336287/1‑149 (MQ=255)                                                                    cTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGcc                                                                                   <  3:190661/148‑1 (MQ=255)                                                                       gCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCAcccc                                                                                                                                  <  8:77887/99‑1 (MQ=255)                                                                       gCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCAcccc                                                                                                                                  >  7:77887/1‑99 (MQ=255)                                                                        cTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTAc                                                                              <  6:275648/149‑1 (MQ=255)                                                                           aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCAAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGAGA‑GGAGCCGc                                                                                 >  5:373005/1‑143 (MQ=255)                                                                           aaGAGTGAATACTGGAAGGACGTGCAATTATCACTCAGTAGTTTCTGGTTTTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGTCCACCCCGTGCTCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGc                                                                                 <  6:373005/143‑1 (MQ=255)                                                                               gTGAATACTCGAAGGACGTGCACTGATCACGCACTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCCCCCCGGCCGCCTATTCATTTAAAGTAGCACAAGCATCCGTGCGA‑GCAGCCCCTACAtcgccg                                                                       >  6:172613/1‑144 (MQ=255)                                                                                 gAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGgcgc                                                                      >  3:202704/1‑148 (MQ=255)                                                                                   aTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTGGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGcaac                                                                   <  2:275362/149‑1 (MQ=255)                                                                                            ggACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGt                                                           >  1:1990/1‑148 (MQ=255)                                                                                              aCGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAg                                                        >  5:357158/1‑149 (MQ=255)                                                                                               cGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGc                                                       >  7:181683/1‑149 (MQ=255)                                                                                                 tGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGCGCGCGTTTAGATTTAAAGTTGCACAAGCAACCGTGCGA‑GGAGCCGCTACATGGGCGCAAGCAACACGATAGCCt                                                     >  2:326924/1‑149 (MQ=39)                                                                                                   cACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTg                                                    <  6:105039/148‑1 (MQ=255)                                                                                                   cACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGc                                                   >  6:2536/1‑149 (MQ=255)                                                                                                       gATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGACGGAGCCGCTACATGGGCGCAACCAACAGGTTa                                                         <  1:44462/140‑1 (MQ=255)                                                                                                       gATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTa                                                         >  2:44462/1‑139 (MQ=255)                                                                                                       gATCACCCAGTAGTTTCTGGGTGTGGAACGCGCCCGGGGTGCCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTAtt                                               <  2:336287/149‑1 (MQ=255)                                                                                                        aTCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc                                              >  5:293297/1‑149 (MQ=255)                                                                                                               agtagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATccc                                       >  2:8884/1‑149 (MQ=255)                                                                                                               agtagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATccc                                       >  1:11273/1‑149 (MQ=255)                                                                                                                 tagtTTCTGGTTGTTGAACGCGCCCTGGGTGGCCTAGCAGTATTGCCAGGACACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGc                                     <  6:375671/149‑1 (MQ=255)                                                                                                                 tagtTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGc                                     <  1:326924/149‑1 (MQ=255)                                                                                                                     ttCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGAt                                 >  1:278796/1‑149 (MQ=255)                                                                                                                       cTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATc                                <  4:202704/148‑1 (MQ=255)                                                                                                                             gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTAtt                          <  5:2536/148‑1 (MQ=255)                                                                                                                              ttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATtc                         <  1:8884/148‑1 (MQ=255)                                                                                                                                  aCGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACGCCGGGCGCGTTGTGCTTTAAAGTTGCACAAGGATCCGTGCGA‑GGAGCCGCTACATGGGCGCACCCAACCGGTTAGCCTGCTATTGCGCTCCCGCTGATCACGATTCTCATc                    >  1:83814/1‑149 (MQ=255)                                                                                                                                       cccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGtt               <  2:11273/149‑1 (MQ=255)                                                                                                                                        ccGGGGTGGCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGCTCCCGCTGATCACTATTCTCATCGGGTTc              >  2:78837/1‑149 (MQ=38)                                                                                                                                        ccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCAt                     >  3:103450/1‑142 (MQ=255)                                                                                                                                        ccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCAt                     <  4:103450/142‑1 (MQ=255)                                                                                                                                        ccGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTc              >  1:248292/1‑149 (MQ=255)                                                                                                                                          ggggTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGc            >  8:209532/1‑149 (MQ=255)                                                                                                                                           gggTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGcc           <  2:278796/149‑1 (MQ=255)                                                                                                                                             gTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTAATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTg         >  2:143096/1‑149 (MQ=255)                                                                                                                                                  cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCt    >  7:166262/1‑149 (MQ=255)                                                                                                                                                   tAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGAGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTACTCGGGTTCGCCTGGTTCTg   >  1:291024/1‑149 (MQ=255)                                                                                                                                                    aGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCCCTACAAGGGCGCAACCAACACGTTACCCTCCTATTCCGATCCCGCTGATCACAATACTCATCGGGTTCGCCTGGGTCTgg  >  5:350035/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      AACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGA‑GGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTCTGG  >  NC_002947/1506776‑1507071

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