Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,312,733	A→C	7.6%	intergenic (+31/‑131)	PP_2912 → / → hemB	hypothetical protein/5‑aminolevulinate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,312,733	0	A	C	7.6%	187.2 / 4.5	68	intergenic (+31/‑131)	PP_2912/hemB	hypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base C (3/2);  total (36/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 = 8.53e-01

GGGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCCCCT  >  NC_002947/3312584‑3312882                                                                                                                                                       |                                                                                                                                                       gggTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTa                                                                                                                                                        >  2:1926593/1‑150 (MQ=255)         cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGaa                                                                                                                                               <  1:1658721/150‑1 (MQ=255)         cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGaa                                                                                                                                               <  1:213432/150‑1 (MQ=255)         cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGaa                                                                                                                                               >  2:722288/1‑150 (MQ=255)          gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAag                                                                                                                                              <  2:2388704/150‑1 (MQ=255)           tGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAaga                                                                                                                                             <  1:705507/150‑1 (MQ=255)                  gACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAgg                                                                                                                                           <  1:711342/145‑1 (MQ=255)                  gACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAgg                                                                                                                                           >  2:711342/1‑145 (MQ=255)                      ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAg                                                                                                                                  <  2:1262004/150‑1 (MQ=255)                      ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAg                                                                                                                                  <  2:2056059/150‑1 (MQ=255)                       cgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAAGAGGcggt                                                                                                                                       >  2:575004/1‑141 (MQ=255)                       cgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAAGAGGcggt                                                                                                                                       <  1:575004/144‑4 (MQ=255)                          ccGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAaga                                                                                                                                             >  1:969934/1‑135 (MQ=255)                          ccGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGC‑GGGTAACCCCGCGAaga                                                                                                                                             <  2:969934/135‑1 (MQ=255)                                ggACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                        <  1:134286/150‑1 (MQ=255)                                ggACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                        <  1:1926593/150‑1 (MQ=255)                                    gacgCTCGTCACCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCg                                                                                                                    >  1:1336893/1‑150 (MQ=255)                                       gCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCCACACAGGCAAACCACGGCCGTTg                                                                                                                 >  1:125929/1‑150 (MQ=255)                                       gCTCGTCGCCGTGCAGCAC‑GCGCCGAAGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTg                                                                                                                 >  1:1644269/1‑150 (MQ=255)                                         tcgtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAAcaca                                                                                                                                   >  1:1097568/1‑130 (MQ=255)                                         tcgtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAAcaca                                                                                                                                   <  2:1097568/130‑1 (MQ=255)                                            ggcgcGTGCTGCACGGCGACGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGGGGCCAACACAGGCAAACCACGGCCGTTGCCAt                                                                                                             >  2:1160120/5‑150 (MQ=255)                                                             gccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtc                                                                                            >  2:1588423/1‑150 (MQ=255)                                                               cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                          <  1:2237302/150‑1 (MQ=255)                                                               cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGC‑GGGTAACCCCGCGAAGAgg                                                                                                                                           <  2:1376750/101‑1 (MQ=38)                                                               cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGATGCAGGAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAAGAgg                                                                                                                                           >  1:1376750/1‑101 (MQ=38)                                                                gaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCg                                                                                         >  1:345211/1‑150 (MQ=255)                                                                 aggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGt                                                                                        >  1:792059/1‑150 (MQ=255)                                                                  ggccgagCATGGGCCGGCCAAGGCAGAGGAGGGAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCACCACAGTCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc                                                                                       <  1:1160120/150‑1 (MQ=255)                                                                     cgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCAtt                                                                                    <  1:1588423/150‑1 (MQ=255)                                                                      gagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCAttt                                                                                   <  2:1336893/150‑1 (MQ=255)                                                                              gccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtttt                                                                           >  1:2236513/1‑150 (MQ=255)                                                                                cggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                         >  2:153335/1‑118 (MQ=255)                                                                                cggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                         <  1:153335/118‑1 (MQ=255)                                                                                  gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACt                                                                       <  2:125929/150‑1 (MQ=255)                                                                                   ccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTg                                                                      <  2:2236513/150‑1 (MQ=255)                                                                                    cAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                                  <  1:2059876/137‑1 (MQ=255)                                                                                    cAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                                  >  2:2059876/1‑137 (MQ=255)                                                                                       ggCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCAc                                                                        <  2:2178358/144‑1 (MQ=255)                                                                                       ggCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCAGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTAACGCAAGACACATTCTCGTCATTTAGCCTTTTCAc                                                                        >  1:2178358/1‑144 (MQ=255)                                                                                                 gTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc                                                                                       <  2:385094/119‑1 (MQ=255)                                                                                                 gTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTc                                                                                       >  1:385094/1‑119 (MQ=255)                                                                                                 gTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGAAAGACACATTCTCGTCATTTAGCCTTGTCACTGGCTCCCGACAGGCa                                                        >  2:590595/1‑150 (MQ=255)                                                                                                     cgccgcCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                                    <  2:345211/150‑1 (MQ=255)                                                                                                       ccgccCCGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTTGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                                  <  1:20964/150‑1 (MQ=255)                                                                                                         gccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGtat                                                >  1:1156743/1‑150 (MQ=255)                                                                                                           ccTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATAc                                              >  1:520544/1‑150 (MQ=255)                                                                                                            cTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACg                                             >  1:2279267/1‑150 (MQ=255)                                                                                                               gcgcgagtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGAAACATAATAGTACATGAGCCTATTCCAGGGAGCTCGACGGGCATAGGGTGTAGGCCg                                          >  1:1667276/6‑150 (MQ=255)                                                                                                                 tGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                                  >  2:256869/1‑108 (MQ=255)                                                                                                                 tGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                                  <  1:256869/108‑1 (MQ=255)                                                                                                                   cAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCa                                      >  1:351969/1‑150 (MQ=255)                                                                                                                      tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAAt                                   >  1:960919/1‑150 (MQ=255)                                                                                                                      tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAAAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAAt                                   >  2:1677822/1‑150 (MQ=255)                                                                                                                         tttATTCCTGTACCGCCTCTTCGCGGGGTTG‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGTCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                            >  2:1514984/1‑122 (MQ=255)                                                                                                                         tttATTCCTGTACCGCCTCTTCGCGGGGTTG‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGTCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                            <  1:1514984/122‑1 (MQ=255)                                                                                                                            aTTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTg                             <  2:960919/150‑1 (MQ=255)                                                                                                                                cTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGc                         >  2:2062492/1‑150 (MQ=255)                                                                                                                                cTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGc                         >  2:1443550/1‑150 (MQ=255)                                                                                                                                 tGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                        <  2:1284235/150‑1 (MQ=255)                                                                                                                                   tACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt                      <  1:722288/150‑1 (MQ=255)                                                                                                                                     ccGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCa                    >  1:1522431/1‑150 (MQ=255)                                                                                                                                         ctctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTc                >  1:1169475/1‑150 (MQ=255)                                                                                                                                         ctctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTc                <  2:2277379/150‑1 (MQ=255)                                                                                                                                            ttCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa             >  1:1349703/1‑150 (MQ=255)                                                                                                                                               gcgGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacca          >  2:218769/1‑150 (MQ=255)                                                                                                                                               gcgGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacca          >  2:1663484/1‑150 (MQ=255)                                                                                                                                               gcgGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaacca          >  2:1284847/1‑150 (MQ=255)                                                                                                                                                  gggTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAAtt       <  2:792059/150‑1 (MQ=255)                                                                                                                                                       a‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCCCCt  <  1:1284847/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       GGGTGATGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCCCCT  >  NC_002947/3312584‑3312882

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