Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,376,322	G→C	6.4%	intergenic (+31/+18)	kup → / ← PP_1201	potassium transport system protein Kup/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,376,322	0	G	C	6.4%	197.2 / 3.9	63	intergenic (+31/+18)	kup/PP_1201	potassium transport system protein Kup/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base C (2/2);  total (33/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 = 9.00e-01

TCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTG  >  NC_002947/1376175‑1376467                                                                                                                                                     |                                                                                                                                                  tCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGAAGGCTCAGGGCCGCaga                                                                                                                                                  <  6:329149/149‑1 (MQ=255)         gCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                          <  2:396046/149‑1 (MQ=255)           ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        <  3:311822/149‑1 (MQ=255)           ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                          >  3:271114/1‑147 (MQ=255)           ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  6:1726/149‑1 (MQ=255)            gggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAg                                                                                                                                       <  3:58941/149‑1 (MQ=255)            gggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAg                                                                                                                                       >  2:310122/1‑149 (MQ=255)             ggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGc                                                                                                                                      <  3:259370/149‑1 (MQ=255)               cAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCt                                                                                                                                    >  7:13887/1‑149 (MQ=255)                aaCCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGACGCAGAGCGGCCATGACCCtt                                                                                                                                   >  8:370521/1‑149 (MQ=255)                 aCCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        >  3:147459/1‑143 (MQ=255)                 aCCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        <  4:147459/143‑1 (MQ=255)                           gTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGCGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTg                                                                                                                         <  1:310122/149‑1 (MQ=255)                               ttGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAg                                                                                                                    >  1:168698/1‑149 (MQ=255)                                 gCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGc                                                                                                                   >  7:94181/1‑148 (MQ=255)                                    gaagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACg                                                                                                               >  2:257069/1‑149 (MQ=255)                                             cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTg                                                                                                      >  6:370956/1‑149 (MQ=255)                                             cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTg                                                                                                      >  7:236311/1‑149 (MQ=255)                                             cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTg                                                                                                      <  5:95975/149‑1 (MQ=255)                                              aaCGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCa                                                                                                                               <  2:282327/123‑1 (MQ=255)                                              aaCGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCa                                                                                                                               >  1:282327/1‑123 (MQ=255)                                              aaCGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGc                                                                                                     >  3:337621/1‑149 (MQ=255)                                               aCGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGccc                                                                                                                                           >  2:258852/1‑110 (MQ=255)                                               aCGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGccc                                                                                                                                           <  1:258852/110‑1 (MQ=255)                                                   cAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTgg                                                                                                <  2:168698/149‑1 (MQ=255)                                                       ttGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGcgcg                                                                                            >  4:95617/1‑149 (MQ=255)                                                                  cAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtctc                                                                                  <  3:95617/148‑1 (MQ=255)                                                                  cAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtctc                                                                                  <  4:271114/148‑1 (MQ=255)                                                                  cAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtc                                                                                    <  2:331584/146‑1 (MQ=255)                                                                     ccTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGcc                                                                               >  7:198063/1‑148 (MQ=255)                                                                         ccGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGATATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATc                                                                           <  7:370521/148‑1 (MQ=255)                                                                            ttGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        >  4:11014/1‑84 (MQ=255)                                                                            ttGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  3:11014/84‑1 (MQ=255)                                                                                  cGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTgg                                                                 <  1:257069/149‑1 (MQ=255)                                                                                  cGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTgg                                                                 >  2:329512/1‑149 (MQ=255)                                                                                   gggTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGc                                                                <  8:13887/149‑1 (MQ=255)                                                                                      tGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCACGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAg                                                             <  6:286075/149‑1 (MQ=255)                                                                                           gAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCtt                                                                                                  <  7:187467/107‑1 (MQ=255)                                                                                           gAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCtt                                                                                                  >  8:187467/1‑107 (MQ=255)                                                                                           gAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCg                                                        >  5:332633/1‑149 (MQ=255)                                                                                            aGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                       >  7:64351/1‑149 (MQ=255)                                                                                            aGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                       <  8:94181/149‑1 (MQ=255)                                                                                            aGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                       >  3:120456/1‑149 (MQ=255)                                                                                                  ggACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGtt                                                 <  4:173088/149‑1 (MQ=255)                                                                                                   gACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTAAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACCCTGGACTTGCCTTGGCGCGTTTCGATCTAGCCAATCAACCCCTTGGCGAGTGCCGGGTAGTtc                                                >  3:118222/1‑149 (MQ=255)                                                                                                     cccAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                              >  7:174122/1‑149 (MQ=255)                                                                                                         ggTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgtcg                                          <  8:64351/149‑1 (MQ=255)                                                                                                         ggTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgtcg                                          <  8:174122/149‑1 (MQ=255)                                                                                                            cGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAg                                       >  3:333904/1‑149 (MQ=255)                                                                                                             gagaTCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAt                                                                                     <  5:4771/102‑1 (MQ=255)                                                                                                             gagaTCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAt                                                                                     >  6:4771/1‑102 (MQ=255)                                                                                                                     aaTG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtggtg                                  <  5:370956/145‑1 (MQ=255)                                                                                                                             ttATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCa                       <  4:333904/149‑1 (MQ=255)                                                                                                                                 aCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGtt                    >  5:369348/1‑148 (MQ=255)                                                                                                                                  cagcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                  <  8:198063/149‑1 (MQ=255)                                                                                                                                   agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGttt                   >  7:100506/1‑147 (MQ=255)                                                                                                                                   agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATATAGCCAATAAACCGCTTGGCGACTGCCGGGTAGTTATAGTCGAAGTGGTGCCCGCCACGCAGTTTc                  >  8:237068/1‑148 (MQ=255)                                                                                                                                      agGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgc              >  2:256738/1‑149 (MQ=255)                                                                                                                                        gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgcgc            >  5:24511/1‑149 (MQ=255)                                                                                                                                         cTCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt           >  6:202223/1‑149 (MQ=255)                                                                                                                                              ggCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCAcc       >  3:340713/1‑148 (MQ=255)                                                                                                                                                ccGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg                <  5:346376/137‑1 (MQ=255)                                                                                                                                                ccGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg                >  6:346376/1‑137 (MQ=255)                                                                                                                                                  gCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTg  <  8:100506/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATG‑GGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTG  >  NC_002947/1376175‑1376467

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