Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,657,956	G→A	7.2%	R4H (CGT→CAT)	sahR →	methionine metabolism transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,657,956	0	G	A	7.2%	155.5 / 2.2	55	R4H (CGT→CAT)	sahR	methionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (30/21);  new base A (2/2);  total (32/23)
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.60e-01

GAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCA  >  NC_002947/5657812‑5658100                                                                                                                                                  |                                                                                                                                                  gaagaaGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgcg                                                                                                                                               >  5:279575/1‑149 (MQ=255)  aagatgCAGAGCGCGACTCGGGTCTCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTTGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgc                                                                                                                                                <  6:40661/147‑1 (MQ=255)  aagaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgcg                                                                                                                                               >  3:7470/1‑148 (MQ=255)  aagaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgcg                                                                                                                                               >  2:257114/1‑148 (MQ=255)                  cTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGc                                                                                                                               >  3:273514/1‑148 (MQ=255)                   tCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCGTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGc                                                                                                                               >  1:78086/1‑148 (MQ=255)                   tCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCaa                                                                                                                             >  6:86499/1‑149 (MQ=255)                     ggggCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcg                                                                                                                           >  7:93486/1‑149 (MQ=255)                         cgcCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGaa                                                                                                                       <  4:115151/149‑1 (MQ=255)                          gcCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCACCGAGCAACGCGaaa                                                                                                                      >  1:270468/1‑149 (MQ=255)                            cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACa                                                                                                                    >  4:363374/1‑149 (MQ=255)                             ccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgc                                                                                                                                                <  6:32544/120‑1 (MQ=255)                             ccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgc                                                                                                                                                >  5:32544/1‑120 (MQ=255)                                         aCTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgt                                                                                                       <  5:172565/149‑1 (MQ=255)                                           ttGAGGCGTTCGGCTTGCTGC‑TCATCGCCAATATCAATAGTTTTTGTTTGTGGCCTTGCTAGGGCATCGATGCCTGTCTATACTGCCGCCCCATGTATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTTGCAGCCCTGTGc                                                                                                     <  4:221546/149‑1 (MQ=255)                                            tGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGc                                                                                                     >  2:102855/1‑148 (MQ=255)                                             gAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCa                                                                                                    >  7:21417/1‑148 (MQ=255)                                             gAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCATGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCa                                                                                                    >  3:348068/1‑148 (MQ=255)                                             gAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCATGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCa                                                                                                    >  3:353574/1‑148 (MQ=255)                                                gCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGc                                                                                                <  5:86499/149‑1 (MQ=255)                                                       gCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcccc                                                                                                                                   <  4:373293/107‑1 (MQ=255)                                                       gCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcccc                                                                                                                                   >  3:373293/1‑107 (MQ=255)                                                       gCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGc                                                                                         >  1:5566/1‑149 (MQ=255)                                                        cTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCg                                                                                        <  6:279575/149‑1 (MQ=255)                                                            cAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtg                                                                                     <  2:5566/148‑1 (MQ=255)                                                                  cATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCATGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgctgc                                                                               <  4:348068/149‑1 (MQ=255)                                                                   aTCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCATGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgctgc                                                                               <  4:353574/148‑1 (MQ=255)                                                                    tCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgctgcgc                                                                             <  4:7470/149‑1 (MQ=255)                                                                    tCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgctgcgc                                                                             >  7:29308/1‑149 (MQ=255)                                                                      gCCAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCTGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt                                                                           >  5:269273/1‑149 (MQ=255)                                                                          atatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAAc                                                                       <  3:363374/149‑1 (MQ=255)                                                                              cAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAt                                                                   >  5:308931/1‑149 (MQ=255)                                                                              cAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAt                                                                   >  8:176666/1‑149 (MQ=255)                                                                               aaaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAtt                                                                  >  2:314159/1‑149 (MQ=255)                                                                                aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg                                                                 >  2:18365/1‑149 (MQ=255)                                                                                aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg                                                                 >  6:303829/1‑149 (MQ=255)                                                                                aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg                                                                 >  5:245620/1‑149 (MQ=255)                                                                                aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg                                                                 >  1:249471/1‑149 (MQ=255)                                                                                          atatTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTg                                                       <  8:93486/148‑1 (MQ=255)                                                                                           tatTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGc                                                     <  8:29308/149‑1 (MQ=255)                                                                                                gCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGTAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAAt                                                 >  2:1684/1‑148 (MQ=255)                                                                                                 ccTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGa                                               >  4:7150/1‑149 (MQ=255)                                                                                                  cTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGAc                                              >  6:105414/1‑149 (MQ=255)                                                                                                      cTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGt                                          >  4:172784/1‑149 (MQ=255)                                                                                                         gggCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCg                                       <  7:176666/149‑1 (MQ=255)                                                                                                           gCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGc                                     <  4:273514/149‑1 (MQ=255)                                                                                                             aTCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGt                                   <  2:249471/149‑1 (MQ=255)                                                                                                                gATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGTAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCt                                >  7:103260/1‑149 (MQ=255)                                                                                                                gATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGc                                 <  6:245620/148‑1 (MQ=255)                                                                                                                 aTGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTg                               >  1:130377/1‑149 (MQ=255)                                                                                                                  tGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTgg                              <  8:21417/149‑1 (MQ=255)                                                                                                                            gACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCAGCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTgg                        <  1:257114/146‑1 (MQ=255)                                                                                                                            gACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGgc                       <  6:308931/146‑1 (MQ=255)                                                                                                                             aCTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAAAATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGgcg                      >  5:175987/1‑146 (MQ=255)                                                                                                                              cTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGa                  <  1:314159/149‑1 (MQ=255)                                                                                                                               tGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGgcg                      >  5:150256/1‑144 (MQ=255)                                                                                                                               tGCCGCCCAATGAATCTGCGTGCGCAATAGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑TCCTAGGCCAATGACTCGTACGGCGTGCTGGAACTGCCGCAGa                  >  5:194892/1‑148 (MQ=255)                                                                                                                                      ccATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGc                                                     >  7:189053/1‑106 (MQ=255)                                                                                                                                      ccATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGc                                                     <  8:189053/106‑1 (MQ=255)                                                                                                                                             tCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTc                                        >  7:221352/1‑112 (MQ=255)                                                                                                                                             tCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTc                                        <  8:221352/112‑1 (MQ=255)                                                                                                                                                 cgtgcgCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCa  <  1:102855/146‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGC‑TCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGC‑GCCTTGGCCAATGACTCGTTCGGCGTGCTGGAACTGGCGCAGATCTTCGACATCGGCCA  >  NC_002947/5657812‑5658100

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