Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,971,954	G→T	8.8%	intergenic (‑389/+28)	PP_4379 ← / ← flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,971,954	0	G	T	8.8%	179.2 / 7.0	57	intergenic (‑389/+28)	PP_4379/flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base T (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

CGGGCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTA  >  NC_002947/4971806‑4972093                                                                                                                                                     |                                                                                                                                             cgggaACTGGAGCCCGAATGATCAGCTTCGCAGCAGTTATCGATGTGTCCGCATCGCGGCCCTTTCGATACTCAAGGCCGCTCACACGAATACGCCACCTGCAGTAGAGATGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTtcg                                                                                                                                              <  4:51650/149‑1 (MQ=255)    gcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                           <  5:263851/149‑1 (MQ=255)    gcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                           >  6:263851/1‑149 (MQ=255)          gACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAG‑ccc                                                                                                                                    >  8:71839/1‑149 (MQ=255)            ccccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑ccc                                                                                                                                    >  7:40181/1‑147 (MQ=255)             cccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAg                                                                                                                                 >  4:211568/1‑149 (MQ=255)              ccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAG‑ccctatc                                                                                                                                <  7:71839/149‑5 (MQ=255)              ccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAg                                                                                                                                 >  3:99610/1‑148 (MQ=255)               cGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAG‑ccctatcg                                                                                                                               <  8:239704/149‑6 (MQ=255)                gAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGAt                                                                                                                              <  7:168401/149‑1 (MQ=255)                 aaTGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATc                                                                                                                             <  5:356372/149‑1 (MQ=255)                   tGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACAACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcag                                                                                                                           <  6:255305/149‑1 (MQ=255)                      tCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggca                                                                                                                        <  5:3511/149‑1 (MQ=255)                        aGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcagg                                                                                                                      >  5:50342/1‑149 (MQ=255)                        aGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcag                                                                                                                       <  3:211568/148‑1 (MQ=255)                         gCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcag                                                                                                                       <  1:23018/147‑1 (MQ=255)                         gCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcag                                                                                                                       >  2:23018/1‑147 (MQ=255)                              gcagcaGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGAt                                                                                                                 >  6:2024/1‑148 (MQ=255)                              gcagcaGATATCGATGGGGCCGCATCACGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGAt                                                                                                                 <  2:324829/148‑1 (MQ=255)                                 gcagATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGa                                                                                                             <  8:40181/149‑1 (MQ=255)                                        tCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGt                                                                                                       <  5:2024/148‑1 (MQ=255)                                                  cGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCa                                                                                                                                  <  2:15678/111‑1 (MQ=255)                                                  cGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCt                                                                                             >  7:200885/1‑148 (MQ=255)                                                  cGCATCGCGGCCCTTTCGCGACGCAAGGACGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCa                                                                                                                                  >  1:15678/1‑111 (MQ=255)                                                   gCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTa                                                                                           <  3:70496/149‑1 (MQ=255)                                                    cATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAg                                                                                          <  4:99610/149‑1 (MQ=255)                                                        gcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAg                                                                                          <  1:248474/146‑1 (MQ=255)                                                        gcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAg                                                                                          >  2:248474/1‑145 (MQ=255)                                                            cccTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTgcg                                                                                                                                              <  6:65927/90‑1 (MQ=255)                                                            cccTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTgcg                                                                                                                                              >  5:65927/1‑90 (MQ=255)                                                            cccTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGaaa                                                                                  >  6:175276/1‑149 (MQ=255)                                                              cTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            <  3:55387/105‑1 (MQ=255)                                                              cTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            >  4:55387/1‑105 (MQ=255)                                                              cTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATg                                                                                >  5:182102/1‑149 (MQ=255)                                                                  cgcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGTGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            >  6:187700/1‑101 (MQ=255)                                                                  cgcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGTGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            <  5:187700/101‑1 (MQ=255)                                                                   gcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                            >  1:147793/1‑148 (MQ=255)                                                                           aGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGc                                                                     >  6:295793/1‑147 (MQ=255)                                                                                   cccACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGaa                                                           <  6:182102/149‑1 (MQ=255)                                                                                    ccACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCg                                                             >  1:74770/1‑146 (MQ=255)                                                                                    ccACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCg                                                             <  2:74770/146‑1 (MQ=255)                                                                                     cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGTCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAgg                                                                                               <  5:216631/111‑1 (MQ=255)                                                                                     cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGTCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAgg                                                                                               >  6:216631/1‑111 (MQ=255)                                                                                      aCGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                            <  5:275952/129‑1 (MQ=255)                                                                                      aCGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                            >  6:275952/1‑129 (MQ=255)                                                                                              gCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCAt                                                <  7:194187/149‑1 (MQ=255)                                                                                                   cTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCa                                                                   <  3:149532/125‑1 (MQ=255)                                                                                                   cTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCa                                                                   >  4:149532/1‑125 (MQ=255)                                                                                                      cAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctg                                          <  5:175276/147‑1 (MQ=255)                                                                                                                     ccTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGaa                         <  7:319218/149‑1 (MQ=255)                                                                                                                     ccTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGaa                         >  8:18980/1‑149 (MQ=255)                                                                                                                      cTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                    <  7:73608/121‑1 (MQ=255)                                                                                                                      cTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                    >  8:73608/1‑121 (MQ=255)                                                                                                                        gtgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAg                                                                                                              <  3:63537/61‑1 (MQ=255)                                                                                                                        gtgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAg                                                                                                              >  4:63537/1‑61 (MQ=255)                                                                                                                                cTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAg              >  2:314002/1‑149 (MQ=255)                                                                                                                                cTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAg              >  6:111856/1‑149 (MQ=255)                                                                                                                                    tgtCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGcc           <  6:50342/148‑1 (MQ=255)                                                                                                                                         cgATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCa     <  7:18980/149‑1 (MQ=255)                                                                                                                                            tggggCTGCTGCGCAG‑CCCTAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTa  >  7:32222/3‑149 (MQ=255)                                                                                                                                            tggggCTGCTGCGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTa  >  5:66232/3‑149 (MQ=255)                                                                                                                                                  tGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctg                                          >  2:48259/1‑103 (MQ=255)                                                                                                                                                     |                                                                                                                                             CGGGCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTA  >  NC_002947/4971806‑4972093

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑