Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,911,840	C→G	6.5%	R129P (CGT→CCT)	ccmH ←	cytochrome c heme lyase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,911,840	0	C	G	6.5%	194.4 / 5.1	62	R129P (CGT→CCT)	ccmH	cytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base G (2/2);  total (33/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.37e-01

CGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATCTGC  >  NC_002947/4911691‑4911989                                                                                                                                                       |                                                                                                                                                      cGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCAc                                                                                                                                                       <  2:1624586/150‑1 (MQ=255) cGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCAc                                                                                                                                                       <  2:1281929/150‑1 (MQ=255)  gggAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACg                                                                                                                                                      <  1:305820/150‑1 (MQ=255)   ggAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGCCGGCACGCTCGCGTTCTTAGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACgg                                                                                                                                                     >  2:387452/1‑150 (MQ=255)             aaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCCTCCGCAGACAACTCATCGGCGCTCCGT‑TCGGCCGTAGCACGCTGCCTGCGtg                                                                                                                                           <  2:602443/150‑2 (MQ=255)                    aGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCAc                                                                                                                                    <  1:1256502/150‑1 (MQ=255)                     gggCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCAcc                                                                                                                                   >  1:221323/1‑150 (MQ=255)                       gCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGc                                                                                                                                 >  2:224327/1‑150 (MQ=255)                       gCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTCTATCCCGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGc                                                                                                                                 >  2:222563/1‑150 (MQ=255)                         aGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCa                                                                                                                               >  1:113424/1‑150 (MQ=255)                            gAGCAGCAGAAGGCCCGCGCTAAGCCAGAAGTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGca                                                                                                                            >  1:965084/1‑150 (MQ=255)                                 gcagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACg                                                                                                                       >  1:996709/1‑150 (MQ=255)                                  cagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGGCGcc                                                                                                                                                 >  1:716211/1‑123 (MQ=38)                                  cagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGGCGcc                                                                                                                                                 <  2:716211/123‑1 (MQ=38)                                   agAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                                     <  1:696616/150‑1 (MQ=255)                                            gcgcTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCACCATCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                                            >  1:1577486/1‑150 (MQ=255)                                                    ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCAc                                                                                                                <  2:293421/138‑1 (MQ=255)                                                    ccAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCAc                                                                                                                >  1:293421/1‑138 (MQ=255)                                                       gAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATcc                                                                                                 >  2:1331501/1‑150 (MQ=255)                                                         aTTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATcccc                                                                                               <  2:1624594/150‑1 (MQ=255)                                                             aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAACCGCCCGGCCATCAAGATCCCCGGGc                                                                                           <  2:965084/150‑1 (MQ=255)                                                              aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATc                                                                                                  <  2:889859/142‑1 (MQ=255)                                                              aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATc                                                                                                  >  1:889859/1‑142 (MQ=255)                                                              aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGcc                                                                                          >  2:240419/1‑150 (MQ=255)                                                               tCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCa                                                                                         >  1:796516/1‑150 (MQ=255)                                                                     gTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCa                                                                                                                               <  2:1525730/106‑1 (MQ=255)                                                                     gTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCa                                                                                                                               >  1:1525730/1‑106 (MQ=255)                                                                         gTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                               >  2:534091/1‑150 (MQ=255)                                                                         gTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                               >  1:1120893/1‑150 (MQ=255)                                                                          ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagcag                                                                              >  1:91318/1‑150 (MQ=255)                                                                            cTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGcc                                                                            >  1:759028/1‑150 (MQ=255)                                                                             tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAAGTCAGCGGCGCCCTGC‑ACGGCCGTAGCAAGGCGCCGGCGTACGATCACCGCCAGCACGACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCa                                                                           >  1:785466/1‑150 (MQ=255)                                                                             tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAg                                                                                                                              <  1:938115/99‑1 (MQ=255)                                                                             tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAg                                                                                                                              >  2:938115/1‑99 (MQ=255)                                                                             tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                               >  2:1039759/1‑146 (MQ=255)                                                                             tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                               <  1:1039759/146‑1 (MQ=255)                                                                                ttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCCTAGCACGGCGCCGGCGTCCGATCACCTCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGt                                                                        <  2:108013/150‑1 (MQ=255)                                                                                   cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAg                                                                                                         <  1:95381/114‑1 (MQ=255)                                                                                   cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAg                                                                                                         >  2:95381/1‑114 (MQ=255)                                                                                   cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCaa                                                                                                      >  1:1092514/1‑117 (MQ=255)                                                                                   cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCaa                                                                                                      <  2:1092514/117‑1 (MQ=255)                                                                                               gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                                            <  2:917805/99‑1 (MQ=255)                                                                                               gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                                            >  1:917805/1‑99 (MQ=255)                                                                                               gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                                               >  2:604949/1‑144 (MQ=255)                                                                                               gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                                               <  1:604949/144‑1 (MQ=255)                                                                                                      tCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTg                                                  <  2:104397/150‑1 (MQ=255)                                                                                                        gcgTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                                                <  2:796516/150‑1 (MQ=255)                                                                                                          gTTCTTCGGCAGACAACTCATCGGCGACCTGC‑CCGGCCCTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAgc                                              <  2:1365095/150‑1 (MQ=255)                                                                                                                    aGACAACTCATCGGCGCCCTGC‑ACGGGCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTc                                    <  2:785466/150‑1 (MQ=255)                                                                                                                    aGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAACCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTGGCGCACGAAGTc                                    <  2:221323/150‑1 (MQ=255)                                                                                                                     gACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                               <  2:959020/106‑1 (MQ=255)                                                                                                                     gACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                               >  1:959020/1‑106 (MQ=255)                                                                                                                        aaCTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCa                                >  2:1463699/1‑150 (MQ=255)                                                                                                                          cTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATa                              >  2:439185/1‑150 (MQ=255)                                                                                                                           tCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg                             <  2:996709/150‑1 (MQ=255)                                                                                                                              tCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCgg                          <  1:439185/150‑1 (MQ=255)                                                                                                                               cGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                                             <  1:606370/82‑1 (MQ=25)                                                                                                                               cGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                                             >  2:606370/1‑82 (MQ=25)                                                                                                                               cGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt                         <  1:534091/150‑1 (MQ=255)                                                                                                                                   gcCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAc                     <  2:1505332/150‑1 (MQ=255)                                                                                                                                     ccTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCa                   >  1:953857/1‑150 (MQ=255)                                                                                                                                     ccTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCACCAAGGCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTAGACCa                   >  1:928059/1‑150 (MQ=255)                                                                                                                                         c‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGaa               <  1:1463699/150‑1 (MQ=255)                                                                                                                                                     cACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATct    >  1:1108609/1‑150 (MQ=255)                                                                                                                                                       cGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGACGTCGCCATAGCGGTCGACCATGAAGTCGATGATctgc  <  1:1331501/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATCTGC  >  NC_002947/4911691‑4911989

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