Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,850,158	C→T	9.2%	intergenic (+75/‑112)	hemN → / → fnrA	oxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,850,158	0	C	T	9.2%	183.0 / 7.0	65	intergenic (+75/‑112)	hemN/fnrA	oxygen‑independent coproporphyrinogen III dehydrogenase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base T (3/3);  total (34/31)
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.20e-01

GGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACA  >  NC_002947/4850017‑4850298                                                                                                                                               |                                                                                                                                              ggTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGc                                                                                                                                        <  2:280479/149‑1 (MQ=255)     tGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCgg                                                                                                                                     <  3:346117/148‑1 (MQ=255)       cATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                   <  2:200576/148‑2 (MQ=255)       cATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGTGGCCTGCCggga                                                                                                                                   <  4:274168/148‑2 (MQ=255)       cATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCgggg                                                                                                                                   <  7:212310/148‑1 (MQ=255)          ggTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTc                                                                                                                               >  4:141748/1‑149 (MQ=255)                    gcgGATCTGGCGACGCCCAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTTTGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAgcgcgc                                                                                                                    <  2:169212/149‑1 (MQ=255)                          cTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACcat                                                                                                               >  1:330864/1‑148 (MQ=255)                             gCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACcat                                                                                                               <  2:217179/145‑1 (MQ=255)                               gATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATcgc                                                                                                         >  5:29754/1‑149 (MQ=255)                               gATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATcgc                                                                                                         >  2:270433/1‑149 (MQ=255)                                aTGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATcgcg                                                                                                        <  2:31087/149‑1 (MQ=255)                                    acaACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGcaca                                                                                                    >  6:347862/1‑149 (MQ=255)                                    acaACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGcaca                                                                                                    >  4:184564/1‑149 (MQ=255)                                         cGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGc                                                                                               <  3:184564/149‑1 (MQ=255)                                          gCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGcc                                                                                              >  4:193877/1‑149 (MQ=255)                                          gCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGcc                                                                                              >  7:54879/1‑149 (MQ=255)                                           ccAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCt                                                                                             <  3:141748/149‑1 (MQ=255)                                             aGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTa                                                                                            >  8:197114/1‑148 (MQ=255)                                                     tcGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGGACACAATATGAGCCAACAAGGAAACCCGGCAGGCCACATCCCGGGGTTACG‑CCACCGCGCACCAACATCGGGCACAATGGCCTACACCAtga                                                                                    >  1:101475/1‑148 (MQ=255)                                                     tcGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCtgtg                                                                                   >  3:250988/1‑149 (MQ=255)                                                         gAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGt                                                                               <  5:154169/149‑1 (MQ=255)                                                         gAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCACACCATCATCGCGCACACTGGCCTACACCCTGTGTCGt                                                                               >  2:304520/1‑149 (MQ=255)                                                          aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGcaca                                                                                                    >  5:222149/1‑127 (MQ=255)                                                          aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGcaca                                                                                                    <  6:222149/127‑1 (MQ=255)                                                          aGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTg                                                                              <  8:54879/149‑1 (MQ=255)                                                           gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                             <  2:330864/149‑1 (MQ=255)                                                           gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                             <  5:347862/149‑1 (MQ=255)                                                           gTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGa                                                                             <  7:202100/149‑1 (MQ=255)                                                            tGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CGAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCTTGa                                                                             <  2:6175/148‑1 (MQ=255)                                                              aTCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGt                                                                           <  7:197114/148‑1 (MQ=255)                                                                      aaCAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCACGCCACCTGCCGGGGTTTCG‑CCAGCGCTCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTAc                                                                  >  1:19649/1‑149 (MQ=255)                                                                        cAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGa                                                                <  6:29754/149‑1 (MQ=255)                                                                        cAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGa                                                                >  2:303159/1‑149 (MQ=255)                                                                                   aCCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGc                                                     >  3:275426/1‑149 (MQ=255)                                                                                       cATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTc                                                 <  1:77910/149‑1 (MQ=255)                                                                                           gACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGAGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCcac                                             <  2:101475/149‑1 (MQ=255)                                                                                            aCAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAt                                                           >  6:369267/1‑134 (MQ=255)                                                                                            aCAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTAt                                                           <  5:369267/134‑1 (MQ=255)                                                                                             cAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCcac                                             >  3:19390/1‑147 (MQ=255)                                                                                             cAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACaa                                           <  5:129417/149‑1 (MQ=255)                                                                                              aaGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAg                                          >  8:32775/1‑149 (MQ=255)                                                                                                  ttCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGCCTTATGTGTGCTTTCCCACAAGGATc                                      >  8:186357/1‑149 (MQ=255)                                                                                                    ccTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCg                                     <  1:303159/148‑1 (MQ=255)                                                                                                      ttGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCGACCAGCGCGCACCATCATAGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTTTGCTTTCCCACAAGGATCGAt                                   >  1:285755/1‑149 (MQ=255)                                                                                                        gTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGAtt                                  <  7:32775/147‑1 (MQ=255)                                                                                                        gTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTc                                 <  3:193877/148‑1 (MQ=255)                                                                                                         tCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCa                                >  8:262013/1‑148 (MQ=255)                                                                                                         tCAGC‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTa                                                                   <  2:199233/113‑1 (MQ=255)                                                                                                            gg‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATg                            >  5:239599/1‑149 (MQ=255)                                                                                                                  acTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTcc                         >  1:50214/1‑147 (MQ=255)                                                                                                                         gCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTc                <  6:239599/149‑1 (MQ=255)                                                                                                                         gCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGCCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCtatt                                                 >  6:283813/1‑113 (MQ=255)                                                                                                                            cacaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCaaa             <  4:250988/149‑1 (MQ=255)                                                                                                                             acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAAc            >  6:44071/1‑149 (MQ=255)                                                                                                                             acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAAc            >  5:210738/1‑149 (MQ=255)                                                                                                                             acaAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAAc            >  1:245635/1‑149 (MQ=255)                                                                                                                              caAGGACTCCCGGCAGGCCAGCTGCCGGGGTTTCG‑CCAGCGCGCACCATCTTTCTGCACTCTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGAGTTATGTGTGCTTTCCCTCAAGGATCGATTCAAATGTCCGAGCCAGTCAAACt           <  1:302171/149‑1 (MQ=255)                                                                                                                                  aaaccccGGCAGGTGGCCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCAACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTgcgc       >  3:313459/5‑149 (MQ=255)                                                                                                                                     tCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCca    >  6:288182/1‑149 (MQ=255)                                                                                                                                      cccGGCAGGTGGCCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCca    >  5:40634/1‑148 (MQ=255)                                                                                                                                      cccGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcac   >  8:151558/1‑149 (MQ=255)                                                                                                                                       ccGGCAGGTGGCCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCcaca  >  2:309665/1‑149 (MQ=255)                                                                                                                                         ggCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGcc      >  7:64542/1‑143 (MQ=255)                                                                                                                                         ggCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGcc      <  8:64542/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              GGTATGCATGGTGTTCGACGCGTATCTGGCGATGCACAACCGCCAGCGTTTCTCGCGAGTGATCTGAAAAACAAGGCGTTCAACCGCATGGACAAGCTTCCTTGTCAGG‑CACACTATGAGCACACAAGGACTCCCGGCAGGCCACCTGCCGGGGTTTCG‑CCAGCGCGCACCATCATCGCGCACACTGGCCTACACCCTGTGTCGTGAGTTACCCTTACGACTTATGTGTGCTTTCCCACAAGGATCGATTCAAATGTCCGAGCCAGTCAAACTGCGCCCACA  >  NC_002947/4850017‑4850298

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