Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,681,795	T→A	5.3%	L29M (TTG→ATG) ‡	mltD →	membrane‑bound lytic murein transglycosylase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,681,795	0	T	A	5.3%	240.0 / 2.7	76	L29M (TTG→ATG) ‡	mltD	membrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base T (43/29);  new base A (2/2);  total (45/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 = 5.25e-01

CCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681650‑4681941                                                                                                                                                    |                                                                                                                                                   ccAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTtgg                                                                                                                                                 <  2:68930/149‑1 (MQ=255)    ggCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTtggtgg                                                                                                                                              >  4:273920/1‑149 (MQ=255)       ccGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCt                                                                                                                                           <  8:124602/149‑1 (MQ=255)       ccGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCt                                                                                                                                           <  2:50983/149‑1 (MQ=255)        cGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCGCAATCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTg                                                                                                                                          <  8:181685/149‑1 (MQ=255)               gTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCagag                                                                                                                                    >  3:278016/1‑148 (MQ=255)                    tAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCg                                                                                                                              >  1:116247/1‑149 (MQ=255)                       aaTCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                             >  5:180964/1‑147 (MQ=255)                       aaTCGCCGAACTTTTTCGCCCCGATACCTGTATCAGCCCAAGACTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGACCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCAGCCAGACCACCCGTCa                                                                                                                           >  6:282822/1‑149 (MQ=255)                             cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGTGCTGCCAGAGCACCCGTCAGCTcgc                                                                                                                     >  8:301107/1‑148 (MQ=255)                             cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                             >  4:149280/1‑141 (MQ=255)                             cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                             <  3:149280/141‑1 (MQ=255)                               aaCTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    >  1:170874/1‑148 (MQ=255)                                 cTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGGACCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    >  1:114192/1‑148 (MQ=255)                                        gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCGGAGCACCCGTCAGCTCGACGAATCCGa                                                                                                            >  2:70004/1‑147 (MQ=255)                                        gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGa                                                                                                            >  7:162010/1‑147 (MQ=255)                                        gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATa                                                                                                          <  2:174461/149‑1 (MQ=255)                                         cccGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAg                                                                                                         <  3:129768/149‑1 (MQ=255)                                              aTACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTc                                                                                                    <  5:310931/149‑1 (MQ=255)                                                aCCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                                  >  2:119371/1‑149 (MQ=255)                                                aCCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                                  >  5:321979/1‑149 (MQ=255)                                                 ccTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcg                                                                                                 >  4:34298/1‑149 (MQ=255)                                                  cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                >  6:344940/1‑149 (MQ=255)                                                  cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                                  >  4:61377/1‑147 (MQ=255)                                                  cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                                  <  3:61377/147‑1 (MQ=255)                                                   tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgcg                                                                                               >  8:259659/1‑149 (MQ=255)                                                   tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgcg                                                                                               >  1:207620/1‑149 (MQ=255)                                                   tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCAGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgcg                                                                                               >  6:96154/1‑149 (MQ=255)                                                       tgaGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCTCCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcac                                                                                            <  5:282822/146‑1 (MQ=255)                                                       tCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcac                                                                                            >  3:160561/1‑148 (MQ=255)                                                           ccGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTa                                                                                       >  5:41918/1‑149 (MQ=255)                                                            cGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAc                                                                                      >  5:39815/1‑149 (MQ=255)                                                            cGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGCCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAc                                                                                      <  5:344940/149‑1 (MQ=255)                                                                  cTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCa                                                                                >  6:216740/1‑149 (MQ=255)                                                                     gccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCa                                                                                <  5:182276/145‑1 (MQ=255)                                                                       cGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAt                                                                           <  6:39815/149‑1 (MQ=255)                                                                                  aaCCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAc                                                                >  1:60439/1‑149 (MQ=255)                                                                                      tctcATTCCGTCTCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGATCGACGAAACCGATACCGTTCGCGCGCACAACTACCACCCACGGATCAAGCACAAACCTTc                                                            >  6:284996/1‑149 (MQ=255)                                                                                      tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                            >  4:3927/1‑149 (MQ=255)                                                                                       ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                            >  4:219047/1‑148 (MQ=255)                                                                                       ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                            <  3:219047/148‑1 (MQ=255)                                                                                       ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTA‑CAGGCACGGATCAAGCACAAACCTTccc                                                          >  1:299320/1‑149 (MQ=255)                                                                                        tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTccc                                                          >  5:259303/1‑149 (MQ=255)                                                                                          aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCg                                                        >  2:38662/1‑149 (MQ=255)                                                                                          aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCACTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCg                                                        >  6:320900/1‑149 (MQ=255)                                                                                              cGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgctgc                                                    >  4:175457/1‑149 (MQ=255)                                                                                              cGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgctgc                                                    >  7:243920/1‑149 (MQ=255)                                                                                                 cGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGa                                                                                                            <  8:91007/90‑1 (MQ=255)                                                                                                    ccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTc                                               >  2:91215/1‑148 (MQ=255)                                                                                                      tGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACaaa                                                                 <  8:127433/128‑1 (MQ=255)                                                                                                      tGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACaaa                                                                 >  7:127433/1‑128 (MQ=255)                                                                                                            gcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggcc                                      <  3:175457/149‑1 (MQ=255)                                                                                                             gcTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggcc                                      >  4:25410/2‑148 (MQ=255)                                                                                                             ccTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggccg                                     <  6:321979/149‑1 (MQ=255)                                                                                                                 gCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGaa                                   <  2:116247/147‑1 (MQ=255)                                                                                                                 gCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGaa                                   >  1:216034/1‑147 (MQ=255)                                                                                                                 gCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGaa                                   <  3:34298/147‑1 (MQ=255)                                                                                                                 gCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACa                                 <  4:160561/149‑1 (MQ=255)                                                                                                                 gCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACa                                 <  3:25410/149‑1 (MQ=255)                                                                                                                  cccAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAg                                <  4:278016/149‑1 (MQ=255)                                                                                                                        tCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgccg                          <  6:180964/149‑1 (MQ=255)                                                                                                                          aGTGCGCTGGACCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACACGTCAGCTCGAGGAATCCGATAGCGTTCGCGCGCACAACTACCATGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgccg                          >  4:107348/1‑147 (MQ=255)                                                                                                                                   gCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTc                 <  7:259659/147‑1 (MQ=255)                                                                                                                                    ctcTGGCAGC‑CCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACaaa                                                                 <  4:104536/96‑1 (MQ=37)                                                                                                                                    ctcTGGCAGC‑CCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACaaa                                                                 >  3:104536/3‑98 (MQ=37)                                                                                                                                    cccTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAATAGGCGCCGCAGGATGTCTgg              <  1:119371/149‑1 (MQ=255)                                                                                                                                    cccTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              <  3:3927/149‑1 (MQ=255)                                                                                                                                    cccTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              <  6:41918/149‑1 (MQ=255)                                                                                                                                      cTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTgg              >  7:1887/1‑147 (MQ=255)                                                                                                                                       tGGCAGC‑CCACCAAGGT‑GGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCa                                                                                    <  3:217004/76‑1 (MQ=25)                                                                                                                                       tGGCAGC‑CCACCAAGGT‑GGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCa                                                                                    >  4:217004/1‑76 (MQ=25)                                                                                                                                               ccACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATg     >  4:243130/1‑149 (MQ=255)                                                                                                                                               ccACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATg     >  5:318274/1‑149 (MQ=255)                                                                                                                                                 aCCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcg   >  7:45460/1‑149 (MQ=255)                                                                                                                                                  ccTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCGGCTGGTTGTCTGGGAGCGCATgcgc  >  3:196163/1‑149 (MQ=255)                                                                                                                                                  ccTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcgc  <  5:96154/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681650‑4681941

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