Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,024,924	T→G	6.0%	P585P (CCA→CCC)	PP_2641 ←	iron‑sulfur cluster‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,024,924	0	T	G	6.0%	191.4 / 3.0	66	P585P (CCA→CCC)	PP_2641	iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base G (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 = 9.09e-01

GGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTCCTGGCCGTACAACA  >  NC_002947/3024778‑3025072                                                                                                                                                   |                                                                                                                                                       ttccttcATGCACTTGGCGCATGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGc                                                                                                                                                     <  2:746941/144‑1 (MQ=255)  gccttcATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCa                                                                                                                                                    >  1:746941/6‑150 (MQ=255)    cctatctcacaAGGCGCAGTCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATCGTAACCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGGTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGc                                                                                                                                                  <  1:725870/143‑1 (MQ=255)       cATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGcccc                                                                                                                                               >  2:955065/1‑150 (MQ=255)         tGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCg                                                                                                                                             <  2:517547/150‑1 (MQ=255)          gCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCAGGGCAGcttgata                                                                                                                                            <  1:326862/150‑7 (MQ=255)            acaAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAATGGCCGATGCGGTATCAATCTGCCCTGGCAGccagatcgg                                                                                                                                          >  2:725870/1‑143 (MQ=255)            acaAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTgg                                                                                                                                                      >  1:759815/1‑138 (MQ=255)            acaAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTgg                                                                                                                                                      <  2:759815/138‑1 (MQ=255)               aGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGc                                                                                                                                       <  1:1171779/150‑1 (MQ=255)                   gcAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTc                                                                                                                                   <  2:1066729/150‑1 (MQ=255)                                aGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGt                                                                                                                      >  2:1187814/1‑150 (MQ=255)                                  gtcgtcAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCt                                                                                                                    >  1:968594/1‑150 (MQ=255)                                   tcgtcAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGATGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCtt                                                                                                                   <  2:1032492/150‑1 (MQ=255)                                      tcAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAg                                                                                                                           >  2:709896/1‑139 (MQ=255)                                      tcAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAg                                                                                                                           <  1:709896/139‑1 (MQ=255)                                           ggTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAAt                                                                                                           >  2:1125745/1‑150 (MQ=255)                                            gTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATg                                                                                                          >  2:1285222/1‑150 (MQ=255)                                               aGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGccc                                                                                                       <  1:926668/150‑1 (MQ=255)                                               aGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGccc                                                                                                       >  2:1220829/1‑150 (MQ=255)                                               aGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTCGACACGACGCTCTTAAGCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGccc                                                                                                       <  2:748812/150‑1 (MQ=255)                                                   tgcgtACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGAc                                                                                                   >  1:895457/1‑150 (MQ=255)                                                   tgcgtACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGAc                                                                                                   <  1:1187814/150‑1 (MQ=255)                                                       tACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGt                                                                                               >  1:165884/1‑150 (MQ=255)                                                              aTGGTAATCGTTCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCtt                                                                                                                                         <  1:27731/101‑1 (MQ=255)                                                              aTGGTAATCGTTCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCtt                                                                                                                                         >  2:27731/1‑101 (MQ=255)                                                              aTGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTg                                                                                                                                                       >  2:1312369/1‑87 (MQ=255)                                                              aTGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTg                                                                                                                                                       <  1:1312369/87‑1 (MQ=255)                                                               tGGTAATCGTCCTCGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAgg                                                                                       <  1:1003366/150‑1 (MQ=255)                                                                ggTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAggg                                                                                      <  2:165884/150‑1 (MQ=255)                                                                 gTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGt                                                                                     <  1:955065/150‑1 (MQ=255)                                                                   aaTCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAggg                                                                                      <  2:336668/147‑1 (MQ=255)                                                                   aaTCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAggg                                                                                      >  1:336668/1‑147 (MQ=255)                                                                     tCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGca                                                                                 <  1:1125745/150‑1 (MQ=255)                                                                         ccTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGc                                                                                                                                 <  2:1215521/98‑1 (MQ=255)                                                                         ccTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGc                                                                                                                                 >  1:1215521/1‑98 (MQ=255)                                                                          cTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCAGGGCAGctt                                                                                                                                                >  2:1199319/1‑80 (MQ=255)                                                                          cTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCAGGGCAGctt                                                                                                                                                <  1:1199319/82‑3 (MQ=255)                                                                                 ggTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTAtctc                                                                   <  2:968594/150‑1 (MQ=255)                                                                                      ccTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCa                                                              >  2:1319146/1‑150 (MQ=255)                                                                                        ttCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCAcc                                                            >  2:1132154/1‑150 (MQ=255)                                                                                        ttCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCAcc                                                            >  2:615454/1‑150 (MQ=255)                                                                                                cccACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTTCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCACCACCCAcaccac                                                    >  2:238992/1‑150 (MQ=255)                                                                                                cccACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGTTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCAcaccac                                                    <  1:625754/150‑1 (MQ=255)                                                                                                 ccACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACaccacc                                                   >  1:807043/1‑150 (MQ=255)                                                                                                  cACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACct                                                  <  1:1319146/150‑1 (MQ=255)                                                                                                   aCTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATATCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACctc                                                 >  2:755935/1‑150 (MQ=255)                                                                                                           cTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCt                                                                          >  1:527932/1‑117 (MQ=255)                                                                                                           cTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCt                                                                          <  2:527932/117‑1 (MQ=255)                                                                                                                          gggCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCcac                                                         >  1:908568/1‑119 (MQ=255)                                                                                                                          gggCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCcac                                                         <  2:908568/119‑1 (MQ=255)                                                                                                                                 tGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCtt                   >  1:292053/1‑150 (MQ=255)                                                                                                                                     gTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTcc               <  1:615454/150‑1 (MQ=255)                                                                                                                                     gTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTGCGGTAAGGCATCAGTGAACTGGCTTTCC‑‑CACCTTCACCACCGCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTcc               >  1:167869/1‑150 (MQ=255)                                                                                                                                      tATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTc                                                                                             <  2:135337/73‑1 (MQ=255)                                                                                                                                      tATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTc                                                                                             >  1:135337/1‑73 (MQ=255)                                                                                                                                      tATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCAcaccac                                                    >  1:34621/1‑112 (MQ=255)                                                                                                                                      tATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCAcaccac                                                    <  2:34621/112‑1 (MQ=255)                                                                                                                                      tATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTCCt              >  1:558621/1‑150 (MQ=255)                                                                                                                                       cgggggCTGCCAGGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTCCTg             >  2:1236201/6‑150 (MQ=255)                                                                                                                                         cAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTCCTGGc           <  1:1285222/150‑1 (MQ=255)                                                                                                                                          aaGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAg                                                                            >  2:1228121/1‑84 (MQ=255)                                                                                                                                          aaGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAg                                                                            <  1:1228121/84‑1 (MQ=255)                                                                                                                                           aGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACctcgctc                                             >  1:1009342/1‑114 (MQ=255)                                                                                                                                           aGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACctcgctc                                             <  2:1009342/114‑1 (MQ=255)                                                                                                                                                  cTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCAAGCAACGACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTCCTGGCCGTacaaca  >  1:1236800/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GGCATCCATGCACAAGGCGCAGGCTTGGGGGAGGTCGTCAGGGGTGAGCGTGCGTACTTGCATGGTAATCGTCCTTGTCGGGTGTCCTTCTGTTGCCCACTTTAACCTGTTACGGGTTCAAGGGCCGATGCGGTATCAAGCTGCCCTGGCAGCCCCCGGCTTGCGTTCGCAATCAGGATGTCTTCGAAGAATGCCCCGACAGGTTCGGTAGGGTGGCA‑‑CAGCTGTATCTCCAGCACCCACACCACCTCGCTCGGGACAATCTCCACCTCGGCCAGCTTGTCCTGGCCGTACAACA  >  NC_002947/3024778‑3025072

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