Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	834,278	G→A	W403* (TGG→TGA)	dinG →	ATP‑dependent DNA helicase DinG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	834,278	0	G	A	98.5%	223.6 / NA	67	W403* (TGG→TGA)	dinG	ATP‑dependent DNA helicase DinG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (37/29);  minor base C (1/0);  total (38/29)
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.93e-01

TCAATCTCTACATGCCTGCCGGGCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAGGAGATGAGTGGTCTGAAAGAGAAAGCGGGCGACCGTT  >  NC_000913/834034‑834522                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     tcAATCTCTACATGCCTGCCGGGCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACggctg                                                                                                                                                                                                                                                >  1:194117/1‑251 (MQ=255)          aCATGCCTGCCGGGCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTg                                                                                                                                                                                                                                       >  2:145059/1‑251 (MQ=255)            aTGCCTGCCGGGCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCCTGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCNGgc                                                                                                                                                                                                                                     >  2:405474/1‑251 (MQ=255)                   ccGGGCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCt                                                                                                                                                                                                                              <  1:405474/251‑1 (MQ=255)                   ccGGGCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCt                                                                                                                                                                                                                              >  2:113898/1‑251 (MQ=255)                     gggCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTc                                                                                                                                                                                                                            <  1:393203/251‑1 (MQ=255)                         aggaggCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGt                                                                                                                                                                                                                        >  1:457723/1‑251 (MQ=255)                         aggaggCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGt                                                                                                                                                                                                                        >  1:287211/1‑251 (MQ=255)                          ggaggCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTg                                                                                                                                                                                                                       >  1:313/1‑251 (MQ=255)                           gaggCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGCCGCAAAGCAAACTCTGACGGCTGTCTTCGCTGgcgacatactccggtgc                                                                                                                                                                                                                      >  2:292139/1‑237 (MQ=255)                                agagCCCCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTCTTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCgg                                                                                                                                                                                                                 <  1:127725/251‑1 (MQ=255)                                agagCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCgg                                                                                                                                                                                                                 >  2:77944/1‑251 (MQ=255)                                      ccGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGAcc                                                                                                                                                                                                            >  1:138766/1‑250 (MQ=255)                                       cGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCa                                                                                                                                                                                                           >  1:367007/1‑250 (MQ=255)                                                  gggCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGAcgc                                                                                                                                                                                               >  1:212112/1‑251 (MQ=255)                                                                 tGAAGTGATGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTATTAACGATTTAAGGGAGAAAACTGGCAGCCATGACATTGTACGTCTGCATTGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGATCGAGGCGCAAAGAAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGTACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGa                                                                                                                                                                                  >  1:340896/1‑249 (MQ=255)                                                                            agaTCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTacac                                                                                                                                                                     >  2:120185/1‑251 (MQ=255)                                                                                  gCCAGAGGCAGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTATGACGGCTGGCTTCGCTGGCGCAATCTTCCGGGGCACCGGAGACCACATGGGCGCCGCGGGAAGAGCGCGAAGGGGAGCTACACCTCTgg                                                                                                                                                               >  2:260461/1‑251 (MQ=255)                                                                                    cAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGtt                                                                                                                                                             <  1:29352/251‑1 (MQ=255)                                                                                       cGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTcc                                                                                                                                                          >  2:155896/1‑250 (MQ=255)                                                                                              cAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTg                                                                                                                                                   <  2:313/251‑1 (MQ=255)                                                                                               aaaaCTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTgg                                                                                                                                                  >  2:138325/1‑251 (MQ=255)                                                                                                 aacgcACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGACAGCCATGACATTGTCCGTCTGCATCGGTTGATGTTGCAGATGTACCGCGCGTTGTTGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGTAAGAGCGCGAATGCCAGCTACACCTCTGGTTTCACTGCGTgtga                                                                                                                                                <  2:115462/247‑4 (MQ=255)                                                                                                     cACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTTACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATa                                                                                                                                             <  1:284510/250‑1 (MQ=255)                                                                                                     cACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATAc                                                                                                                                            <  1:280933/251‑1 (MQ=255)                                                                                                         gagaTGCTGCGTGGCCTGGGTGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACgtgt                                                                                                                                        <  2:287211/251‑1 (MQ=255)                                                                                                         gagaTGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCGGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACgtgt                                                                                                                                        <  2:91046/251‑1 (MQ=255)                                                                                                           gaTNCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCa                                                                                                                                      <  2:408813/251‑1 (MQ=255)                                                                                                            aTGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAg                                                                                                                                     <  2:263293/251‑1 (MQ=255)                                                                                                                 gcGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGAt                                                                                                                                 <  1:113898/250‑1 (MQ=255)                                                                                                                     ggCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCt                                                                                                                            >  2:202523/1‑251 (MQ=255)                                                                                                                          ggCGGTGTTTTTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGCTGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAACCTTTGACGGCTGGCTTCGCTGGCGCAATCTGCCGGTGCTCCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTTGGAATACGTGTCAGCGATCAGCTGGaaa                                                                                                                       <  2:232817/251‑1 (MQ=255)                                                                                                                          ggCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTAGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGaa                                                                                                                        <  1:357178/250‑1 (MQ=255)                                                                                                                            cGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAgg                                                                                                                     <  1:155692/251‑1 (MQ=255)                                                                                                                             ggAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGGATACGTGTCAGCGATAAGCTGGAAAGgc                                                                                                                    >  2:321241/1‑251 (MQ=255)                                                                                                                               aGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGgctg                                                                                                                  >  2:20520/1‑251 (MQ=255)                                                                                                                                 ttATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTACGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCGGCTACACCTCTGGTTTCACTGc                                                                                                                                                      >  1:180861/1‑213 (MQ=255)                                                                                                                                 ttATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTGGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCCCTGc                                                                                                                                                      <  2:180861/213‑1 (MQ=255)                                                                                                                                        ttAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgcgc                                                                                                         >  1:149670/1‑251 (MQ=255)                                                                                                                                         tAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACggctgg                                                                                                                                                                                                                                               >  1:325041/1‑116 (MQ=255)                                                                                                                                         tAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACcgctgg                                                                                                                                                                                                                                               <  2:325041/116‑1 (MQ=255)                                                                                                                                                 tAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGCCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGCAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAggatcccg                                                                                                >  2:350095/1‑244 (MQ=255)                                                                                                                                                  aaGTGATAAAACCGGCAGCCATGACATTTTACGTCTTCATCGGTTGATTTTGCCGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTTGCGCAATCTTCCGGTGCACCGGTGACCGAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTTAGCGATCAGCGGGAAAGTCTGCTGTGGCGCAGTATGCCg                                                                                                <  2:252809/250‑1 (MQ=255)                                                                                                                                                  aaGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCg                                                                                                <  1:20520/250‑1 (MQ=255)                                                                                                                                                   aGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGc                                                                                               >  1:207810/1‑250 (MQ=255)                                                                                                                                                       agaAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCACATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACatt                                                                                          >  1:137850/1‑251 (MQ=255)                                                                                                                                                       agaAAACCGGCAGCCATGACATTGTACGGCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACatt                                                                                          <  2:254662/251‑1 (MQ=255)                                                                                                                                                         aaaaCCGGAAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTATCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACattat                                                                                        >  2:23219/1‑251 (MQ=255)                                                                                                                                                          aaaCCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGATGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACattatt                                                                                       >  1:51967/1‑251 (MQ=255)                                                                                                                                                              cGGCAGCCATGACATTGGCCGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTTGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCa                                                                                   <  1:155896/251‑1 (MQ=255)                                                                                                                                                                      aTGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCa                                                                           <  2:367007/251‑1 (MQ=255)                                                                                                                                                                                aCGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGtt                                                                 <  2:250890/251‑1 (MQ=255)                                                                                                                                                                                        aTCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAAc                                                         >  2:8442/1‑251 (MQ=255)                                                                                                                                                                                        aTCGGATGATTTTGCAGATGAACCCCGCGTTGGGGATGTTCGACACGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTCACCAAATGGGCGCCGCGGGCCGAGCGCGAAGCGCAGCTACACCCCTGGCTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAGGGTCGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGCTCgatgcac                                                         >  2:91405/1‑245 (MQ=255)                                                                                                                                                                                         tCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCGGGCGCAATCTTCCGGTGCACCGGTGAACAAGTGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCGGGAAAGGCTGGTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTGCGCTGAAc                                                         >  1:207479/1‑250 (MQ=255)                                                                                                                                                                                          cGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAAc                                                         >  2:435542/1‑249 (MQ=255)                                                                                                                                                                                             ttGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGttt                                                    >  1:423658/1‑251 (MQ=255)                                                                                                                                                                                             ttGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGttt                                                    >  2:198714/1‑251 (MQ=255)                                                                                                                                                                                                     gCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGt                                             <  2:352371/250‑1 (MQ=255)                                                                                                                                                                                                         aTGAACCGCGCGTTGGGGAGGTTCGAGGGGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCa                                        <  2:297574/251‑1 (MQ=255)                                                                                                                                                                                                           gAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAgg                                      <  2:162590/251‑1 (MQ=255)                                                                                                                                                                                                              ccGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTATGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAgg                                      >  2:397615/1‑248 (MQ=255)                                                                                                                                                                                                                                   ggCGCAAAGCAAACTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAGGAGATGAGTGGTCTGAAAGAGaa                >  1:48431/1‑249 (MQ=255)                                                                                                                                                                                                                                      gcaaagcaaaCTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGAGTTTGCAGGAGATGAGTGGTCTGAAAGAGAAAGCgg           >  2:178826/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaagcaaaCTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAGGAGATGAGTGGTCTGAAAGAGAAAGCGGGc         <  1:242741/251‑1 (MQ=255)                                                                                                                                                                                                                                            caaaCTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAGGAGATGAGTGGTCTGAAAGAGAAAGCGGGCGAcc     <  1:128722/251‑1 (MQ=255)                                                                                                                                                                                                                                               aCTCTGACGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAGGAGATGAGTGGTCTGAAAGAGAAAGCGGGCGACCGtt  >  1:95026/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     TCAATCTCTACATGCCTGCCGGGCAGGAGGCAGAGCACCGTTTTGCGATGGGCGAACTGCCAGATGAAGTGCTGGAGATCTGCCAGCGGCTGGCAAAACTCACCGAGATGCTGCGTGGCCTGGCGGAGTTATTTCTTAACGATTTAAGTGAGAAAACCGGCAGCCATGACATTGTACGTCTGCATCGGTTGATTTTGCAGATGAACCGCGCGTTGGGGATGTTCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCGCAACCTTGCGTTCGCTGAACAGTTTTTCGCGTTTGCAGGAGATGAGTGGTCTGAAAGAGAAAGCGGGCGACCGTT  >  NC_000913/834034‑834522

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