Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,137,960	A→T	7.1%	E20V (GAA→GTA)	yegH →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,137,960	0	A	T	7.1%	193.3 / 5.8	56	E20V (GAA→GTA)	yegH	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (38/14);  new base T (2/2);  total (40/16)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

TGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAGACAGCAATAA  >  NC_000913/2137717‑2138248                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                 tgtCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCTGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCACg                                                                                                                                                                                                                                                      >  1:44751/1‑288 (MQ=255)                                   tattattaCCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCgcgc                                                                                                                                                                                                                    >  1:399713/1‑288 (MQ=255)                                          aCCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCt                                                                                                                                                                                                                                                                                         <  2:91806/212‑1 (MQ=255)                                          aCCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCt                                                                                                                                                                                                                                                                                         >  1:91806/1‑212 (MQ=255)                                                                    aGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGTACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAAt                                                                                                                                                                                   >  1:399688/1‑288 (MQ=255)                                                                         attaCCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATgcg                                                                                                                                                                                >  2:587658/1‑286 (MQ=255)                                                                                    ataataTGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCCATGTTCATGCGCCTGTTCctgccg                                                                                                                                                                   >  1:236925/1‑286 (MQ=255)                                                                                              aTAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAAtct                                                                                                                                                         >  1:24949/1‑288 (MQ=255)                                                                                                  aTAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCctg                                                                                                                                                     >  1:108707/1‑288 (MQ=255)                                                                                                      tgctgcAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAag                                                                                                                                                                                                                                  <  1:10071/207‑1 (MQ=255)                                                                                                      tgctgcAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAag                                                                                                                                                                                                                                  >  2:10071/1‑207 (MQ=255)                                                                                                        ctgcAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTc                                                                                                                                              >  1:446212/1‑289 (MQ=255)                                                                                                                     tttttGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCaa                                                                                                                                  >  1:661446/1‑288 (MQ=255)                                                                                                                       tttGCCCCGATTTTTCCATTATCGAATTGGTACATTTGTCATACAAGCAAGAGGTTTTGTACTTACTATGGAATGAATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCg                                                                                                                               <  1:374955/289‑1 (MQ=255)                                                                                                                          gCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGTACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAAcc                                                                                                                                >  1:152660/1‑285 (MQ=255)                                                                                                                                     tCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGGGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCCATCTCCTGGCTTGTCACCCTTACTCCACCGCGTTTCAGGTTcc                                                                                                                  >  1:642670/1‑288 (MQ=255)                                                                                                                                          gATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGt                                                                                                            >  1:656013/1‑289 (MQ=255)                                                                                                                                            tCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGTACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGtt                                                                                                           <  2:152660/288‑1 (MQ=255)                                                                                                                                            tCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGtt                                                                                                           <  2:661446/288‑1 (MQ=255)                                                                                                                                                      aCATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAgc                                                                                                 >  2:183480/1‑288 (MQ=255)                                                                                                                                                      aCATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGAGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCAGCTCGTTTACCTTTAgc                                                                                                 >  2:461347/1‑288 (MQ=255)                                                                                                                                                                             ttttGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGtt                                                                          <  1:183480/288‑1 (MQ=255)                                                                                                                                                                               ttGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCt                                                                                                                                         >  1:529310/1‑223 (MQ=255)                                                                                                                                                                               ttGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCt                                                                                                                                         <  2:529310/223‑1 (MQ=255)                                                                                                                                                                               ttGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTACCGCCCGCGACTTACACCTGCTGTTtgg                                                                       >  1:356370/1‑289 (MQ=255)                                                                                                                                                                                tGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTtggt                                                                      >  1:493780/1‑289 (MQ=255)                                                                                                                                                                                tGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTtggt                                                                      >  1:1329/1‑289 (MQ=255)                                                                                                                                                                                tGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGCTACTGCTGGACTCAATCTCCTGGCCGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCTACTTAATCACGCTGTTtggg                                                                      >  1:642964/1‑288 (MQ=255)                                                                                                                                                                                   acttactATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGTACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTtggtgg                                                                    <  2:399688/288‑1 (MQ=255)                                                                                                                                                                                      tactATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCCCGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGtt                                                                  >  2:403523/1‑287 (MQ=255)                                                                                                                                                                                       actATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATg                                                                               <  1:207018/273‑1 (MQ=255)                                                                                                                                                                                       actATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATg                                                                               >  2:207018/1‑273 (MQ=255)                                                                                                                                                                                       actATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTttct                                                               >  2:109568/1‑289 (MQ=255)                                                                                                                                                                                       actATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTttct                                                               <  2:446212/289‑1 (MQ=255)                                                                                                                                                                                       actATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTttct                                                               >  2:395740/1‑289 (MQ=255)                                                                                                                                                                                       actATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTACTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTttct                                                               >  1:214478/1‑289 (MQ=255)                                                                                                                                                                                            ggaatggaTTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCCCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTTGTttctt                                                              >  1:404485/1‑285 (MQ=255)                                                                                                                                                                                               atggaTTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCtgttg                                                       >  2:190279/1‑289 (MQ=255)                                                                                                                                                                                                  gaTTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTgttgtt                                                     <  2:642670/288‑1 (MQ=255)                                                                                                                                                                                                   aTTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCAcc                                                                                                                                           <  1:252223/201‑1 (MQ=255)                                                                                                                                                                                                   aTTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCAcc                                                                                                                                           >  2:252223/1‑201 (MQ=255)                                                                                                                                                                                                     tGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCaaa                                                 >  2:277077/1‑289 (MQ=255)                                                                                                                                                                                                       ccGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCCGTTTGGTGGTTTCTTCCTGTTGTTCAAAg                                                >  1:351424/1‑288 (MQ=255)                                                                                                                                                                                                               tCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATgg                                       >  1:302595/1‑289 (MQ=255)                                                                                                                                                                                                                 tATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAg                                     <  2:656013/289‑1 (MQ=255)                                                                                                                                                                                                                       ggCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTgaac                               >  2:483107/1‑289 (MQ=255)                                                                                                                                                                                                                         ccGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATTGAGCTgaacga                             >  2:342397/1‑289 (MQ=255)                                                                                                                                                                                                                            ggTTAATCACGCTGATTGTGATCGAACTGGTCCTCTGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAATCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCCATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTTACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAAGCATGCTGTTTGGTGCTTTCTTCCTGTTGTTCAACTCTACGATGGAGCTGaacgaccg                          >  1:605107/1‑289 (MQ=255)                                                                                                                                                                                                                             gTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACgg                         >  1:627782/1‑289 (MQ=255)                                                                                                                                                                                                                                      gCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAggg                >  2:116664/1‑289 (MQ=255)                                                                                                                                                                                                                                          aTTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAg            >  1:674572/1‑289 (MQ=255)                                                                                                                                                                                                                                                aTCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAGACAGCa      >  1:345156/1‑289 (MQ=255)                                                                                                                                                                                                                                                aTCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCCCCTGGCTGGTCACACTGACTCAACCGCTGTTCAGCTTCAGCTCGTTTACCTTTAGCGCCCGCGACTTAACCATGCTGTTTGGTGGTTTCTCCCTGTTGTTCAAAGCCTCGATGGAGCGGAACGAACGGCCGGAAGGGATAGACCGCa      >  1:312672/1‑289 (MQ=255)                                                                                                                                                                                                                                                 tCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAGACAGCaa     <  2:493780/289‑1 (MQ=255)                                                                                                                                                                                                                                                 tCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAGACAGCaa     <  2:627782/289‑1 (MQ=255)                                                                                                                                                                                                                                                 tCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAGACAGCaa     <  2:108707/289‑1 (MQ=255)                                                                                                                                                                                                                                                    aaCTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATg                                                                               <  2:113588/212‑1 (MQ=255)                                                                                                                                                                                                                                                    aaCTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATg                                                                               >  1:113588/1‑212 (MQ=255)                                                                                                                                                                                                                                                    aaCTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCCGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAGACAGCaataa  >  2:562128/1‑289 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                 TGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTACCGGGCTGCTGCTGGCAATGTTAATGCGCCTGTTACTGCTGGCGTCAATCTCCTGGCTGGTCACCCTGACTCAACCGCTGTTCAGCTTCCGCTCGTTTACCTTTAGCGCCCGCGACTTAATCATGCTGTTTGGTGGTTTCTTCCTGTTGTTCAAAGCCACGATGGAGCTGAACGAACGGCTGGAAGGGAAAGACAGCAATAA  >  NC_000913/2137717‑2138248

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