Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,312,687	G→T	15.7%	intergenic (+28/‑134)	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,687	0	G	T	15.7%	136.9 / 16.4	57	intergenic (+28/‑134)	PP_2912/hemB	hypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base T (2/7);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATT  >  NC_002947/3312550‑3312834                                                                                                                                           |                                                                                                                                                       gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                             <  1:133629/149‑1 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                             <  8:247839/149‑1 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                             <  3:364859/149‑1 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCcgcg                                                                                                                                               >  2:303536/1‑147 (MQ=255) gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGaa                                                                                                                                             <  4:28116/149‑1 (MQ=255)   gagaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAaga                                                                                                                                           <  7:282010/149‑1 (MQ=255)    agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAagag                                                                                                                                          <  1:199908/149‑1 (MQ=255)     gaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAgg                                                                                                                                         >  2:58871/1‑149 (MQ=255)        cGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                             >  8:248989/1‑142 (MQ=255)        cGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                             <  7:248989/142‑1 (MQ=255)          aCGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGAAGAGGcg                                                                                                                                       <  4:158971/146‑2 (MQ=255)          aCGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTACCCCGCGAAGAGGcggc                                                                                                                                     >  3:158971/1‑145 (MQ=255)              cgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGt                                                                                                                                                       <  6:319976/126‑1 (MQ=255)              cgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGt                                                                                                                                                       >  5:319976/1‑126 (MQ=255)              cgccCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAg                                                                                                                                >  2:320221/1‑149 (MQ=255)                 ccGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca                                                                                                                                 <  3:307722/145‑1 (MQ=255)                 ccGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCAGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAAcaca                                                                                                                                 >  4:307722/1‑145 (MQ=255)                        gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                             <  6:313517/126‑1 (MQ=255)                        gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAACCCCGCGaa                                                                                                                                             >  5:313517/1‑126 (MQ=255)                        gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                      <  4:319359/149‑1 (MQ=255)                        gacgacGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                      >  3:22030/1‑149 (MQ=255)                           gacgCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑cc                                                                                                                >  5:77153/1‑149 (MQ=255)                                tcgtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGc                                                                                                           >  5:150513/1‑149 (MQ=255)                                ctcggcgcGTGCTGCACGGCGACGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTTCCCGCGCAGAGGCCAACACCGGCAAACCCCCG‑‑‑CCGTTg                                                                                                            >  2:167603/8‑149 (MQ=255)                                  gtcgCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCa                                                                                                         <  1:320221/149‑1 (MQ=255)                                    cgccgcGCAGCAC‑GCGCCGAGGCCGAGCAAGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCAGGCTGCAGTAGTTTATTCCTGTACCCCCTCTTCGCGGTGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCAtt                                                                                                       <  2:121483/149‑1 (MQ=255)                                       cGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGAc                                                                                                    <  1:58871/149‑1 (MQ=255)                                             gcaC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGa                                                                                              >  3:224476/1‑149 (MQ=255)                                             gcaC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGa                                                                                              <  5:225418/149‑1 (MQ=255)                                               aC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCATCGCGGGGTTTCCCGCGAAGAGGCCAACACAGGCAAACCCCGG‑‑‑CCGTTGCCATTGACGCAAcaca                                                                                            >  3:355766/1‑149 (MQ=255)                                                c‑gcgCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTACCCGCGAAGAGGGCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGaca                                                                                            >  3:47698/1‑148 (MQ=255)                                                      cgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTG‑CATTGACGCAAGACACATtctc                                                                                     <  4:129861/149‑1 (MQ=255)                                                       gaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATtctc                                                                                     <  7:297287/149‑1 (MQ=255)                                                             gagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCAtt                                                                               <  4:224476/149‑1 (MQ=255)                                                                      ccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtttt                                                                      <  4:22030/149‑1 (MQ=255)                                                                         gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTAGCGGGGTTACCCGCGAAGCGGCCAACACACGCCAACCCCGG‑‑‑CCGTTGGCATTGACGCAAGACACATTCTAGTCATTTAGCCttttgcc                                                                   >  6:215773/1‑146 (MQ=255)                                                                               gCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTACCCGCCAAGAGGCCCACACAGGCAAACCACGG‑‑‑CCGTTGCCAATGCCGCAAGACACCTTGTCGTCATTTAGCCTTTTCACTGgggc                                                             >  4:249580/1‑149 (MQ=255)                                                                               gCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTAGCGGGGTTTACCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCCTTGCCATTGACGCAAGACCCATTATAGGGATTTATCCTTTTAACTGgcgg                                                             >  2:152995/1‑148 (MQ=255)                                                                                     gaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGa                                                         >  8:53253/1‑147 (MQ=255)                                                                                     gaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                       <  1:167603/149‑1 (MQ=255)                                                                                             gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                               <  3:249580/149‑1 (MQ=255)                                                                                             gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                               <  4:198717/149‑1 (MQ=255)                                                                                             gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                               <  6:150513/149‑1 (MQ=255)                                                                                                gccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACATGCAAACCACGGGACCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAg                                               <  1:303536/149‑1 (MQ=255)                                                                                                                  tATTCCTGTACCGCCGATTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTcc                           <  2:86224/148‑1 (MQ=255)                                                                                                                    ttCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCt                         <  4:47698/148‑1 (MQ=255)                                                                                                                      ggatctccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGa                      <  5:129266/143‑1 (MQ=255)                                                                                                                       cTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAg                     >  2:261972/1‑149 (MQ=255)                                                                                                                       cTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGGTGCCATTGACGCAAGACACATTCTCCTCATTTAGCCTTTTCACTGGCGCCCGACAGGCATAAGGTATAAGCCGTTCAAATTCCCTg                        >  1:66974/1‑146 (MQ=255)                                                                                                                        tGTACCGCCTCTTCCCGGGGTTACCCGCGAAGAGCCCAACACAGGCAAACCACGG‑‑‑CCGTTTCCAATGACCCAAGACACAATCTCGTCATTTAGACTTTTCACTGGCGAGCCACAAGCATAAGGTATAAGCCGATCAAATTCCATgtac                     >  8:327621/1‑145 (MQ=255)                                                                                                                         gTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCt                   >  5:313066/1‑149 (MQ=255)                                                                                                                            ccGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc                <  1:261972/149‑1 (MQ=255)                                                                                                                             cGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCAGTTCAAATTCCCTGGAGCTGTCa               >  4:238526/1‑149 (MQ=255)                                                                                                                              gCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt                 <  2:123486/146‑1 (MQ=255)                                                                                                                              gCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGAGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGt                 >  1:123486/1‑146 (MQ=255)                                                                                                                                ctctTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGt            <  1:152995/149‑1 (MQ=255)                                                                                                                                      gcgGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaac       <  5:215773/148‑1 (MQ=255)                                                                                                                                          ggTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGGCAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAAtt  >  1:330448/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                       GTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCAC‑GCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTACCCGCGAAGAGGCCAACACAGGCAAACCACGG‑‑‑CCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATT  >  NC_002947/3312550‑3312834

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