Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,911,836	Δ1 bp	6.1%	coding (390/477 nt)	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,836	0	A	.	6.1%	201.7 / 7.3	66	coding (390/477 nt)	ccmH	cytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base . (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

TCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAA  >  NC_002947/4911690‑4911976                                                                                                                                                    |                                                                                                                                             tCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCa                                                                                                                                           >  1:1505332/1‑150 (MQ=255)  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                                                                                                                                    <  2:716211/123‑1 (MQ=38)                                   cagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGGCGcc                                                                                                                                    >  1:716211/1‑123 (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                                                                                     >  1:889859/1‑142 (MQ=255)                                                               aTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATc                                                                                     <  2:889859/142‑1 (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)                                                                         tGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAg                                                                                                                                             <  1:564218/76‑1 (MQ=255)                                                                         tGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAg                                                                                                                                             >  2:564218/1‑76 (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                                                                                                                 >  2:938115/1‑99 (MQ=255)                                                                              tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAg                                                                                                                 <  1:938115/99‑1 (MQ=255)                                                                              tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                  <  1:1039759/146‑1 (MQ=255)                                                                              tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                                  >  2:1039759/1‑146 (MQ=255)                                                                                 ttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCCTAGCACGGCGCCGGCGTCCGATCACCTCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGt                                                           <  2:108013/150‑1 (MQ=255)                                                                                    cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAg                                                                                            >  2:95381/1‑114 (MQ=255)                                                                                    cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAg                                                                                            <  1:95381/114‑1 (MQ=255)                                                                                    cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCaa                                                                                         >  1:1092514/1‑117 (MQ=255)                                                                                    cagcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCaa                                                                                         <  2:1092514/117‑1 (MQ=255)                                                                                                gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                               >  1:917805/1‑99 (MQ=255)                                                                                                gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGc                                                                                               <  2:917805/99‑1 (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                                                                                >  2:606370/1‑82 (MQ=25)                                                                                                                                cGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                                <  1:606370/82‑1 (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)                                                                                                                                                    |                                                                                                                                             TCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAA  >  NC_002947/4911690‑4911976

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