Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,819,672	A→C	6.3%	E306A (GAA→GCA) ‡	rpoS →	RNA polymerase sigma S factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,819,672	0	A	C	6.3%	198.1 / 4.5	64	E306A (GAA→GCA) ‡	rpoS	RNA polymerase sigma S factor
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base C (2/2);  total (33/31)
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.28e-01

ACCTGTCGCAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCG  >  NC_002947/1819534‑1819815                                                                                                                                           |                                                                                                                                                 acCTGTCGCAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCgg                                                                                                                                        <  5:221242/149‑1 (MQ=255) acCTGTCGCAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCgg                                                                                                                                        <  7:208562/149‑1 (MQ=255)      tCGCAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCaaa                                                                                                                                    <  2:271212/148‑1 (MQ=255)        gCAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATc                                                                                                                                  >  4:313489/1‑148 (MQ=255)        gCAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATc                                                                                                                                  >  5:299164/1‑148 (MQ=255)         cAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATc                                                                                                                                  >  3:173141/1‑147 (MQ=255)         cAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCa                                                                                                                                <  3:313489/149‑1 (MQ=255)           aaGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAgg                                                                                                                              <  6:197067/149‑1 (MQ=255)           aaGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAgg                                                                                                                              <  4:173141/149‑1 (MQ=255)           aaGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAgg                                                                                                                              <  3:55223/149‑1 (MQ=255)                  ggtCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAgg                                                                                                                       <  8:197662/147‑1 (MQ=255)                   aTCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGa                                                                                                                                                  <  5:67584/121‑1 (MQ=255)                   aTCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGa                                                                                                                                                  >  6:67584/1‑121 (MQ=255)                      aGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTc                                                                                                                   <  2:199014/149‑1 (MQ=255)                          gTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAgc                                                                                                               >  2:13603/1‑149 (MQ=255)                          gTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAgc                                                                                                               >  2:262759/1‑149 (MQ=255)                              ggTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCtt                                                                                                           >  4:191328/1‑149 (MQ=255)                                tGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGACTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTgc                                                                                                         >  3:22278/1‑149 (MQ=255)                                   gTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTgcgc                                                                                                       >  6:267701/1‑148 (MQ=255)                                           aCAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCgag                                                                                                    <  5:329332/143‑1 (MQ=255)                                           aCAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCgag                                                                                                    >  6:329332/1‑143 (MQ=255)                                           aCAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCtt                                                                                              >  7:24912/1‑149 (MQ=255)                                           aCAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCtt                                                                                              <  3:191328/149‑1 (MQ=255)                                                cagcGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAg                                                                                                                               >  1:222625/1‑111 (MQ=255)                                                cagcGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAg                                                                                                                               <  2:222625/111‑1 (MQ=255)                                                cagcGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGaa                                                                                           >  5:238682/1‑147 (MQ=255)                                                  gcGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAg                                                                                        >  4:70209/1‑148 (MQ=255)                                                   cGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAg                                                                                        <  4:60773/147‑1 (MQ=255)                                                   cGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGaa                                                                                      >  3:60773/1‑149 (MQ=255)                                                       gtgtcgTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGt                                                                                  <  1:8183/144‑1 (MQ=255)                                                       aaGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGt                                                                                  <  7:214823/149‑1 (MQ=255)                                                            gTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTTCACGGGtaagg                                                                                                                                        >  1:273207/1‑86 (MQ=255)                                                            gTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTTCACGGGtaagg                                                                                                                                        <  2:273207/90‑5 (MQ=255)                                                                  cgccgcTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGc                                                                         >  7:153468/1‑147 (MQ=255)                                                                   gccgcTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAg                                                                      >  1:72550/1‑149 (MQ=255)                                                                       cTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCAt                                                                  <  2:225327/149‑1 (MQ=255)                                                                            gCCTTCGCGGCCATGAAAGCAGTACCCTTGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCATGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTc                                                             <  8:24912/149‑1 (MQ=255)                                                                            gCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTc                                                             <  1:262759/149‑1 (MQ=255)                                                                              cTTCGCGGCCATGAAAGCAGTACCCTGGAGGAGGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCa                                                           <  2:127345/149‑1 (MQ=255)                                                                               ttCGCGGTCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGACCGGGTGCGGCAAATCCAGGTGGAAGGGCGCAAGCGCTTGCGCGAGATCCTGGAAAAGAATTGTCGGTCCAGCGAGTCATGGTTACAg                                                          >  6:253687/1‑149 (MQ=255)                                                                                  gcgGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGtag                                                       <  1:13603/149‑1 (MQ=255)                                                                                  gcgGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGta                                                        <  3:70209/148‑1 (MQ=255)                                                                                      ccATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGtagctag                                                   <  4:22278/149‑1 (MQ=255)                                                                                      ccATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGtagctag                                                   >  2:267852/1‑149 (MQ=255)                                                                                          gAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCt                                               >  7:182378/1‑149 (MQ=255)                                                                                          gAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGAGAGATCCTTGAAAATAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCt                                               >  5:325527/1‑149 (MQ=255)                                                                                             tgCAGTACCGTTGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACTGGTGCGGCATATCCAGGTGTAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAATGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTa                                            <  5:141092/148‑1 (MQ=255)                                                                                                  tACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACaa                                        >  3:174225/1‑148 (MQ=255)                                                                                                    cccTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACaaaa                                      >  1:290325/1‑148 (MQ=255)                                                                                                       tGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGt                                                               <  3:91329/120‑1 (MQ=255)                                                                                                       tGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGt                                                               >  4:91329/1‑120 (MQ=255)                                                                                                       tGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGtagcta                                                    >  4:113622/1‑131 (MQ=255)                                                                                                       tGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGtagcta                                                    <  3:113622/131‑1 (MQ=255)                                                                                                              tGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGt                           >  1:226606/1‑149 (MQ=255)                                                                                                                     cTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCg                   >  7:44607/1‑149 (MQ=255)                                                                                                                        gagaTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCggg                 <  8:182378/148‑1 (MQ=255)                                                                                                                        gagaTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCgg                  <  6:299164/147‑1 (MQ=255)                                                                                                                           aTCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTATCAAAAAACGCCCCGGTGATC‑TCGGGGCGt             >  3:46410/1‑149 (MQ=255)                                                                                                                            tCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGtt            >  3:89189/1‑149 (MQ=255)                                                                                                                              ggCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGtttt          >  1:170342/1‑149 (MQ=255)                                                                                                                               gCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCC‑GAGATCACCGGGGCGttttt         <  6:29623/149‑1 (MQ=255)                                                                                                                                ccgcaCCCGTTCACGGGTGCGGCAAATCCAGGTGTAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGttt           >  1:49660/5‑146 (MQ=255)                                                                                                                                ccgcaCCCGTTCACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGttt           <  2:49660/142‑1 (MQ=255)                                                                                                                                    aCCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTtgtg      <  4:174225/147‑1 (MQ=255)                                                                                                                                      ccGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTgcg  >  1:327581/1‑149 (MQ=255)                                                                                                                                      ccGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGGTTTTGTGTgcg  >  4:168496/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                 ACCTGTCGCAAAGCATCGATCAGTGGTTGGGTGAGTTGACCGACAAGCAGCGTGAAGTGGTTGTACGCCGCTTTGGCCTTCGCGGCCATGAAAGCAGTACCCTGGAGGATGTTGGGCTGGAGATCGGCCTTACCCGTGAACGGGTGCGGCAAATCCAGGTGGAAGGGCTCAAGCGCTTGCGCGAGATCCTTGAAAAGAATGGTCTGTCCAGCGAGTCATTGTTCCAGTAGCTAGCGCTGTAACAAAAAACGCCCCGGTGATC‑TCGGGGCGTTTTTGTGTGCG  >  NC_002947/1819534‑1819815

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