Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,940,404	T→C	5.1%	K411R (AAA→AGA)	mepM ←	murein DD‑endopeptidase, space‑maker hydrolase, septation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,940,404	0	T	C	5.1%	243.9 / 2.5	78	K411R (AAA→AGA)	mepM	murein DD‑endopeptidase, space‑maker hydrolase, septation protein
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base C (2/2);  total (41/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.34e-01

CACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGT  >  NC_000913/1940259‑1940551                                                                                                                                                  |                                                                                                                                                    cACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTg                                                                                                                                                  >  2:470129/1‑149 (MQ=255) cACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGGTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGGCAGCCCTTCGGTACGCGGCAGTTTTg                                                                                                                                                  <  6:14141/149‑1 (MQ=255)    cAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCg                                                                                                                                               >  3:420112/1‑149 (MQ=255)     aaaTAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGt                                                                                                                                              >  6:99403/1‑149 (MQ=255)     aaaTAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGt                                                                                                                                              >  3:152064/1‑149 (MQ=255)                   ggTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACgg                                                                                                                                <  1:470129/149‑1 (MQ=255)                        tgactgCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCctgc                                                                                                                           >  1:88709/1‑149 (MQ=255)                          actgCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCctgctg                                                                                                                         >  3:67555/1‑149 (MQ=255)                              cGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCGTTCGGTACGCGGCTGGTTTGCCGTCAGCGGGTTTACGGCCTGCTGGttt                                                                                                                     >  4:247893/1‑149 (MQ=255)                              cGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGttt                                                                                                                     <  2:182887/149‑1 (MQ=255)                                tACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGc                                                                                                                               >  5:162205/1‑137 (MQ=255)                                tACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGc                                                                                                                               <  6:162205/137‑1 (MQ=255)                                  ccGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATcc                                                                                                                 <  6:119621/149‑1 (MQ=255)                                   cGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCa                                                                                                                <  4:420112/149‑1 (MQ=255)                                       tGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACt                                                                                                            >  6:114983/1‑149 (MQ=255)                                            aTGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCa                                                                                                         >  3:311831/1‑147 (MQ=255)                                             tGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAg                                                                                                      >  7:181932/1‑149 (MQ=255)                                               gATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTg                                                                                                    >  5:294262/1‑149 (MQ=255)                                                aTGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGc                                                                                                   <  5:99403/149‑1 (MQ=255)                                                aTGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGc                                                                                                   <  2:88709/149‑1 (MQ=255)                                                  gTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTCTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATTCTTCATAGTg                                                                                                    >  2:181611/1‑146 (MQ=255)                                                  gTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTCTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATTCTTCATAGTg                                                                                                    >  2:175157/1‑146 (MQ=255)                                                   ttaattaaTCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGa                                                                                                <  4:311831/149‑1 (MQ=255)                                                     aattaatCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATg                                                                                              >  3:26983/1‑149 (MQ=255)                                                     aattaatCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACAATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACg                                                                                                                                 >  1:113821/1‑114 (MQ=255)                                                     aattaatCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACAATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACg                                                                                                                                 <  2:113821/114‑1 (MQ=255)                                                      attaatCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATc                                                                                                                  <  5:281137/128‑1 (MQ=255)                                                      attaatCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATc                                                                                                                  >  6:281137/1‑128 (MQ=255)                                                         aatCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATa                                                                                                       <  7:99402/136‑1 (MQ=255)                                                         aatCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATa                                                                                                       >  8:99402/1‑136 (MQ=255)                                                          atCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGcc                                                                                         >  1:67296/1‑149 (MQ=255)                                                           tCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGc                                                                                                                                                 <  2:235773/92‑1 (MQ=255)                                                           tCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGccc                                                                                        >  4:194124/1‑149 (MQ=255)                                                            cAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCg                                                                                       >  3:396125/1‑149 (MQ=255)                                                              aaCCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGTTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTTTc                                                                                                                  <  5:225223/120‑1 (MQ=255)                                                              aaCCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATc                                                                                                                  >  6:225223/1‑120 (MQ=255)                                                              aaCCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGt                                                                                     >  1:209485/1‑149 (MQ=255)                                                                 cGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGGCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGa                                                                                  <  4:67555/149‑1 (MQ=255)                                                                   tAGCTGCGGGACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGTCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAAc                                                                                <  5:60488/149‑1 (MQ=255)                                                                      cTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTc                                                                             <  4:152064/149‑1 (MQ=255)                                                                        gCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCg                                                                           >  4:98152/1‑149 (MQ=255)                                                                          ggCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGt                                                                         <  1:384725/149‑1 (MQ=255)                                                                           gCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTa                                                                        <  4:396125/149‑1 (MQ=255)                                                                            cacaATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTAt                                                                       >  7:134257/1‑149 (MQ=255)                                                                              caATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGGACGCGGCAGTTTTGCCGGCAGCGGGTTTCCGGCCGGCTGGTTTATCCATACTTCAGAGTGCAGATGCGGCCCGGTGGCACGTCCGGTATTa                                                                     >  4:123288/1‑149 (MQ=255)                                                                               aaTCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTAc                                                                    >  4:264257/1‑149 (MQ=255)                                                                                       tGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTGCGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGaaa                                                              <  8:455617/147‑1 (MQ=255)                                                                                       tGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGaaa                                                              <  2:209485/147‑1 (MQ=255)                                                                                       tGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAg                                                             <  2:67296/148‑1 (MQ=255)                                                                                         gCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAgcgc                                                          >  7:280399/1‑149 (MQ=255)                                                                                                 cAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGt                                                  >  2:22186/1‑149 (MQ=255)                                                                                                 cAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGt                                                  >  7:436367/1‑149 (MQ=255)                                                                                                 cAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGt                                                  >  8:168553/1‑149 (MQ=255)                                                                                                  agtaaTTCGCGACGCTCGGAGCCGGTCCGCCCTTCGGTCCGCGGCAGTTTTGCCGTCAGCGGGTGTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGACCGTCCGGTAGTACCGGAAAGCGCGATACGGTc                                                 <  3:123288/146‑1 (MQ=255)                                                                                                  aGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTc                                                 <  3:264257/149‑1 (MQ=255)                                                                                                    gAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAgcg                                                           >  5:294173/1‑137 (MQ=255)                                                                                                    gAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAgcg                                                           <  6:294173/137‑1 (MQ=255)                                                                                                     aaTTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGcc                                              <  5:114983/149‑1 (MQ=255)                                                                                                      aTTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGAGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATTCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCa                                             >  2:432559/1‑149 (MQ=255)                                                                                                          gcgACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGt                                          <  3:98152/148‑1 (MQ=255)                                                                                                            gacgaTCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGTCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTc                                       <  1:22186/149‑1 (MQ=255)                                                                                                              cgaTCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCAc                                     <  7:168553/149‑1 (MQ=255)                                                                                                                    gAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTg                               >  4:417827/1‑149 (MQ=255)                                                                                                                     aGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGt                              <  3:417827/149‑1 (MQ=255)                                                                                                                        cGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTccc                           <  8:436367/149‑1 (MQ=255)                                                                                                                        cGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTccc                           >  7:285656/1‑149 (MQ=255)                                                                                                                        cGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTccc                           >  4:333074/1‑149 (MQ=255)                                                                                                                         ggTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCg                          <  4:26983/149‑1 (MQ=255)                                                                                                                           tCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGt                        >  7:241063/1‑149 (MQ=255)                                                                                                                           tCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGt                        >  6:107914/1‑149 (MQ=255)                                                                                                                               cccTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCa                    <  8:241063/149‑1 (MQ=255)                                                                                                                                ccTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCAc                   >  3:93376/1‑149 (MQ=255)                                                                                                                                    cGGTACGCGGCAGTCTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATTCTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCa                 <  1:181611/147‑1 (MQ=255)                                                                                                                                    cGGTACGCGGCAGTCTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATTCTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCa                 <  1:175157/147‑1 (MQ=255)                                                                                                                                      gTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAAt             >  4:46333/1‑149 (MQ=255)                                                                                                                                      gTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAAt             >  4:167877/1‑149 (MQ=255)                                                                                                                                       tACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATc            <  8:280399/149‑1 (MQ=255)                                                                                                                                          gcgGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTg         <  3:167877/149‑1 (MQ=255)                                                                                                                                                 ttttGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGt  >  1:72494/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGGCACAATCTCTTTGGCCTGCGCCAGGAATTCGCGACGATCGGAGCCGGTCAGCCCTTCGGTACGCGGCAGTTTTGCCGTCAGCGGGTTTACGGCCTGCTGGTTTATCCATACTTCATAGTGCAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGT  >  NC_000913/1940259‑1940551

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑