Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,681,796	T→A	10.0%	L29* (TTG→TAG) ‡	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,796	0	T	A	10.0%	167.7 / 8.8	60	L29* (TTG→TAG) ‡	mltD	membrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base A (3/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

CAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681651‑4681941                                                                                                                                                   |                                                                                                                                                  cAGGACCGGAGCGGTTTCTAGAATCGCAGAACTTTTTCGCCCGGATACATGTATCAGCCGTTGTCTTCCCGTAGCCGCAGAACCTCTCATTCGGTCGCCCTGTCGCGCCTGGCCCAAATCACTGCGCTGGCCCTGGCAG‑CCACCTtggt                                                                                                                                                <  6:100330/149‑1 (MQ=255)   tgCCCGGAGGGGTTTCTAGTATCGCGGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTtggtgg                                                                                                                                              <  6:211428/148‑1 (MQ=255)          aGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCa                                                                                                                                       <  5:291341/149‑1 (MQ=255)            cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCaga                                                                                                                                     <  2:163047/149‑1 (MQ=255)              gTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGT‑GGCTGCCAGGGc                                                                                                                                   <  7:7211/149‑1 (MQ=255)                    aGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGT‑GGCTGCCAGGGccag                                                                                                                                <  5:12573/146‑4 (MQ=38)                    aGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGT‑GGCTGCCAGGGccag                                                                                                                                >  6:12573/1‑143 (MQ=38)                    aGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGt                                                                                                                             <  5:209314/149‑1 (MQ=255)                     gAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCTCTGGCAGCCCACCAAGGT‑GGCTGCCAGGGccagc                                                                                                                               <  6:223279/146‑5 (MQ=38)                     gAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                               <  2:76005/146‑1 (MQ=255)                        tCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTGTCATTCCGTCGCACTGACGCGCCTGGCACAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGc                                                                                                                         >  5:204034/1‑149 (MQ=255)                             gAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    <  5:313289/149‑1 (MQ=255)                             gAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    <  4:85926/149‑1 (MQ=255)                              aaCTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                    >  4:61842/1‑148 (MQ=255)                                cTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgacga                                                                                                                  >  3:238146/1‑148 (MQ=255)                                 tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                                >  1:255835/1‑149 (MQ=255)                                 tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                                >  7:195419/1‑149 (MQ=255)                                 tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                                >  7:24035/1‑149 (MQ=255)                                   tttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                                <  2:255835/147‑1 (MQ=255)                                   tttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATc                                                                                                               <  5:146776/148‑1 (MQ=255)                                    ttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCATCCGTCAGCTCGACGAATcc                                                                                                              >  5:392893/1‑148 (MQ=255)                                      cGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGa                                                                                                            <  2:200420/148‑1 (MQ=255)                                       gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGa                                                                                                            >  6:136378/1‑147 (MQ=255)                                       gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGAt                                                                                                           <  4:238146/148‑1 (MQ=255)                                         ccGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGc                                                                                                        >  8:32932/1‑149 (MQ=255)                                             aTACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTc                                                                                                    <  4:228025/149‑1 (MQ=255)                                                 cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgcgc                                                                                                >  1:145218/1‑149 (MQ=255)                                                  tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                               <  1:276297/117‑1 (MQ=255)                                                  tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                               >  2:276297/1‑117 (MQ=255)                                                      tCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGcaca                                                                                           <  1:264398/149‑1 (MQ=255)                                                        aGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAAc                                                                                         <  5:332758/149‑1 (MQ=255)                                                         gCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAAc                                                                                         >  3:276200/1‑148 (MQ=255)                                                            gATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CACCCTTGGTGGGCTGCCAGAGCACCCGACAGCTCGACGAATCCGATAGCGTTCTCGCGCACAACTAca                                                                                     >  4:189293/1‑148 (MQ=255)                                                                  ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTGGACCAATCCGATAGCGTTCGCGCGAACAACTAGCAGGCAc                                                                               >  8:372983/1‑149 (MQ=255)                                                                     ccGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACgg                                                                             <  3:61842/148‑1 (MQ=255)                                                                       gTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATc                                                                          <  8:195419/149‑1 (MQ=255)                                                                        tAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCa                                                                         <  8:238883/149‑1 (MQ=255)                                                                          gCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCCCAACTACCAGGCACGGATCAAg                                                                       >  6:151528/1‑149 (MQ=255)                                                                           ccGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGc                                                                      <  5:51768/149‑1 (MQ=255)                                                                                   ccTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCt                                                              >  7:228128/1‑149 (MQ=255)                                                                                      ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                           >  5:215800/1‑149 (MQ=255)                                                                                      ctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcc                                                           >  6:284216/1‑149 (MQ=255)                                                                                         aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcccc                                                         <  6:122258/148‑1 (MQ=255)                                                                                         aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcccc                                                         <  5:136378/148‑1 (MQ=255)                                                                                         aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGACAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACACCTACCAGGCACGGATCAAGCACAAACCTTCCCCg                                                        >  4:133371/1‑149 (MQ=255)                                                                                                 gCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCCGCTCGACCAATCCCATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCgcagatgct                                                >  8:26623/1‑142 (MQ=255)                                                                                                         gcgcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccgg                                        <  6:204034/149‑1 (MQ=255)                                                                                                          cgcCTGGCCCAAATCCGTGCGCTGGCCCTCGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggc                                       <  8:135564/149‑1 (MQ=255)                                                                                                           gcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggcc                                      <  7:337192/149‑1 (MQ=255)                                                                                                              tGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGaa                                   >  6:367398/1‑149 (MQ=255)                                                                                                                  ccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                               >  6:167167/1‑149 (MQ=255)                                                                                                                    aaaTCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgcg                             >  3:346044/1‑149 (MQ=255)                                                                                                                     aaTCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGCGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGgcgc                            >  3:239835/1‑149 (MQ=255)                                                                                                                      aTCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGAGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgcc                           >  4:264468/1‑149 (MQ=255)                                                                                                                           tGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAg                       >  5:32214/1‑148 (MQ=255)                                                                                                                               cTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGt                  <  7:32932/149‑1 (MQ=255)                                                                                                                                    tcTGGCAGCCCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTggg             >  5:250616/2‑149 (MQ=255)                                                                                                                                    tcTGGCAGCCCACCAAGGT‑GGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTggg             >  3:381184/2‑149 (MQ=255)                                                                                                                                     cTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa            <  5:367398/149‑1 (MQ=255)                                                                                                                                              cACCTTGGTGGGCTGCCAGAGCATCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATTTCTGGGAGCGCATgc    <  6:392893/149‑1 (MQ=255)                                                                                                                                              cACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgc    >  2:63574/1‑149 (MQ=255)                                                                                                                                                ccTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATgcgc  <  2:145218/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGC  >  NC_002947/4681651‑4681941

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