Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,115,851	A→C	8.8%	S16A (TCG→GCG) ‡	wcaM ←	colanic acid biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,115,851	0	A	C	8.8%	234.4 / 13.9	80	S16A (TCG→GCG) ‡	wcaM	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base C (2/5);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-02

ACCTGGCTGCCGTCCTGTAAGACAAATCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTATTCCAC  >  NC_000913/2115586‑2116118                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                            aCCTGGCTGCCGTCCTGTAAGACAAATCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGt                                                                                                                                                                                                                                                                    >  2:420286/1‑275 (MQ=255)                          aTCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGGTGGGAGGGtttt                                                                                                                                                                                                                                           >  2:107215/1‑275 (MQ=255)                                   ttccGTTGCCGCGCAAGCTTCCAAGAATGTGTTACGTTTTACCGGGAGGGATAAAGATTCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGc                                                                                                                                                                                                                                  <  2:142562/273‑1 (MQ=255)                                    gCCGTGGCCGCTGAAGCGTTCAAGAGTGTGTAACGTAGTACCGGGAGGGATAAAGATGCCGGGGTCGAGATATGCACGAACCAATCCGTCAGGGACGTGGACTGGTTGACCTTAGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTGGTGAATGTTAACGCCGTGTCGGGCGGGAAGCGGGCGGGCGAAAGGGGTATGGGGGAAGGCTAGCGCCGGGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGc                                                                                                                                                                                                                                  <  4:339209/274‑1 (MQ=255)                                     ccGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTTTTTTCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGGGGGTTGTAGTCGTTACTGTTAACGATTGGTCGGGCGGGAAGCCCGGGGGCGAAAGGGGAATGGAGGAAGGAAAGCGCCGAACTTGCCGTCAGGAACGTGCGACGGGCAGGTTTTTTAAATGGCAt                                                                                                                                                                                                                                >  4:235470/1‑275 (MQ=255)                                                                    aCGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtat                                                                                                                                                                                                                  <  2:16267/258‑1 (MQ=255)                                                                               ggagggaTAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                                                                                                                                                        >  2:121076/1‑273 (MQ=255)                                                                                gagtgaTAAAGAGGCCTGTGTTGAGATTGTCACAAACCAATCCGGCAGGCAGGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                                                                                                      <  2:229693/274‑1 (MQ=255)                                                                                 gtggaTAAATGTGCCGGTGTGGATATTGTCAAAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGa                                                                                                                                                                                    <  1:69259/273‑1 (MQ=255)                                                                                 agggaTAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCACTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGa                                                                                                                                                                                    >  2:98057/1‑275 (MQ=255)                                                                                   ggaAAAAGATTCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                                                                                                                                                  <  2:403816/275‑1 (MQ=255)                                                                                             gCCGGTGGGGAGATGGTCACAAACCAATCCGGGAGGAAAGACGACTGTTTGCCCTTCGCTGAAGGATTTGTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAAGGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTTCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                                                                                                         <  4:420281/274‑1 (MQ=255)                                                                                                                      aaTCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTCGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGa                                                                                                                                                                                                                                                                             <  8:389622/149‑1 (MQ=255)                                                                                                                                gcaCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCTCGGCGCTTGCCTTc                                                                                                                                                                                                                                                                   <  6:303051/149‑1 (MQ=255)                                                                                                                                 caCGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGACGGCAAGCTCGGCGCTTGCCTTcc                                                                                                                                                                                                                                                                  <  8:288212/149‑2 (MQ=255)                                                                                                                                     acgacTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTga                                                                                                                                <  2:185540/275‑1 (MQ=255)                                                                                                                                          tGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGc                                                                                                                                                                                                                                                         >  6:225584/1‑149 (MQ=255)                                                                                                                                            tttGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAGAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCANGAACGGGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAACCTGCAGCAAGCTGTCGAGTTCTCGATTTNAACCCTGCTGGTTAAAATCGCGTTCNACTTTTTCGTGCGCTCCttttcagcccgg                                                                                                                         >  3:93066/1‑266 (MQ=255)                                                                                                                                               gCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCgg                                                                                                                                                                                                                                                    <  7:150987/149‑1 (MQ=255)                                                                                                                                                cccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                                                                                                                                                   <  6:455648/149‑1 (MQ=255)                                                                                                                                                cccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCggg                                                                                                                                                                                                                                                   >  7:267550/1‑149 (MQ=255)                                                                                                                                                 ccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGga                                                                                                                                                                                                                                                  <  8:242456/149‑1 (MQ=255)                                                                                                                                                 ccTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGga                                                                                                                                                                                                                                                  <  7:27532/149‑1 (MQ=255)                                                                                                                                                    tCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTTCAGCAAGCTGGCGAGTtct                                                                                                                                                                                       >  2:81484/1‑205 (MQ=255)                                                                                                                                                       cTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttt                                                                                                                                                                                                                                            <  6:257080/149‑1 (MQ=255)                                                                                                                                                        tGAAGGCTTGTTTAAATGAGGCGATCCAGGCGTGTGGGTTGTAGTCGTTAATGGTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATtcgt                                                                                                              <  1:48952/274‑1 (MQ=255)                                                                                                                                                         gAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGttttt                                                                                                                                                                                                                                          <  8:420934/149‑1 (MQ=255)                                                                                                                                                          aaGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGtttttt                                                                                                                                                                                                                                         >  5:130733/1‑149 (MQ=255)                                                                                                                                                             gCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTaaa                                                                                                                                                                                                                                      >  6:462239/1‑149 (MQ=255)                                                                                                                                                              cTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAAt                                                                                                                                                                                                                                     <  5:462239/149‑1 (MQ=255)                                                                                                                                                                           ggCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGCAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTtctc                                                                                                                                                                                                                        <  7:50040/149‑1 (MQ=255)                                                                                                                                                                            gCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                                                                                                                                                        >  1:24840/1‑180 (MQ=255)                                                                                                                                                                            gCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAAGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAAcg                                                                                         >  2:86145/1‑275 (MQ=255)                                                                                                                                                                                tCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTCTAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAAcgccgc                                                                                     >  2:29681/1‑275 (MQ=255)                                                                                                                                                                                tCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCCTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTAGGGAGAGTTTTTTAAATTGCATGCGTTCTCCTCta                                                                                                                                                                                                                   >  8:288751/1‑149 (MQ=255)                                                                                                                                                                                tCCAGGTGTGTGGGTTGTAGTCGTTAATGTTAACGCTGTGTNGGGCGGNAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGa                                                                                                                                                                                    <  3:24840/180‑1 (MQ=255)                                                                                                                                                                                        tgtgGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCt                                                                                                                                                                                                           <  8:267550/149‑1 (MQ=255)                                                                                                                                                                                             gTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTctcctc                                                                                                                                                                                                                     >  3:3004/1‑134 (MQ=255)                                                                                                                                                                                              ttGTGGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                                                                                                                                                 <  7:207816/137‑1 (MQ=255)                                                                                                                                                                                              ttGTGGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                                                                                                                                                 >  8:207816/1‑137 (MQ=255)                                                                                                                                                                                               tGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCaa                                                                                                                                                                                                    <  5:271628/149‑1 (MQ=255)                                                                                                                                                                                                 tAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCtata                                                                                                                                                                                                                 <  1:3004/134‑1 (MQ=255)                                                                                                                                                                                                   gTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTg                                                                                                                                                                                                >  6:162461/1‑149 (MQ=255)                                                                                                                                                                                                   gTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTg                                                                                                                                                                                                >  6:61880/1‑149 (MQ=255)                                                                                                                                                                                                      gTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCa                                                                                                                                                                             >  2:328926/1‑165 (MQ=255)                                                                                                                                                                                                       ttAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGa                                                                                                                                                                                            >  5:45219/1‑149 (MQ=255)                                                                                                                                                                                                       ttAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGa                                                                                                                                                                                            >  7:405831/1‑149 (MQ=255)                                                                                                                                                                                                          aTGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                                                                                                         <  4:328922/165‑1 (MQ=255)                                                                                                                                                                                                                 cgcttgggggggggggAAGCGCGGGGTGGGATGGGGGATGGAGGAAGGCAAGCGCGTGGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGAATGCGTTCTGCGCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATTCCCTGTTGGTTAAACTCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTtct                                                     <  4:272989/265‑1 (MQ=255)                                                                                                                                                                                                                 cGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGAtt                                                                                                                                                                                  >  5:95078/1‑149 (MQ=255)                                                                                                                                                                                                                    tttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                                                                                                                                               <  5:61880/149‑1 (MQ=255)                                                                                                                                                                                                                    tttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGAt                                                                                                                                                                               <  8:189754/149‑1 (MQ=255)                                                                                                                                                                                                                      tGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAACCGCCGCCAGTCGTTGCGCCAGTTCGCGAGCATCGTTCTCAGGCa                                               >  2:297862/1‑275 (MQ=255)                                                                                                                                                                                                                      tGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCa                                               >  2:319321/1‑275 (MQ=255)                                                                                                                                                                                                                             cgggAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgctg                                                                                                                                                                      <  8:405831/149‑1 (MQ=255)                                                                                                                                                                                                                                  aGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCgg                                   >  2:151124/1‑275 (MQ=255)                                                                                                                                                                                                                                  aGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCAAGCCg                                    >  2:479263/1‑274 (MQ=255)                                                                                                                                                                                                                                         gggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                                                                                                                                          <  5:225584/149‑1 (MQ=255)                                                                                                                                                                                                                                          ggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtgt                                                                                                                                                         <  6:95078/149‑1 (MQ=255)                                                                                                                                                                                                                                          ggCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTCCGCCAGTGCGCGAGCATCGCTTTCAGGCACCAGCCAGCCGGGTTTGtg                            >  2:436794/1‑273 (MQ=255)                                                                                                                                                                                                                                             gAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTCAAAATCGTGTTCa                                                                                                                                                      >  8:301687/1‑149 (MQ=255)                                                                                                                                                                                                                                                   ggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTccac                   <  1:93066/274‑1 (MQ=255)                                                                                                                                                                                                                                                     tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                                                                                                                                              <  8:4903/149‑1 (MQ=255)                                                                                                                                                                                                                                                     tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                                                                                                                                              <  5:162461/149‑1 (MQ=255)                                                                                                                                                                                                                                                     tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTc                                                                                                                                              >  7:222668/1‑149 (MQ=255)                                                                                                                                                                                                                                                        ggaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcg                                                                                                                                           >  5:221454/1‑149 (MQ=255)                                                                                                                                                                                                                                                         gaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGGGTCGGGAGAGTTTTTTAATTGGCATGCATTCTCCCCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTc             <  2:335947/274‑1 (MQ=255)                                                                                                                                                                                                                                                         gaggaAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGGTCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgc                                                                                                                                          >  8:36864/1‑149 (MQ=255)                                                                                                                                                                                                                                                         gagcaAGGAAAGCGCGGCGCGGGCCGTCAGGAGCGTGCGTTGGTAGAGTTTTTTAGATGGCGTCCGTTCTCCTATATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTGTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTc             <  4:189341/274‑1 (MQ=255)                                                                                                                                                                                                                                                           tgacggCAAGCTCGGCGCTTGCCTTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACtttt                                                                                                                                                >  7:317627/5‑141 (MQ=255)                                                                                                                                                                                                                                                           tgacggCAAGCTCGGCGCTTGCCTTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACtttt                                                                                                                                                <  8:317627/137‑1 (MQ=255)                                                                                                                                                                                                                                                           ggagtgCAAGGGCCGTGGTTGGCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTGTAAAGCCGGCAGGAAGCTGGCGAGTTCTCGATTGATCACCGGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCg           <  4:319317/269‑1 (MQ=255)                                                                                                                                                                                                                                                             aaGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTccacca                 >  2:189343/1‑266 (MQ=255)                                                                                                                                                                                                                                                              aGGCGAGCGCCGACCTTGCCGTGCGGAGAGTGGATAGGGAGAGATGGTTAAATGGCTGTCTCTGTCCTCTACAAAACGTTCGGTAAGTCTGCAACGTTTCGTGTGATCACCTTCTGCTTAAAGTTGTGTTTGACTTTTACGCGCTAGAGGTTGACAACAGGTGCCGAGTCGTCGGTGTCCAGTTGTGTAAAGGCCGCCAGTCGTTTAGACAGTGCGCGGGCATCGTTCTCAGGCACCAGCCAGCCGGATTGGTCAGCCTCCACCAGTTCCGGTa        <  4:436789/274‑1 (MQ=255)                                                                                                                                                                                                                                                               ggCAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGTAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGt                                                                                                                                    >  5:389623/1‑149 (MQ=255)                                                                                                                                                                                                                                                               ggCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCACTTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTAtt      >  2:456296/1‑275 (MQ=255)                                                                                                                                                                                                                                                               ggCAAGCGCCGAGCTTGCAGCCAGGAACGTGCGTCGGGAGAGGTTTTTTAATGGCATGCGTTCTTCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGCTCACCTGCTGGTTAAACCCGGGTTCAACTCTTTCGCGCGCGCGTTTGACCACCGGAGCCACTTCGTCGGTTTCCAGTTCGCTAAGCCCCCCCAGCCCTTGCGACCGTGCGCGAGCATCGGTCTCAGGCCCCATCCCACCGGATTTTTCCCCCCCCACCCGCTCTGGCAtt      >  4:248320/1‑275 (MQ=255)                                                                                                                                                                                                                                                                 cAAGCTCGGCGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGttt                                                                                                                                  <  6:389623/149‑1 (MQ=255)                                                                                                                                                                                                                                                                   aGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAACTCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGGTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTATTCCAc  >  2:160028/1‑275 (MQ=255)                                                                                                                                                                                                                                                                      gcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTgacga                                                                                                                             >  6:83775/1‑149 (MQ=255)                                                                                                                                                                                                                                                                        cGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGAcc                                                                                                                           >  5:120907/1‑149 (MQ=255)                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                            ACCTGGCTGCCGTCCTGTAAGACAAATCGCCCTCTGCCGTTGCCGCGCAAGCTTCCAAGAATGTGTAACGTTTTACCGGGAGGGATAAAGATGCCGGTGTTGATATTGTCACAAACCAATCCGGCAGGCACGACGACTGTTTGCCCTTCGCTGAAGGCTTGTTTAAATGAGGCGATCCAGTCGTGTGGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCGTTGCGCCAGTGCGCGAGCATCGTTCTCAGGCACCAGCCAGCCGGATTTGTCAGCCTCCACCAGTTCCGGTATTCCAC  >  NC_000913/2115586‑2116118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑