Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,681,789:1	+C	10.0%	coding (79/1431 nt)	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,789	1	.	C	10.0%	216.7 / 21.0	60	coding (79/1431 nt)	mltD	membrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base . (27/27);  new base C (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 = 6.39e-01

AAAGGTTGACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGA  >  NC_002947/4681641‑4681931                                                                                                                                                      |                                                                                                                                               aaaGGTTGACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAg                                                                                                                                                 >  3:228025/1‑149 (MQ=255)   aGGTTGACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑c                                                                                                                                               >  6:313289/1‑148 (MQ=255)     gTTGACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCAc                                                                                                                                            >  5:122258/1‑149 (MQ=255)           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                                                                                                      >  7:195419/1‑149 (MQ=255)                                           tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                      >  7:24035/1‑149 (MQ=255)                                           tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                      >  1:255835/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)                                                                                                                                                      |                                                                                                                                               AAAGGTTGACCAGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGA  >  NC_002947/4681641‑4681931

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