Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,489,782	G→T	7.3%	M88I (ATG→ATT)	PP_4825 →	MarC family membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,489,782	0	G	T	7.3%	210.6 / 3.7	66	M88I (ATG→ATT)	PP_4825	MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base G (24/37);  new base T (3/2);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489637‑5489924                                                                                                                                                   |                                                                                                                                               tGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGc                                                                                                                                              >  5:16673/1‑148 (MQ=255)   gCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGccc                                                                                                                                            <  2:345916/148‑1 (MQ=255)    ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCCCC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAg                                                                                                                                          <  5:33468/149‑1 (MQ=255)     cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         <  1:171159/149‑1 (MQ=255)      tGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        >  1:196629/1‑149 (MQ=255)      tGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        >  3:195568/1‑149 (MQ=255)      tGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         >  1:319054/1‑148 (MQ=255)        gAAGGTCGCGTTCGCAGCCTTTGTCGCCAGCGTGTTGCTGTACCTGTTCGGGGGGGCGATCTTCGGAATTTTCGTCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGc                                                                                                                                       <  8:340412/148‑1 (MQ=255)        gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        <  3:393814/147‑1 (MQ=255)        gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        <  3:374700/147‑1 (MQ=255)        gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         <  2:319054/146‑1 (MQ=255)         aaGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCaa                                                                                                                                     <  6:289581/149‑1 (MQ=255)            gTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGccgaaa                                                                                                                                  <  6:411228/149‑6 (MQ=255)              cgcgTTCGCCGCCTTGGTCGCCTGCGTGGTGCTGAACCTATTCGGGCGGGCGATCTTCGTAATTTTCGGCAGCACC‑GCCTATGCCTTCCGCCTCGGGGCTGGCAGCGTGGTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcgg                                                                                                                                <  6:236494/149‑1 (MQ=255)              cgcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcgg                                                                                                                                >  4:345555/1‑149 (MQ=255)               gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGccgaaa                                                                                                                                  >  5:2733/1‑141 (MQ=255)               gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAGGGccgaaa                                                                                                                                  <  6:2733/146‑6 (MQ=255)                cgTTCGCAGCCTTGGGCGCCAGCGGGTTGCCGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCGGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                              <  3:236352/149‑1 (MQ=255)                      cAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAg                                                                                                                        <  3:345555/149‑1 (MQ=255)                      aaGCCTTGGGCGACAGGGTTTTTCTGTACCTGTTCGGGCGGGCTATCTTGGTAATTTTCGGCATCACA‑GCGGATGCCTTCCGAAGCGGGGGTGGCAGCGTGCTGTTCATTTCGCCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAg                                                                                                                        <  6:16673/148‑1 (MQ=255)                       aGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAgg                                                                                                                       <  1:58066/149‑1 (MQ=255)                       aGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCGGGGCATGGCCCAGGGCAAGCCGGCGGTACAgg                                                                                                                       <  5:168721/149‑1 (MQ=255)                        gCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTAc                                                                                                                          >  3:325824/1‑145 (MQ=255)                        gCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTAc                                                                                                                          <  4:325824/145‑1 (MQ=255)                          ctgggTCGCCAGGGTGTTGCTGTACCTGTTCGGGCGGGCTATCTTCGGAATTTTCGGCATCACC‑GCCTATGCCTTCCGCATCGGTGGTGGCAGCGGGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAg                                                                                                                    <  4:195568/146‑1 (MQ=255)                              gTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATgg                                                                                                                                               <  8:10741/118‑1 (MQ=255)                              gTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATgg                                                                                                                                               >  7:10741/1‑118 (MQ=255)                               tCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCAGGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAAc                                                                                                               <  6:349905/149‑1 (MQ=255)                                 gCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACga                                                                                                             >  7:100877/1‑148 (MQ=255)                                    aGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTca                                                                                                          <  1:173231/149‑1 (MQ=255)                                       gtgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagca                                                                                                       <  6:290142/149‑1 (MQ=255)                                             cTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAtgt                                                                                                 >  1:148365/1‑149 (MQ=255)                                                tACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGAc                                                                                              <  4:290993/149‑1 (MQ=255)                                                 aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGAc                                                                                              <  5:322562/148‑1 (MQ=255)                                                  ccTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGAcc                                                                                             <  4:393940/148‑1 (MQ=255)                                                  ccTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCa                                                                                            >  2:228945/1‑149 (MQ=255)                                                   cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCAt                                                                                           >  7:99431/1‑149 (MQ=255)                                                      ttCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GACGATGCCTTCCGCATCGGTGCTGGCATCGTGCTGTTGATTTCGGCCCTGGGAATTGCCCAGGGCAAGCCGGCGGTCCCGGCAGCCAACGTTCAGAAGGATGTGACCAGCg                                                                                         >  4:213857/1‑148 (MQ=255)                                                           gcgggcgATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGc                                                                                   >  3:232748/1‑149 (MQ=255)                                                            cgggctATCTTCGGTATTTTCGG‑ATTACCGGCCGTTGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCTTTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGcc                                                                                  <  4:232748/149‑2 (MQ=255)                                                                 gATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCa                                                                             >  5:416526/1‑149 (MQ=255)                                                                      tCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCg                                                                        >  6:135968/1‑149 (MQ=255)                                                                        ggAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt                                                                                                   >  8:69374/1‑120 (MQ=255)                                                                        ggAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt                                                                                                   <  7:69374/120‑1 (MQ=255)                                                                        ggAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                      >  4:28678/1‑149 (MQ=255)                                                                         gAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCAtt                                                                           >  5:319533/1‑143 (MQ=255)                                                                         gAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCAtt                                                                           <  6:319533/143‑1 (MQ=255)                                                                         gAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                      >  1:124800/1‑148 (MQ=255)                                                                              ttCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGt                                                                >  3:272220/1‑149 (MQ=255)                                                                                              gATGCCTTCCGCATGGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCAt                                                 <  8:338903/149‑1 (MQ=255)                                                                                                 gCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCg                                               <  2:148365/148‑1 (MQ=255)                                                                                                    ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAcc                                                                                     >  6:373178/1‑107 (MQ=255)                                                                                                    ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAcc                                                                                     <  5:373178/107‑1 (MQ=255)                                                                                                        gCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctg                                       >  6:19160/1‑149 (MQ=255)                                                                                                               tGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTgg                                                             <  6:178515/120‑1 (MQ=255)                                                                                                               tGCTGGCAGCGTGCTGTTCATTACGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTgg                                                             >  5:178515/1‑120 (MQ=255)                                                                                                                     cAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATggggg                          <  6:187416/149‑1 (MQ=255)                                                                                                                     cAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATggggg                          <  6:187413/149‑1 (MQ=255)                                                                                                                             tGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagccag                  >  4:146270/1‑149 (MQ=255)                                                                                                                               ttCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagccagcc                >  6:30619/1‑149 (MQ=255)                                                                                                                                  ctttCGGCCCTGGGCATGTCCCAGTGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGa                                <  4:242134/129‑1 (MQ=255)                                                                                                                                     tCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACt           <  2:124800/148‑1 (MQ=255)                                                                                                                                       ggCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTggg        <  8:100877/149‑1 (MQ=255)                                                                                                                                       ggCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTgg         <  3:333419/148‑1 (MQ=255)                                                                                                                                        gCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGa       <  3:146270/149‑1 (MQ=255)                                                                                                                                          ccTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGtt                                                               >  8:285563/1‑91 (MQ=255)                                                                                                                                            tGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGa    >  1:131926/1‑148 (MQ=255)                                                                                                                                              ggCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaaa  <  4:272220/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               TGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACC‑GCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489637‑5489924

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