Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,807,984	A→G	7.9%	intergenic (+119/+36)	mrcA → / ← maeB	penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,807,984	0	A	G	7.9%	168.5 / 5.1	63	intergenic (+119/+36)	mrcA/maeB	penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base G (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

GAAGCGACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGA  >  NC_002947/5807839‑5808132                                                                                                                                                  |                                                                                                                                                      gAAGCGACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGc                                                                                                                                                    >  2:93445/1‑149 (MQ=255)         gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCggg                                                                                                                                            >  3:262628/1‑149 (MQ=255)          ccGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTTGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCgggg                                                                                                                                           <  2:232799/149‑1 (MQ=255)          ccGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg                                                                                                                                           <  2:158331/149‑1 (MQ=255)          ccGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg                                                                                                                                           <  2:153070/149‑1 (MQ=255)             aTGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑Ctt                                                                                                                                       >  4:265833/1‑149 (MQ=255)                 aTCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑Ctttttt                                                                                                                                   <  1:93445/149‑1 (MQ=255)                    tGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGa                                                                                                                                                  <  5:8863/132‑1 (MQ=255)                    tGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCAACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGa                                                                                                                                                  >  6:8863/1‑132 (MQ=255)                      ttCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTc                                                                                                                              >  3:318933/1‑149 (MQ=255)                           aGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTtg                                                                                                                         >  2:305142/1‑149 (MQ=255)                            ggCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg                                                                                                                                           <  3:33043/131‑1 (MQ=255)                            ggCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg                                                                                                                                           >  4:33043/1‑131 (MQ=255)                            ggCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTtgc                                                                                                                        >  5:18536/1‑149 (MQ=255)                                          cTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGtt                                                                                                          >  3:98322/1‑149 (MQ=255)                                          cTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGtt                                                                                                          >  4:68900/1‑149 (MQ=255)                                           tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGtt                                                                                                          >  4:47502/1‑148 (MQ=255)                                           tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGtt                                                                                                          >  1:147839/1‑148 (MQ=255)                                           tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTCACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCTCTTGCTGAAATCAGCCttt                                                                                                          >  8:160857/1‑148 (MQ=255)                                             aGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAg                                                                                                               <  6:242057/141‑1 (MQ=255)                                             aGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAg                                                                                                               >  5:242057/1‑141 (MQ=255)                                             aGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTg                                                                                                         <  4:318933/147‑1 (MQ=255)                                                ttttGCGTGGGCGCTCGATCTTATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAACGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAaca                                                                                                     <  7:160857/148‑1 (MQ=255)                                                ttttGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCGACATAAAAAAGCCCCGGCTCGCATGAGCCGGGG‑c                                                                                                                                         <  3:171792/112‑1 (MQ=37)                                                ttttGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCGACATAAAAAAGCCCCGGCTCGCATGAGCCGGGG‑c                                                                                                                                         >  4:171792/1‑112 (MQ=37)                                                           cgcTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGc                                                                                         <  6:18536/149‑1 (MQ=255)                                                                      aTAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGc                                                                                 <  7:12886/146‑1 (MQ=255)                                                                      aTAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGc                                                                                 <  3:68900/146‑1 (MQ=255)                                                                      aTAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGc                                                                                 <  3:265833/146‑1 (MQ=255)                                                                          gcgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAAt                                                                          >  3:242203/1‑149 (MQ=255)                                                                               acaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTc                                                                     <  2:147839/149‑1 (MQ=255)                                                                                 aaGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCAt                                                                                               >  1:45593/1‑121 (MQ=255)                                                                                 aaGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCAt                                                                                               <  2:45593/121‑1 (MQ=255)                                                                                      cccTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAg                                                              <  4:242203/149‑1 (MQ=255)                                                                                      cccTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAg                                                              <  1:305142/149‑1 (MQ=255)                                                                                       ccTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAgg                                                             >  5:93875/1‑149 (MQ=255)                                                                                           aCCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGt                                                                                                           >  4:103637/1‑99 (MQ=255)                                                                                           aCCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAgg                                                         <  4:202760/149‑1 (MQ=255)                                                                                             cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGAGTTTTTTATTTCGATTGATGCAATCCGCCGTTGAACAGTACATCCACGATGGTCAGCGGGTATTGATTCGGATAAGGAAggg                                                        >  1:105141/1‑149 (MQ=255)                                                                                                         gCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGa                                           >  1:191152/1‑149 (MQ=255)                                                                                                          cccTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAt                                          >  7:184783/1‑149 (MQ=255)                                                                                                          cccTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAt                                          >  2:184846/1‑149 (MQ=255)                                                                                                           ccTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAtt                                         <  8:184783/149‑1 (MQ=255)                                                                                                            cTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTc                                        >  3:253081/1‑149 (MQ=255)                                                                                                               acacATTCCCAAGCCCCATAAAAACGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGa                                      <  4:262628/148‑1 (MQ=255)                                                                                                               acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGa                                      <  3:47502/148‑1 (MQ=255)                                                                                                               acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGAt                                     <  6:93875/149‑1 (MQ=255)                                                                                                                   aTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAg                                                          >  4:64734/1‑124 (MQ=255)                                                                                                                   aTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAg                                                          <  3:64734/124‑1 (MQ=255)                                                                                                                       ccAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTc                                        >  5:270592/1‑138 (MQ=255)                                                                                                                       ccAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTc                                        <  6:270592/138‑1 (MQ=255)                                                                                                                       ccAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCt                             >  6:226716/1‑149 (MQ=255)                                                                                                                        cAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCtt                            >  7:265122/1‑149 (MQ=255)                                                                                                                         aaGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTg                           <  4:253081/149‑1 (MQ=255)                                                                                                                              ccATAAAAAAGCCCCGTCTCATGCGAGACGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTaa                                                                           <  1:146719/96‑1 (MQ=255)                                                                                                                              ccATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTTCAATCAGCCGTTGAACAATTCATCCACGCTGGTCACCGGGTaa                                                                           >  2:146719/1‑96 (MQ=255)                                                                                                                              ccATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCcac                      >  1:332922/1‑149 (MQ=255)                                                                                                                                aTAAAAAAGCCCCGGCTCGCATGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAtt                                         <  6:305650/128‑1 (MQ=255)                                                                                                                                aTAAAAAAGCCCCGGCTCGCATGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAtt                                         >  5:305650/1‑128 (MQ=255)                                                                                                                                aTAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCtt                            >  6:178119/1‑141 (MQ=255)                                                                                                                                aTAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCtt                            <  5:178119/141‑1 (MQ=255)                                                                                                                                         ccccGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGaaa           <  4:98322/149‑1 (MQ=255)                                                                                                                                              gCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCATCCTTGGCCACAGCATCGGAAacgacg      <  1:184846/149‑1 (MQ=255)                                                                                                                                                  aTGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCAGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGa  >  8:36350/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GAAGCGACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGG‑CTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGA  >  NC_002947/5807839‑5808132

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