Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,914,141	T→C	9.7%	intergenic (+33/+120)	spuI → / ← argA	glutamylpolyamine synthetase/amino acid N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,914,141	0	T	C	9.7%	161.7 / 6.4	61	intergenic (+33/+120)	spuI/argA	glutamylpolyamine synthetase/amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base C (4/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAAT  >  NC_002947/5913993‑5914281                                                                                                                                                      |                                                                                                                                             aGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGttttt                                                                                                                                              <  8:197374/149‑1 (MQ=255)     cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                          <  5:80155/149‑1 (MQ=255)     cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                          >  5:153606/1‑149 (MQ=255)     cATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgtt                                                                                                                                           <  7:188933/148‑1 (MQ=255)        cGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                          >  3:289443/1‑146 (MQ=255)         gATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCTTGAGGGGCGGCGTTttcgtttaca                                                                                                                                       <  3:29135/148‑4 (MQ=39)         gATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTgaga                                                                                                                                      <  1:157957/149‑1 (MQ=255)           tCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGAcc                                                                                                                                    >  6:144072/1‑149 (MQ=255)             cGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCg                                                                                                                                   <  5:1483/148‑1 (MQ=255)                 gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                          <  6:190124/137‑1 (MQ=255)                 gTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTtttgttt                                                                                                                                          >  5:190124/1‑137 (MQ=255)                      tCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCct                                                                                                                         <  5:144072/149‑1 (MQ=255)                         aCATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTAGCAACTGGTACATGAATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCTTTTTTGTTTGAGACCGGTGTCGCCctct                                                                                                                       >  1:275325/1‑149 (MQ=255)                                           aGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTcc                                                                                                    >  4:274216/1‑149 (MQ=255)                                           aGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTcc                                                                                                    >  2:188387/1‑149 (MQ=255)                                            ggAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTccc                                                                                                   >  7:131724/1‑149 (MQ=255)                                            ggAAAGCGAGCTGGAAGAGTTCGAGCACTAGATCTCCGAACTCGAGTA‑CAACTGGTACATGCATACCGAGTAAAAGAAAACGCCGCCCCTCAAGGCCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCTCCCACGCCCGCAcc                                                                                                    >  8:200147/1‑148 (MQ=255)                                              aaaTCGATCTGGAAGTGTTCGATCACTCGATATCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGAGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCg                                                                                                  <  7:200147/148‑1 (MQ=255)                                               aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgc                                                                                                                   >  7:207381/1‑130 (MQ=255)                                               aaGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgc                                                                                                                   <  8:207381/130‑1 (MQ=255)                                                 gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                                              >  5:149639/1‑149 (MQ=255)                                                 gCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAgg                                                                                              <  6:197707/149‑1 (MQ=255)                                                        ggAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAAttt                                                                                       >  7:198001/1‑149 (MQ=255)                                                         gAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGccc                                                                                                         >  5:54526/1‑130 (MQ=255)                                                         gAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGccc                                                                                                         <  6:54526/130‑1 (MQ=255)                                                           agagTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCg                                                                                     >  4:284977/1‑148 (MQ=255)                                                              gTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAgg                                                                                 <  6:153606/149‑1 (MQ=255)                                                                  gatcACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGATTACGGCGCCCCTGTAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCGTCTTCGCGGGCACGCCCGCTCCCGc                                                                                                 <  4:96858/126‑1 (MQ=255)                                                                  gAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACAGTCTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTCTTTGAGACCGGTGTCGCCCTCTACGCGGGCACGCCCGCTCCCGc                                                                                                 >  3:96858/1‑129 (MQ=255)                                                                    gCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAg                                                                           >  3:205030/1‑149 (MQ=255)                                                                     cACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGc                                                                          <  2:217878/149‑1 (MQ=255)                                                                       cTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgcg                                                                                                                  <  4:28498/107‑1 (MQ=255)                                                                       cTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTcgcg                                                                                                                  >  3:28498/1‑107 (MQ=255)                                                                       cTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCt                                                                        >  4:61787/1‑149 (MQ=255)                                                                        tCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCACCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCGGAAGCCtg                                                                       >  7:232610/1‑149 (MQ=255)                                                                          gATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                      <  4:205030/148‑1 (MQ=255)                                                                               ccGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                                    <  3:270999/145‑1 (MQ=255)                                                                               ccGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                                    >  4:270999/1‑145 (MQ=255)                                                                                cGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGa                                                                    <  3:8724/144‑1 (MQ=255)                                                                                cGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGAc                                                                   <  8:131724/145‑1 (MQ=255)                                                                                cGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACtt                                                                 <  6:121695/147‑1 (MQ=255)                                                                                    cTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                                           <  3:274216/149‑1 (MQ=255)                                                                                    cTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgt                                                           <  6:149639/149‑1 (MQ=255)                                                                                     tCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                                          >  6:341728/1‑149 (MQ=255)                                                                                     tCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCtgtg                                                          >  7:249462/1‑149 (MQ=255)                                                                                             cAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCCCTACCGCAGGGTAATTTGGAAGGCTGAAGCCGGTGTCTTACCGGTGGGAGCgg                                                   >  7:244713/1‑149 (MQ=255)                                                                                              aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTGGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCggg                                                  >  5:139030/1‑149 (MQ=255)                                                                                              aaCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGAGGGAGCggg                                                  <  7:16623/149‑1 (MQ=255)                                                                                               aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                                                 <  4:316255/149‑1 (MQ=255)                                                                                               aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                                                 <  7:5378/149‑1 (MQ=255)                                                                                               aCTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                                                 <  4:289443/149‑1 (MQ=255)                                                                                                   gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                      >  4:28562/1‑124 (MQ=255)                                                                                                   gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCtgt                                                                      <  3:28562/124‑1 (MQ=255)                                                                                                   gTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGt                                               <  1:188387/147‑1 (MQ=255)                                                                                                        tGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCcgc                                         >  1:263344/1‑148 (MQ=255)                                                                                                                   caaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                                                 <  4:234908/128‑1 (MQ=38)                                                                                                                   caaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGc                                                 >  3:234908/2‑129 (MQ=38)                                                                                                                    aaaCAAAAACGCCGCCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTTGGAGCGGGCGTGCCCGCGATGACGGCGGcc                            >  6:89534/1‑149 (MQ=38)                                                                                                                        gAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCGGTGAGGTACCTGTGGGAGCGGGCGTGCCCGCGTAGAGGGCGGCCCAgg                        >  8:61722/1‑149 (MQ=255)                                                                                                                        gAAAACGCCGCCCCTCAAGGGCAGAGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCACCCGCAGGGGAAATTAGAAGGGGGAAGGCGGTGACATACCGGGGGGAGCGGGGGTGCCCGCGGAGAGGGGGGCCCAgg                        >  7:52185/1‑149 (MQ=255)                                                                                                                           aaCGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCt                      >  5:146569/1‑148 (MQ=255)                                                                                                                           aaCGCCGCCCCCCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCt                      >  2:58476/1‑148 (MQ=255)                                                                                                                             cgccgcCCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTa                    >  4:257882/1‑148 (MQ=255)                                                                                                                             cgccgcCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCAATTACCCGGCAACCCCCCAACCCCAAGGGAATTTAGAAAGCATAAACCATTGACTTTACTGTTGGAGCGGGCGTTCCCGCCAAAAAGGCCGCCCAAGCATa                    >  2:64330/1‑148 (MQ=255)                                                                                                                              gccgccCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                                               <  6:150419/104‑1 (MQ=38)                                                                                                                              gccgccCTTGAGGGGCGGCGTTTTCGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTAc                                                               >  5:150419/1‑104 (MQ=38)                                                                                                                                              ggcgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAAt  <  5:341728/149‑1 (MQ=255)                                                                                                                                               gcgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                         <  4:178886/93‑1 (MQ=255)                                                                                                                                               gcgTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTgg                                                         >  3:178886/1‑93 (MQ=255)                                                                                                                                                      |                                                                                                                                             AGTACATCGATCCGAAGTACATCGACATCTTCGTCGCGTGCAAGGAAAGCGAGCTGGAAGAGTTCGAGCACTCGATCTCCGACCTCGAGTA‑CAACTGGTACCTGCATACCGTGTAAACGAAAACGCCGCCCCTCAAGGGCGGCGTTTTTGTTTGAGACCGGTGTCGCCCTCTTCGCGGGCACGCCCGCTCCCGCAGGGTAATTTCGAAGGCTGAAGCCTGTGACTTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGCCCAGGCTTACAGGGGCTTCTCGAAAAT  >  NC_002947/5913993‑5914281

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