Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,984,835	T→C	7.1%	E95E (GAA→GAG)	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,984,835	0	T	C	7.1%	292.3 / 8.7	98	E95E (GAA→GAG)	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (49/42);  new base C (2/5);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGC  >  NC_002947/5984688‑5984978                                                                                                                                                    |                                                                                                                                                     gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCTGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCtt                                                                                                                                                     >  2:339172/1‑149 (MQ=255)   aGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg                                                                                                                                                   >  7:94940/1‑149 (MQ=255)   aGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCg                                                                                                                                                   >  1:109564/1‑149 (MQ=255)    ggTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGt                                                                                                                                                  >  8:263503/1‑149 (MQ=255)     gTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt                                                                                                                                                 >  1:152754/1‑149 (MQ=255)     gTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt                                                                                                                                                 >  7:4218/1‑149 (MQ=255)     gTGACCGATGCCTCGGGCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGtt                                                                                                                                                 >  7:27213/1‑149 (MQ=255)         ccGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTcc                                                                                                                                             <  7:128521/149‑1 (MQ=255)                aTCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGt                                                                                                                                      >  3:123113/1‑149 (MQ=255)                aTCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTAACGAAGTCTTCGTTGTCCCAGGCGt                                                                                                                                      >  3:373132/1‑149 (MQ=255)                    gTCACCACGAACACTTCAAAACCTTCCGCCCGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAt                                                                                                                                  <  8:344454/149‑1 (MQ=255)                       accacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCTGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTg                                                                                                                               <  1:339172/149‑1 (MQ=255)                        ccacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTgg                                                                                                                              <  5:141962/149‑1 (MQ=255)                        ccacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTgg                                                                                                                              >  8:161544/1‑149 (MQ=255)                        ccacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTgg                                                                                                                              <  5:64055/149‑1 (MQ=255)                         cacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACTCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGc                                                                                                                             >  7:81594/1‑149 (MQ=255)                         cacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGc                                                                                                                             >  7:272183/1‑149 (MQ=255)                           cGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCg                                                                                                                           <  8:4218/149‑1 (MQ=255)                           cGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCg                                                                                                                           <  6:75096/149‑1 (MQ=255)                                  ttCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACggg                                                                                                                    >  8:38458/1‑149 (MQ=255)                                  ttCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTAGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACggg                                                                                                                    >  7:9352/1‑149 (MQ=255)                                   tCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt                                                                                                                                     >  8:4238/1‑131 (MQ=255)                                   tCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt                                                                                                                                     <  7:4238/131‑1 (MQ=255)                                      aagccTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAAt                                                                                                                >  7:206223/4‑149 (MQ=255)                                      aaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACgg                                                                                                                     >  6:380333/1‑144 (MQ=255)                                      aaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACgg                                                                                                                     <  5:380333/144‑1 (MQ=255)                                      aaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAAt                                                                                                                >  2:355845/1‑149 (MQ=255)                                          cTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTAc                                                                                                            <  7:263503/149‑1 (MQ=255)                                          cTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTAc                                                                                                            <  7:126159/149‑1 (MQ=255)                                             ccGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACgg                                                                                                          <  7:336195/148‑1 (MQ=255)                                                 cAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCAt                                                                                                     <  1:355845/149‑1 (MQ=255)                                                  aGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATc                                                                                                    >  6:251158/1‑149 (MQ=255)                                                    cgcCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCg                                                                                                  <  4:123113/149‑1 (MQ=255)                                                     gcCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCgg                                                                                                 <  8:206223/149‑1 (MQ=255)                                                        gACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGaa                                                                                              >  1:311527/1‑149 (MQ=255)                                                         aCAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAAc                                                                                             <  3:224220/149‑1 (MQ=255)                                                            gggTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGc                                                                                                                             <  3:328336/114‑1 (MQ=255)                                                            gggTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGc                                                                                                                             >  4:328336/1‑114 (MQ=255)                                                                 gggAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCttt                                                                                     <  8:272183/149‑1 (MQ=255)                                                                  ggAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTc                                                                                    <  8:94940/149‑1 (MQ=255)                                                                   gAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTCGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                   <  8:9352/149‑1 (MQ=255)                                                                    aaGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAg                                                                                  >  1:39401/1‑149 (MQ=255)                                                                       gCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTc                                                                               >  5:125406/1‑149 (MQ=255)                                                                       gCCACGCACACATCGGTCACCACACCGGCGAAGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTc                                                                               >  1:42708/1‑149 (MQ=255)                                                                        ccACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCg                                                                              <  5:251158/149‑1 (MQ=255)                                                                         cACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCgg                                                                             >  1:258316/1‑149 (MQ=255)                                                                          aCGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCATGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCggg                                                                            <  5:31321/149‑1 (MQ=255)                                                                             cacacaTCGGTCACCACTCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CAc                                                                      <  8:81594/149‑1 (MQ=255)                                                                               cacaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTCACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCa                                                                    <  4:373132/149‑1 (MQ=255)                                                                                acaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                   <  6:333138/136‑1 (MQ=255)                                                                                acaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATCTACGGCGCATCCGGGAACATCTCTTTCa                                                                                   >  5:333138/1‑136 (MQ=255)                                                                                    cGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCgg‑g                                                             <  2:389177/149‑1 (MQ=255)                                                                                    cGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCgg‑g                                                             >  6:305306/1‑149 (MQ=255)                                                                                      gTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTc                                                                                    <  6:364263/129‑1 (MQ=255)                                                                                      gTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTc                                                                                    >  5:364263/1‑129 (MQ=255)                                                                                         accacACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCAtt                                                        >  3:381649/1‑149 (MQ=255)                                                                                           cacaCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCg                                                      >  3:52780/1‑149 (MQ=255)                                                                                            acacCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCgg                                                     <  2:42708/149‑1 (MQ=255)                                                                                                cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCt                                                 >  4:35237/1‑149 (MQ=255)                                                                                                cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCt                                                 >  3:205370/1‑149 (MQ=255)                                                                                                cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCt                                                 >  5:140999/1‑149 (MQ=255)                                                                                                   cGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGc                                              <  8:27213/149‑1 (MQ=255)                                                                                                   cGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGc                                              <  7:161544/149‑1 (MQ=255)                                                                                                    gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCt                                             >  6:260696/1‑149 (MQ=255)                                                                                                        atCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCg                                           >  5:182692/1‑147 (MQ=255)                                                                                                         tCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAg                                        >  4:266994/1‑149 (MQ=255)                                                                                                         tCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAg                                        >  4:261754/1‑149 (MQ=255)                                                                                                           aGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCt                                      >  6:167555/1‑149 (MQ=255)                                                                                                            gTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtg                                     >  4:107213/1‑149 (MQ=255)                                                                                                            gTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtg                                     >  2:368668/1‑149 (MQ=255)                                                                                                            gTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtg                                     >  3:10022/1‑149 (MQ=255)                                                                                                                cttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtggtgg                                 <  5:305306/149‑1 (MQ=255)                                                                                                                 ttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTggaggt                                >  7:208326/1‑149 (MQ=255)                                                                                                                    cGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                              <  4:205370/148‑1 (MQ=255)                                                                                                                    cGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                              <  2:2631/148‑1 (MQ=255)                                                                                                                     ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAgg                            <  4:381649/149‑1 (MQ=255)                                                                                                                     ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAg                             >  1:334892/1‑148 (MQ=255)                                                                                                                      gCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                           <  1:368668/149‑1 (MQ=255)                                                                                                                       cccGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                          <  5:260696/149‑1 (MQ=255)                                                                                                                         cGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATgg                        >  7:315786/1‑149 (MQ=255)                                                                                                                                 tcaccgCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtgg                                    <  2:79639/124‑1 (MQ=11)                                                                                                                                 tcaccgCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCtgg                                    >  1:79639/6‑129 (MQ=11)                                                                                                                                  gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                 <  7:38458/147‑1 (MQ=255)                                                                                                                                  gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                 <  2:39401/147‑1 (MQ=255)                                                                                                                                   accgCCTAGACGAAATCCTCCTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGcc              <  8:261008/146‑1 (MQ=11)                                                                                                                                   aTTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCACGCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                      <  1:81514/142‑1 (MQ=255)                                                                                                                                   aTTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                      >  2:81514/1‑141 (MQ=255)                                                                                                                                       ccTTGACGAAATCCTCGTTGTCCCCGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAg          <  6:214536/149‑1 (MQ=14)                                                                                                                                        cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAg          >  3:270087/1‑148 (MQ=255)                                                                                                                                           gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGc                                                   <  3:53795/104‑1 (MQ=14)                                                                                                                                           gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGc                                                   >  4:53795/1‑104 (MQ=14)                                                                                                                                            aCGAACGCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACtt     <  2:152754/149‑1 (MQ=255)                                                                                                                                             cGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCGTGc                                              >  4:330258/1‑107 (MQ=255)                                                                                                                                             cGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg    <  2:52431/149‑1 (MQ=255)                                                                                                                                              gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTGGTGTTAGAGGA‑CCGGAGGCAATAGGGCCTTGATCGAAGCTGGTGGTGAGGATGTTCGGCAGGCCAAAGAAc       >  1:324330/1‑149 (MQ=18)                                                                                                                                              gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                       <  6:100416/129‑1 (MQ=255)                                                                                                                                              gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                       >  5:100416/1‑129 (MQ=255)                                                                                                                                               aaGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc  <  6:125406/149‑1 (MQ=255)                                                                                                                                                aGTCTTCGTTGTCACAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc  >  5:129721/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGG‑‑‑CACCA‑GCGG‑GCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGC  >  NC_002947/5984688‑5984978

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