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

AGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGC  >  NC_002947/4681652‑4681933                                                                                                                                                  |                                                                                                                                          aGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTtggtg                                                                                                                                       <  2:325566/149‑1 (MQ=255) aGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTtggtg                                                                                                                                       >  2:43859/1‑149 (MQ=255)  ggCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTtggtgg                                                                                                                                      <  5:105195/149‑1 (MQ=255)   gCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTggg                                                                                                                                     <  7:187159/149‑1 (MQ=255)     ccGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCt                                                                                                                                   >  2:208683/1‑149 (MQ=255)       ggAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGc                                                                                                                                 <  5:231204/149‑1 (MQ=255)       ggAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGc                                                                                                                                 >  7:99586/1‑149 (MQ=255)        gAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGcc                                                                                                                                <  4:280482/149‑1 (MQ=255)         aGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTTTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCGCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCa                                                                                                                               <  5:341328/149‑1 (MQ=255)         aGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGACCCTCGCAG‑CCACCTAGGTGGGCTGCCa                                                                                                                               >  8:293549/1‑149 (MQ=255)          gCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCag                                                                                                                              <  1:209830/149‑1 (MQ=255)           cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCaga                                                                                                                             <  7:150473/149‑1 (MQ=255)           cGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCaga                                                                                                                             <  8:99586/149‑1 (MQ=255)              tttCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCa                                                                                                                          <  3:95460/149‑1 (MQ=255)                    gAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTc                                                                                                                    <  1:43859/149‑1 (MQ=255)                          ccGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcg                                                                                                              >  4:253992/1‑149 (MQ=255)                           cGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTcga                                                                                                             >  1:141493/1‑149 (MQ=255)                             aaCTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGT‑GGCTGCCa                                                                                                                               <  6:19415/129‑1 (MQ=255)                             aaCTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAGCCCACCAAGGT‑GGCTGCCa                                                                                                                               >  5:19415/1‑129 (MQ=255)                                tttttCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAAt                                                                                                        >  2:269138/1‑149 (MQ=255)                                      gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGAt                                                                                                   <  2:34790/148‑1 (MQ=255)                                      gCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGAt                                                                                                   <  4:271893/148‑1 (MQ=255)                                               ccTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCACTGCGCTGACCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                          <  7:293549/148‑1 (MQ=255)                                                cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCa                                                                                                                               <  8:300505/110‑1 (MQ=255)                                                cTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCAGCCa                                                                                                                               >  7:300505/1‑110 (MQ=255)                                                     tCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCaga                                                                                                                             >  4:289057/1‑107 (MQ=255)                                                       aGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAAc                                                                                 >  4:140757/1‑149 (MQ=255)                                                           gATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTAcc                                                                             <  3:253992/149‑1 (MQ=255)                                                               tCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAgg                                                                          >  7:91246/1‑148 (MQ=255)                                                                 ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                       >  1:309365/1‑149 (MQ=255)                                                                 ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                       >  4:11310/1‑149 (MQ=255)                                                                 ttCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCAc                                                                       >  3:154834/1‑149 (MQ=255)                                                                     cGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGAt                                                                   >  5:307358/1‑149 (MQ=255)                                                                         gCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCaa                                                                <  3:140757/148‑1 (MQ=255)                                                                         gCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAg                                                               <  1:269138/149‑1 (MQ=255)                                                                         cgatcttaaCCTCTCATTCCGTCGCCCTGACTCGCCTGGCCCAAATCAGTGCGCTGGCACTGGCAG‑CCACCTTGGTGTTCTGCCAGAGCACCCTTCAGCTCGACGAATCCGATAGCGTTCTCGCGCACAACTACCAGGCACGGATCaa                                                                <  1:39143/141‑1 (MQ=255)                                                                            gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTc                                                                                                                    >  3:34717/1‑93 (MQ=255)                                                                            gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTc                                                                                                                    <  4:34717/93‑1 (MQ=255)                                                                            gCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGca                                                             >  1:15747/1‑148 (MQ=255)                                                                                aaCCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAc                                                        >  7:316730/1‑149 (MQ=255)                                                                                aaCCTCTCATTCCGTCGCCCTGACGCGCCTGGACCAAATCAGTGCGCTGGCAATGGCAG‑CCACATTGGTGGGCTGACAGAGCACCCGTCAGCTCGACGAATCAGATAGCCTTCGCGCCCACAACTACCAGGCACGGATaaa                                                                >  2:39143/1‑141 (MQ=255)                                                                                aaCCTCTCATTCCGTCGCCCTAACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAc                                                        >  8:6680/1‑149 (MQ=255)                                                                                 aCCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAAcc                                                       >  1:274346/1‑149 (MQ=255)                                                                                    tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                    >  6:17590/1‑149 (MQ=255)                                                                                    tctcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTc                                                    >  7:265093/1‑149 (MQ=255)                                                                                      tcATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTccc                                                  >  5:260328/1‑149 (MQ=255)                                                                                        aTTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTccc                                                  <  6:307358/147‑1 (MQ=255)                                                                                                  ccTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTAGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACACACCTTCCCCGCTGCTGGTc                                       >  1:254745/1‑148 (MQ=255)                                                                                                    tGACGCGCCTTGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTcgc                                                                                          <  4:339204/95‑1 (MQ=255)                                                                                                          gcCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAgccggcc                              <  6:260328/149‑1 (MQ=255)                                                                                                               gCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGaa                           >  2:41472/1‑147 (MQ=255)                                                                                                                 ccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                       >  1:112471/1‑149 (MQ=255)                                                                                                                 ccAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGAAGAATCCGATAGCGTTCGCGCGCACAACTACCAGGAACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAgg                       >  7:223672/1‑149 (MQ=255)                                                                                                                    aaTCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGctgctg                                          >  2:280451/1‑127 (MQ=255)                                                                                                                          tGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAg               >  7:240176/1‑148 (MQ=255)                                                                                                                           gcgcTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGa             >  8:233025/1‑149 (MQ=255)                                                                                                                            cgcTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGa             <  2:141493/148‑1 (MQ=255)                                                                                                                                 gCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCt        <  5:17590/148‑1 (MQ=255)                                                                                                                                 gCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCt        <  2:112471/148‑1 (MQ=255)                                                                                                                                 gCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTg       >  8:9617/1‑149 (MQ=255)                                                                                                                                    cTGGCAGCCCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATc                                                                  <  3:173784/87‑1 (MQ=25)                                                                                                                                    cTGGCAGCCCACCAAGGT‑GGCTGCCAGGGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATc                                                                  >  4:173784/1‑87 (MQ=25)                                                                                                                                    cTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGa    <  2:309365/149‑1 (MQ=255)                                                                                                                                       gCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgc  <  1:208683/148‑1 (MQ=255)                                                                                                                                       gCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAgc  <  8:265093/148‑1 (MQ=255)                                                                                                                                              aCCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgccgc                 >  4:76795/1‑127 (MQ=255)                                                                                                                                              aCCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGcgccgc                 <  3:76795/127‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                          AGGCCCGGAGCGGTTTCTAGAATCGCCGAACTTTTTCGCCCGGATACCTGTATCAGCCGATGTCTTCCCGTAGCCGCAGAACCTCTCATTCCGTCGCCCTGACGCGCCTGGCCCAAATCAGTGCGCTGGCCCTGGCAG‑CCACCTTGGTGGGCTGCCAGAGCACCCGTCAGCTCGACGAATCCGATAGCGTTCGCGCGCACAACTACCAGGCACGGATCAAGCACAAACCTTCCCCGCTGCTGGTCAAGCCGGCCGAACAGGCGCCGCAGGATGTCTGGGAGC  >  NC_002947/4681652‑4681933

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