Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,669,674	C→A	8.3%	intergenic (+41/+34)	prpD → / ← acnB	2‑methylcitrate dehydratase/aconitate hydratase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,669,674	0	C	A	8.3%	188.9 / 6.5	60	intergenic (+41/+34)	prpD/acnB	2‑methylcitrate dehydratase/aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base A (2/3);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669528‑2669819                                                                                                                                                      |                                                                                                                                                  ccGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTg                                                                                                                                                 <  7:2103/149‑1 (MQ=255)  cGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGc                                                                                                                                                <  7:66910/149‑1 (MQ=255)   gCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGcc                                                                                                                                               <  4:259790/149‑1 (MQ=255)       tGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCAGCGATGCCgg                                                                                                                                             <  1:241606/147‑1 (MQ=255)                   gtgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCTTAGGGACTTGTTTGTGATATGAAGCGGCCGGCAAAGGCCCAAACAA‑CAAAGCCCCCGCATCGCTGCCGGGGCTTTGTTTa                                                                                                                                  >  1:95135/1‑148 (MQ=255)                   gtgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACaa                                                                                                                               >  7:204933/1‑149 (MQ=255)                    tgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAAc                                                                                                                              <  8:8852/149‑1 (MQ=255)                         gcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGt                                                                                                                          >  3:139602/1‑148 (MQ=255)                                   gAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTa                                                                                                               >  8:32544/1‑149 (MQ=255)                                    aaCTGTGCAAGGACCAGCATAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCTGCCGGCACAGGACCAAACAA‑CAAAGCCCCGGCATCGCTGCCGCGGCTTTGTTTACAACCTGTAACTCAGCTTAc                                                                                                              >  4:27468/1‑149 (MQ=255)                                        gtgCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCt                                                                                                          >  1:189400/1‑149 (MQ=255)                                            aaGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAac                                                                                                      >  2:228999/1‑149 (MQ=255)                                                  cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCAGCgg                                                                                                >  4:160397/1‑149 (MQ=255)                                                  cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCgg                                                                                                >  5:233117/1‑149 (MQ=255)                                                  cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCgg                                                                                                >  4:160401/1‑149 (MQ=255)                                                    gcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGa                                                                                              >  1:136867/1‑149 (MQ=255)                                                        aaGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCtt                                                                                          <  2:189400/149‑1 (MQ=255)                                                              gAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGtt                                                                                    >  1:144109/1‑149 (MQ=255)                                                              gAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGt                                                                                     <  2:136867/148‑1 (MQ=255)                                                                    aCGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGc                                                                              >  6:276486/1‑149 (MQ=255)                                                                    aCGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGc                                                                              <  3:160397/149‑1 (MQ=255)                                                                    aCGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGc                                                                              <  3:160401/149‑1 (MQ=255)                                                                    aCGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGc                                                                              >  5:4070/1‑149 (MQ=255)                                                                            gCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcgc                                                                      >  1:112633/1‑149 (MQ=255)                                                                             cacaGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcgcg                                                                     <  2:95135/149‑1 (MQ=255)                                                                                   ttCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTc                                                               <  8:278760/149‑1 (MQ=255)                                                                                      gT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGG‑CCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcg                                                                       <  8:6372/137‑1 (MQ=255)                                                                                      gT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGG‑CCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCAGatcggaag                                                            >  7:6372/1‑141 (MQ=255)                                                                                           gACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGAt                                                         <  1:228999/149‑1 (MQ=255)                                                                                                gtttgtGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTgg                                                     <  2:144109/148‑1 (MQ=255)                                                                                                gtttgtGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                    >  6:72334/1‑149 (MQ=255)                                                                                                 tttgtGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTg                                                      <  2:68383/146‑1 (MQ=255)                                                                                                 tttgtGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGACATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTg                                                      >  1:68383/1‑146 (MQ=255)                                                                                                       atatGAAGCGGCCGGCACAGGCCCAAACAACCAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGc                                              >  1:26233/1‑149 (MQ=255)                                                                                                        tatGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTc                                            <  5:276486/149‑1 (MQ=255)                                                                                                          tGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAg                                          >  3:48391/1‑149 (MQ=255)                                                                                                           gAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAgg                                         >  5:303472/1‑149 (MQ=255)                                                                                                            aaGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAgg                                         >  1:61939/1‑148 (MQ=255)                                                                                                             aGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAgg                                         >  7:60897/1‑147 (MQ=255)                                                                                                               cggccggcACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGc                                     >  7:175946/1‑149 (MQ=255)                                                                                                                ggccggcACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCg                                    <  8:60897/149‑1 (MQ=255)                                                                                                                    ggcACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCGGGTAGCGGTaa                                <  2:112633/149‑1 (MQ=255)                                                                                                                       acaGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTTGCGGTAAACg                             >  5:306739/1‑149 (MQ=255)                                                                                                                        caGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGt                            <  7:32544/149‑1 (MQ=255)                                                                                                                                 taaa‑caaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                  <  5:154639/132‑1 (MQ=255)                                                                                                                                 taaa‑caaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                  >  6:154639/3‑134 (MQ=255)                                                                                                                                  caa‑caaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                               >  6:54857/1‑120 (MQ=255)                                                                                                                                   aa‑caaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                               >  1:172237/1‑119 (MQ=255)                                                                                                                                   aa‑caaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                               <  2:172237/119‑1 (MQ=255)                                                                                                                                   aa‑caaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCAGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATg                 >  3:219937/1‑149 (MQ=255)                                                                                                                                       aaaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTg              >  5:187725/1‑149 (MQ=255)                                                                                                                                        aaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGt             <  6:4070/149‑1 (MQ=255)                                                                                                                                            cccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGAt         <  3:27468/149‑1 (MQ=255)                                                                                                                                            cccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTc            <  7:275945/146‑1 (MQ=255)                                                                                                                                            cccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTc            >  8:275945/1‑146 (MQ=255)                                                                                                                                             ccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                  >  6:222406/1‑123 (MQ=255)                                                                                                                                             ccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                  <  5:222406/123‑1 (MQ=255)                                                                                                                                             ccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATg        >  5:190442/1‑149 (MQ=255)                                                                                                                                                   tCGCTGCCGGGTCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc  <  4:139602/149‑1 (MQ=255)                                                                                                                                                   tCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc  <  5:72334/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGT‑‑GGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAA‑CAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669528‑2669819

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: ‑