Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,884,962	T→G	6.0%	E139A (GAA→GCA) ‡	yedI ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,884,962	0	T	G	6.0%	210.3 / 3.4	67	E139A (GAA→GCA) ‡	yedI	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base G (2/2);  total (30/37)
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.83e-01

CACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCAACAA  >  NC_002947/4884814‑4885110                                                                                                                                                     |                                                                                                                                                     cacacCAATGGTCATCACCACCGCAATGCCGTAAAGCACTATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATTGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGTCttttt                                                                                                                                                      <  2:346030/149‑1 (MQ=255)  acacCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCttttt                                                                                                                                                      >  6:229798/1‑148 (MQ=255)      cgataagcatcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                 <  1:106918/142‑1 (MQ=255)      cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                 <  5:147406/149‑1 (MQ=255)      cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                 <  3:221992/149‑1 (MQ=255)             catcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGTGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTGAAAAAGccaaga                                                                                                                                          <  2:56778/149‑4 (MQ=255)                caccacCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTACAAAAGCCACCAGGt                                                                                                                                       >  3:261053/1‑149 (MQ=255)                 accaccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGCCTTATTCAAAAGCCACCACGt                                                                                                                                       >  4:230819/1‑148 (MQ=255)                   caccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCaa                                                                                                                                    <  5:229798/149‑1 (MQ=255)                               ggAAAGCACGATGATCTGCTGTGTCATCGGCGAGTCAGCCACGATACCCTAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCgg                                                                                                                        <  8:197190/149‑1 (MQ=255)                                gAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAaca                                                                                                                                  >  4:248091/1‑138 (MQ=255)                                gAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAaca                                                                                                                                  <  3:248091/138‑1 (MQ=255)                                  aaGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCaa                                                                                                                     >  7:32397/1‑149 (MQ=255)                                     cACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACgg                                                                                                                  >  5:13743/1‑149 (MQ=255)                                      aCGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGc                                                                                                                 >  4:211098/1‑149 (MQ=255)                                      aCGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGc                                                                                                                 >  7:93760/1‑149 (MQ=255)                                          tgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                 <  6:44967/113‑1 (MQ=255)                                          tgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                 >  5:44967/1‑113 (MQ=255)                                                gctgGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTgcgc                                                                                                       >  1:176569/1‑149 (MQ=255)                                                 ctgGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCtt                                                                                                           <  1:17525/144‑1 (MQ=255)                                                    ggTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTgcgcgc                                                                                                     <  3:230819/147‑1 (MQ=255)                                                      tCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCg                                                                                                 <  8:235761/149‑1 (MQ=255)                                                          cggcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCt                                                                                             <  6:13743/149‑1 (MQ=255)                                                            gagAGTCAGCCACGATTCCCAAGGTAATGGAGATGTTTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCTTCGGCAACGGCCTCCTTGCGCGCCTCGTGCtgt                                                                                           <  4:261053/147‑1 (MQ=255)                                                              gAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGc                                                                                         <  7:156750/149‑1 (MQ=255)                                                                    gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   >  6:119472/1‑149 (MQ=255)                                                                    gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   <  5:119472/149‑1 (MQ=255)                                                                    gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   <  8:93760/149‑1 (MQ=255)                                                                    gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   <  8:32397/149‑1 (MQ=255)                                                                       aCGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCctc                                                                                >  8:88076/1‑149 (MQ=255)                                                                       aCGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCctc                                                                                >  4:109424/1‑149 (MQ=255)                                                                          aTACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTc                                                                             <  3:211098/149‑1 (MQ=255)                                                                             cccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCtt                                                                          <  2:176569/149‑1 (MQ=255)                                                                              ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCAt                                                                                    >  4:92727/1‑138 (MQ=255)                                                                              ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCAt                                                                                    <  3:92727/138‑1 (MQ=255)                                                                              ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   >  6:347844/1‑139 (MQ=255)                                                                              ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   <  5:347844/139‑1 (MQ=255)                                                                              ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                         >  5:320066/1‑149 (MQ=255)                                                                                       tGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTgcgc                                                                                                       <  6:263705/110‑1 (MQ=255)                                                                                        ggCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAg                                                                >  5:46153/1‑148 (MQ=255)                                                                                         gCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   >  2:38849/1‑128 (MQ=255)                                                                                         gCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                   <  1:38849/128‑1 (MQ=255)                                                                                           gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACt                                                            >  2:101107/1‑149 (MQ=255)                                                                                           gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACt                                                            >  4:304323/1‑149 (MQ=255)                                                                                                tttCTGCCGACAGGATGAAGTCGGTACGCACCGCCGCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgc                                                       <  5:198433/149‑1 (MQ=255)                                                                                                tttCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgc                                                       <  4:328888/149‑1 (MQ=255)                                                                                                 ttCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCATGCACTTGTgcg                                                      >  6:189814/1‑149 (MQ=255)                                                                                                      ccGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCa                                                   <  6:320066/147‑1 (MQ=255)                                                                                                          cAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtct                                             <  6:175621/149‑1 (MQ=255)                                                                                                           aGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtctc                                            <  3:109424/149‑1 (MQ=255)                                                                                                             gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGc                                                 <  6:340677/142‑1 (MQ=255)                                                                                                             gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGc                                                 >  5:340677/1‑142 (MQ=255)                                                                                                               tGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                           >  8:3930/1‑146 (MQ=255)                                                                                                               tGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                           <  7:3930/146‑1 (MQ=255)                                                                                                                  aGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtct                                             >  6:108850/1‑141 (MQ=255)                                                                                                                  aGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtct                                             <  5:108850/141‑1 (MQ=255)                                                                                                                       gTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGACACCCTc                                   >  4:130132/1‑146 (MQ=255)                                                                                                                       gTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAACCCCTc                                   <  3:130132/146‑1 (MQ=255)                                                                                                                        tACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGACAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAg                               >  6:333149/1‑149 (MQ=255)                                                                                                                           gCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGcac                            <  1:101107/149‑1 (MQ=255)                                                                                                                              ccGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAgg                         <  7:88076/149‑1 (MQ=255)                                                                                                                                        gATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCa                <  3:304323/148‑1 (MQ=255)                                                                                                                                          tggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGc                                                 <  8:142081/110‑1 (MQ=255)                                                                                                                                          tggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGc                                                 >  7:142081/4‑113 (MQ=255)                                                                                                                                            ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                                             >  8:121879/1‑99 (MQ=255)                                                                                                                                            ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAAc                                                             <  7:121879/99‑1 (MQ=255)                                                                                                                                            ttGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAg           >  2:8557/1‑149 (MQ=255)                                                                                                                                            gtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAg           >  2:225267/2‑149 (MQ=255)                                                                                                                                                cTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCAt        >  7:303100/1‑148 (MQ=255)                                                                                                                                                 tttttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATc       >  5:292260/1‑148 (MQ=255)                                                                                                                                                 tttttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATc       >  3:332689/1‑148 (MQ=255)                                                                                                                                                 tttttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATc       >  1:340587/1‑148 (MQ=255)                                                                                                                                                     tCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCACCTTGCGCCCCTCGTGCTGTGCCTCATCCTCTACCTTGCTCTGCAAGAACATCTGCGCCAGCATCTCGAAACCCACGAAGCCCAGGAACGCCCCACCCACCATcaacaa  >  7:271059/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CACACCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCAACAA  >  NC_002947/4884814‑4885110

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