Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	313,966	A→T	7.7%	V192E (GTG→GAG)	yrfG ←	purine nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	313,966	0	A	T	7.7%	245.8 / 8.9	91	V192E (GTG→GAG)	yrfG	purine nucleotidase
Reads supporting (aligned to +/- strand):  ref base A (44/39);  new base T (5/2);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

GACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATAG  >  NC_002947/313822‑314112                                                                                                                                                 |                                                                                                                                                     gactgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑Gc                                                                                                                                                 >  8:2417/1‑149 (MQ=255)  actgaTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑Gc                                                                                                                                                 <  2:47638/148‑1 (MQ=255)      aTCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGc                                                                                                                                             >  6:194073/1‑148 (MQ=255)       tCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGc                                                                                                                                             >  8:221482/1‑147 (MQ=255)       tCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCaa                                                                                                                                           <  5:78933/149‑1 (MQ=255)       tCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCaa                                                                                                                                           <  3:185364/149‑1 (MQ=255)       tCGCTGAACTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGAAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCCCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCaa                                                                                                                                           <  7:9715/149‑1 (MQ=255)               gTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGTACC‑GCCAGCAAGTGCCCca                                                                                                                                   >  6:215348/1‑149 (MQ=255)                      cAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCaaa                                                                                                                            >  5:43343/1‑149 (MQ=255)                       aGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAAc                                                                                                                           >  6:76037/1‑149 (MQ=255)                         gCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAAccg                                                                                                                         >  6:77069/1‑149 (MQ=255)                           ctctTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAAcagcc                                                                                                                       >  8:108530/1‑149 (MQ=255)                                cgcgGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCa                                                                                                                  >  6:272015/1‑149 (MQ=255)                                cgcgGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCACACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGt                                                                                                                     <  5:135105/146‑1 (MQ=255)                                 gcgGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCAc                                                                                                                 <  6:151450/149‑1 (MQ=255)                                  cgGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACt                                                                                                                >  4:231462/1‑149 (MQ=255)                                    ggCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTgc                                                                                                              <  3:231462/149‑1 (MQ=255)                                    ggCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTgc                                                                                                              <  8:278270/149‑1 (MQ=255)                                      cgcgcCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTgcgc                                                                                                            <  7:221482/149‑1 (MQ=255)                                       gcgcCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCa                                                                                                           <  5:77069/149‑1 (MQ=255)                                         gcCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTACTCCTCGACCGCCGCACACTCTGCAGGATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCCCC‑GCCAGCACGTGCCCCACACCAAACCGCCGTGCACTGCGCAg                                                                                                          >  5:34339/1‑148 (MQ=255)                                           ccGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGTCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAAt                                                                                                       <  6:172347/149‑1 (MQ=255)                                           ccGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCaa                                                                                                                                           <  2:280822/113‑1 (MQ=255)                                           ccGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCaa                                                                                                                                           >  1:280822/1‑113 (MQ=255)                                                  tGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGc                                                                                                >  3:214933/1‑149 (MQ=255)                                                    gggCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgt                                                                                              <  5:194073/149‑1 (MQ=255)                                                      gCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTg                                                                                                                                        <  2:33789/105‑1 (MQ=255)                                                      gCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTg                                                                                                                                        >  1:33789/1‑105 (MQ=255)                                                           aCAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCa                                                                                                                  >  8:179400/1‑122 (MQ=255)                                                            cAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGa                                                                                     <  6:43343/149‑1 (MQ=255)                                                             aGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCGACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGAt                                                                                    <  7:108530/149‑1 (MQ=255)                                                             aGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGAt                                                                                    <  5:76037/149‑1 (MQ=255)                                                              ggCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATg                                                                                   <  5:272015/149‑1 (MQ=255)                                                                ccTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCAGGCGCACC‑GCCAGCAAGTGCCCCACCCCAAACCGCCGTGCACTGCGAAGAATCGCCAGCGTGTC‑ATCGATGaa                                                                                 >  1:223956/1‑149 (MQ=255)                                                                    ccAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAgg                                                                             >  5:2868/1‑149 (MQ=255)                                                                    ccAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAg                                                                              <  5:224665/148‑1 (MQ=255)                                                                     cAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGc                                                                            >  2:865/1‑149 (MQ=255)                                                                     cAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGc                                                                            >  7:131156/1‑149 (MQ=255)                                                                              ccGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCg                                                                   >  6:282298/1‑149 (MQ=255)                                                                              ccGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCAGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAACCCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTcag                                                                   >  1:62357/1‑147 (MQ=255)                                                                                 gTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGtt                                                                >  3:145065/1‑149 (MQ=255)                                                                                   aaTCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACa                                                                               <  6:18084/132‑1 (MQ=255)                                                                                   aaTCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACa                                                                               >  5:18084/1‑132 (MQ=255)                                                                                   aaTCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTc                                                               <  8:131156/148‑1 (MQ=255)                                                                                     tCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgt                                                                                              >  3:32740/1‑116 (MQ=255)                                                                                     tCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCgt                                                                                              <  4:32740/116‑1 (MQ=255)                                                                                     tCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCg                                                              >  5:57482/1‑147 (MQ=255)                                                                                      ccTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTc                                                                     <  4:236864/139‑1 (MQ=255)                                                                                      ccTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTc                                                                     >  3:236864/1‑139 (MQ=255)                                                                                       cTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCTCC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGc                                                          >  7:46584/1‑149 (MQ=255)                                                                                       cTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGc                                                          >  3:204044/1‑149 (MQ=255)                                                                                              gccgcAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCTCC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATc                                                   <  8:46584/149‑1 (MQ=255)                                                                                                cgcAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCg                                                                   <  5:214781/131‑1 (MQ=255)                                                                                                cgcAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCg                                                                   >  6:214781/1‑131 (MQ=255)                                                                                                   aaaCTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCt                                              >  2:288274/1‑149 (MQ=255)                                                                                                    aaCTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTg                                             <  4:214933/149‑1 (MQ=255)                                                                                                       tctTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGGCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTg                                             >  1:71479/1‑147 (MQ=255)                                                                                                         ttCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAggg                                        >  6:304375/1‑149 (MQ=255)                                                                                                         ttCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAggg                                        >  5:244861/1‑149 (MQ=255)                                                                                                           cAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTc                                                               >  3:129574/1‑124 (MQ=255)                                                                                                           cAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTc                                                               <  4:129574/124‑1 (MQ=255)                                                                                                           cAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCa                                      >  1:293681/1‑149 (MQ=255)                                                                                                             gtctCCCGCGGCCCTGCTTTGCTAGCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCCCACCAAACCGCCGTGCACTGCGCAGACTCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATc                                    <  6:34339/146‑1 (MQ=255)                                                                                                             gTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATc                                    <  4:145065/149‑1 (MQ=255)                                                                                                                    gcgGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAAc                             >  7:144074/1‑149 (MQ=255)                                                                                                                       gCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGc                          <  1:288274/149‑1 (MQ=255)                                                                                                                         ccTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCgg                        <  4:204044/149‑1 (MQ=255)                                                                                                                         ccTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCgg                        <  8:144074/149‑1 (MQ=255)                                                                                                                          cTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCggg                       >  5:214989/1‑149 (MQ=255)                                                                                                                           tGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGc                      >  4:187592/1‑149 (MQ=255)                                                                                                                            gCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCt                     >  4:267811/1‑149 (MQ=255)                                                                                                                                 gCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCt                >  2:158884/1‑149 (MQ=255)                                                                                                                                    cacttgCTGGCGGTGC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcg             >  2:148141/6‑149 (MQ=255)                                                                                                                                     tCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcgc            >  6:52640/1‑149 (MQ=255)                                                                                                                                     tCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTc              <  7:2417/147‑1 (MQ=255)                                                                                                                                      ccGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcgc            <  2:293681/148‑1 (MQ=255)                                                                                                                                       ttgCTGGCGGTGC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGCAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcgcg           >  1:203650/3‑148 (MQ=255)                                                                                                                                       ttgCTGGCGGTGC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCaaa                               >  1:248745/3‑128 (MQ=255)                                                                                                                                       ttgCTGGCGGTGC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCaaa                               <  2:248745/126‑1 (MQ=255)                                                                                                                                       ttgCTGGCGGTGC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcgcg           >  2:235165/3‑148 (MQ=255)                                                                                                                                        ggCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTcgcg           >  8:69252/1‑147 (MQ=255)                                                                                                                                        ggCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTa         <  7:69252/149‑1 (MQ=255)                                                                                                                                             gcgcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGAt                                                                                    <  2:149712/69‑1 (MQ=255)                                                                                                                                             gcgcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGAt                                                                                    <  2:144553/69‑1 (MQ=255)                                                                                                                                             gcgcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGAt                                                                                    >  1:149712/1‑69 (MQ=255)                                                                                                                                             gcgcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGAt                                                                                    >  1:144553/1‑69 (MQ=255)                                                                                                                                             gcgcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCTCGTAGCCAt    <  2:62357/149‑1 (MQ=255)                                                                                                                                              cgcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATa   <  6:2868/149‑1 (MQ=255)                                                                                                                                              cgcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATa   <  5:304375/149‑1 (MQ=255)                                                                                                                                               gcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTCAAGCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCa     >  1:275638/1‑147 (MQ=255)                                                                                                                                               gcACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATAg  >  1:209281/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GACTGATCGCTGAAGTCTGTACAGGCCTCTTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACC‑GCCAGCAAGTGCCCCACACCAAACCGCCGTGCACTGCGCAGAATCGCCAGCGTGTC‑ATCGATGAACAGGCTGCGTTCCGGTTCGAAGCCGATATCCGCCTGCAGGGCATCCCAAAACTGCGGGCTTTCCTTCGCGTAGCCATAG  >  NC_002947/313822‑314112

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