Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,884,963	C→A	8.3%	E139* (GAA→TAA) ‡	yedI ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,884,963	0	C	A	8.3%	226.1 / 8.2	72	E139* (GAA→TAA) ‡	yedI	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/35);  new base A (3/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCAACAAC  >  NC_002947/4884818‑4885111                                                                                                                                                  |                                                                                                                                                     ccAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaa                                                                                                                                                   <  4:328878/149‑1 (MQ=255) ccAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaa                                                                                                                                                   <  4:27739/149‑1 (MQ=255) ccAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaa                                                                                                                                                    >  3:328878/1‑148 (MQ=255)  cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  <  6:372764/149‑1 (MQ=255)  cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  <  2:157079/149‑1 (MQ=255)  cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  <  3:372908/149‑1 (MQ=255)  cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaa                                                                                                                                                    >  7:48129/1‑147 (MQ=255)  cAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaa                                                                                                                                                    >  3:352575/1‑147 (MQ=255)  cAATGGTCATCACCACCGCAATGCCGGAAAGCACCAAGCTCTCCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGAGGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  <  8:255683/149‑1 (MQ=255)  cAATGGTCATCACCACCGCAATCCCGGAAAGCACGATGATCTGATGGGTCAGCGGCTAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCaaaa                                                                                                                                                  <  3:69760/149‑1 (MQ=255)       gTCATCACCACCTCAATGCCGGAAAGCACGATGATCTGCGGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATCGCGATGATTTCTGCCGACAGGATGAAGTCGCTACGCACCGCCACCTTCATCTTGCCTTTTTACAAAGccac                                                                                                                                             >  3:40480/1‑149 (MQ=255)       gTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTGAAAAAGccaa                                                                                                                                             <  5:75784/149‑2 (MQ=255)          atcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAg                                                                                                                                          >  8:25542/1‑149 (MQ=255)          atcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAg                                                                                                                                          <  7:290628/149‑1 (MQ=255)          atcaCCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTGAAAAAGccaaga                                                                                                                                           <  6:105772/148‑4 (MQ=255)               caccGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCaa                                                                                                                                     >  3:2217/1‑149 (MQ=255)                    cAATGCCGGAAAGCACGATGATCTGCGGGGGCAGCGGCTAGGCAGCCACGATACCCAAGGTAATTGCGATTCTTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTg                                                                                                                                <  6:215286/149‑1 (MQ=255)                         ccGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATc                                                                                                                           <  1:158181/149‑1 (MQ=255)                          cGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCg                                                                                                                          >  6:181489/1‑149 (MQ=255)                          cGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCg                                                                                                                          >  3:48119/1‑149 (MQ=255)                                    gatgatCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGcct                                                                                                                >  4:42144/1‑149 (MQ=255)                                           tgctgGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCg                                                                                                                          >  3:119917/1‑132 (MQ=255)                                           tgctgGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCg                                                                                                                          <  4:119917/132‑1 (MQ=255)                                                ggTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCt                                                                                                    <  5:181489/149‑1 (MQ=255)                                                 gTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCCTGCGCGCCTc                                                                                                   <  4:177115/149‑1 (MQ=255)                                                      cggcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCt                                                                                              >  1:218115/1‑149 (MQ=255)                                                        gcgAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCtgt                                                                                            <  3:306135/149‑1 (MQ=255)                                                             tCAGCCACGATTCCCAAGGTAATCGCGATGATTTCTGCCGAGAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCt                                                                                        <  2:373825/148‑1 (MQ=255)                                                                gCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATc                                                                                    <  7:25542/149‑1 (MQ=255)                                                                 ccACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTATGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATcc                                                                                   >  5:255112/1‑149 (MQ=255)                                                                        aCCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTcc                                                                             >  5:384036/1‑148 (MQ=255)                                                                        aCCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTcc                                                                             <  6:384036/148‑1 (MQ=255)                                                                        aCCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCt                                                                            >  7:261591/1‑149 (MQ=255)                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                          >  4:6741/1‑149 (MQ=255)                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                          <  4:2217/149‑1 (MQ=255)                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                          >  1:232439/1‑149 (MQ=255)                                                                          ccAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                          <  2:232439/149‑1 (MQ=255)                                                                          ccAAGGTAATGGCGATGATTTCTGCAGACAGGATGAAGTCGGTACGCACCGCGCGGTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCATCGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                          <  4:40480/149‑1 (MQ=255)                                                                            aaGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGcc                                                                                                     >  1:206372/1‑120 (MQ=255)                                                                            aaGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGcc                                                                                                     <  2:206372/120‑1 (MQ=255)                                                                             aGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCtg                                                                       >  8:310185/1‑149 (MQ=255)                                                                                tAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGc                                                                    <  7:310185/149‑1 (MQ=255)                                                                                 aaTGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGc                                                                                               >  7:226137/1‑121 (MQ=255)                                                                                 aaTGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGc                                                                                               <  8:226137/121‑1 (MQ=255)                                                                                       gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATTTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACt                                                             >  3:244410/1‑149 (MQ=255)                                                                                       gatgatTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACt                                                             >  5:255491/1‑149 (MQ=255)                                                                                            tttCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTgc                                                        <  1:168494/149‑1 (MQ=255)                                                                                                 gCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAg                                                   >  1:53514/1‑149 (MQ=255)                                                                                                  ccGACAGGATGAAGTCGGTATGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGc                                                  <  6:255112/149‑1 (MQ=255)                                                                                                     aCAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTg                                                                                                                                >  4:257259/1‑68 (MQ=255)                                                                                                     aCAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTg                                                                                                                                <  3:257259/68‑1 (MQ=255)                                                                                                     aCAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTtc                                               <  4:352575/149‑1 (MQ=255)                                                                                                        ggATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                            <  3:6741/149‑1 (MQ=255)                                                                                                         gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                            >  8:10757/1‑148 (MQ=255)                                                                                                         gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCg                                            >  7:270092/1‑148 (MQ=255)                                                                                                         gATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGa                                           <  2:218115/149‑1 (MQ=255)                                                                                                         gATGAAGTCGGTAAGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGCCAACGGCCTCATTGCGCGCCACGTGCTTTGCCTCAACCTCTTCATTGCTGTGCAAGAACTTGTGCCCCACCTTCACg                                            >  6:368269/1‑148 (MQ=255)                                                                                                          aTGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGat                                          <  6:255491/149‑2 (MQ=255)                                                                                                              aGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAAccc                                      >  7:123527/1‑149 (MQ=255)                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                          >  5:97562/1‑102 (MQ=255)                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTg                                                                          <  6:97562/102‑1 (MQ=255)                                                                                                                         aCCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAg                           >  1:4012/1‑149 (MQ=255)                                                                                                                           cGCCCCCTTGATCTTGGCTTTTACAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAgg                          >  2:246820/1‑148 (MQ=255)                                                                                                                                      tCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCAcc               <  2:53514/148‑1 (MQ=255)                                                                                                                                       ggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGc                               <  7:124469/129‑1 (MQ=255)                                                                                                                                       ggtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGc                               >  8:124469/3‑131 (MQ=255)                                                                                                                                        gtGGCTTTTGAAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAg            >  5:258820/2‑149 (MQ=255)                                                                                                                                          ggCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCa          >  1:222314/1‑149 (MQ=255)                                                                                                                                           gCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCAt         <  3:42144/149‑1 (MQ=255)                                                                                                                                           gCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCAt         <  4:48119/149‑1 (MQ=255)                                                                                                                                             tttttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATc        >  7:43554/1‑148 (MQ=255)                                                                                                                                                ttCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATcaac     <  8:48129/148‑1 (MQ=255)                                                                                                                                                  cAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCaacaac  <  1:246820/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCAATGGTCATCACCACCGCAATGCCGGAAAGCACGATGATCTGCTGGGTCAGCGGCGAGTCAGCCACGATACCCAAGGTAATGGCGATGATTTCTGCCGACAGGATGAAGTCGGTACGCACCGCCCCCTTGATCTTGGCTTTTTCAAAAGCCACCAGGTCAACACTGGCATCGGCAACGGCCTCCTTGCGCGCCTCGTGCTGTGCCTCATCCTCTTCCTTGCTGTGCAAGAACTTGTGCGCCAGCTTCTCGAAACCCTCGAAGCACAGGAACGCGCCACCCAGCATCAACAAC  >  NC_002947/4884818‑4885111

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