Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,681,810	A→G	8.2%	S34G (AGC→GGC)	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,810	0	A	G	8.2%	164.9 / 5.2	62	S34G (AGC→GGC)	mltD	membrane‑bound lytic murein transglycosylase D
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base G (2/3);  total (31/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 = 4.73e-01

CGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTCGCCCTGC  >  NC_002947/4681662‑4681957                                                                                                                                                      |                                                                                                                                                    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                                                                                                                                    <  4:85926/149‑1 (MQ=255)                  gAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcgac                                                                                                                                    <  5:313289/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                                                                                                                                               >  2:276297/1‑117 (MQ=255)                                       tGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCAccc                                                                                                                                               <  1:276297/117‑1 (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                                                                         <  5:136378/148‑1 (MQ=255)                                                                              aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTcccc                                                                         <  6:122258/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)                                                                                                                                         ggggggCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCa                <  7:372983/144‑1 (MQ=255)                                                                                                                                          gtggGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAgg              <  5:284216/148‑1 (MQ=255)                                                                                                                                          gtggGCTGCCAGAGCACCCGGCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAgg              <  8:90921/148‑1 (MQ=255)                                                                                                                                                  ccAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTCGccc     >  3:280981/1‑149 (MQ=255)                                                                                                                                                     gagCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTCGCCCTGc  <  3:264468/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGCGCATGCGCCAGGGCTTCGCCCTGC  >  NC_002947/4681662‑4681957

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