Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,113,177	A→G	intergenic (‑116/+257)	rfbB ← / ← galF	dTDP‑glucose 4,6‑dehydratase 1/UTP:glucose‑1‑phosphate uridylyltransferase, low activity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,113,177	0	A	G	100.0%	232.5 / NA	62	intergenic (‑116/+257)	rfbB/galF	dTDP‑glucose 4,6‑dehydratase 1/UTP:glucose‑1‑phosphate uridylyltransferase, low activity
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/34);  total (28/34)

GATTCCCGGTTTCCGGCGTACGTTAATTTATCGACATTAACAACACTATCCTGCGTATTATTTATAATGTGACGAACTACAGCTGAACCAATAAATCCTGCGCCACCAGTAACAAGTATTTTCACTTAATTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCATGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATC  >  NC_000913/2112937‑2113419                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   acattacatactcctgtcaggtgtgatccccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTa                                                                                                                                                                                                                                          <  2:148278‑M2/121‑1 (MQ=255)      acatactcctgtcaggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTatgatg                                                                                                                                                                                                                                     >  1:632677‑M2/126‑251 (MQ=255)                   aggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTAc                                                                                                                                                                                                                        >  2:178018‑M2/113‑251 (MQ=255)                       gtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCAc                                                                                                                                                                                                                    >  2:558294‑M2/109‑251 (MQ=255)                       gtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCAc                                                                                                                                                                                                                    <  1:577449‑M2/143‑1 (MQ=255)                         gataaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTc                                                                                                                                                                                                                  <  2:845812‑M2/145‑1 (MQ=255)                          atcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTcg                                                                                                                                                                                                                 <  2:792064‑M2/146‑1 (MQ=255)                            caccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTcgc                                                                                                                                                                                                                <  1:723080‑M2/147‑1 (MQ=255)                                 cgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTaaa                                                                                                                                                                                                          <  1:712111‑M2/153‑1 (MQ=255)                                       atgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTTTCAGTTTGATTAATAGAGGTAGTAATGCACTCTGCCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCGCTCGCGCTAAACCCGaa                                                                                                                                                                                                    <  1:270455‑M2/159‑1 (MQ=255)                                                            cgcaaataagtatacaaattacttattaaagtataaatagcttaaccatgcttatatgtttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCt                                                                                                                                                                               <  2:864151‑M2/180‑1 (MQ=255)                                                              aaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCAATATTCCTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTa                                                                                                                                                                             <  1:434521‑M2/182‑1 (MQ=255)                                                               aaataagtatacaaagtacttcttaaagtataaatagcttatccatgcttatatgcttacgtctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTACTAGAGGTATTAATGCACGCTACCGCCCCCGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCAACTCGCGCTCAACCCGAAAAATTCATACGCTAATTGTCTTAc                                                                                                                                                                            >  2:420425‑M2/69‑251 (MQ=255)                                                                    agtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATc                                                                                                                                                                       <  1:783451‑M2/188‑1 (MQ=255)                                                                    agtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATc                                                                                                                                                                       <  1:619890‑M2/188‑1 (MQ=255)                                                                        tacacagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATACGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCAACTCGAGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACAAATCCgcc                                                                                                                                                                   >  2:246720‑M2/60‑250 (MQ=255)                                                                        tacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGct                                                                                                                                                                   >  1:60307‑M2/60‑251 (MQ=255)                                                                          caaagtacttattaacgaatacatagcttatccatacttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTAAAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACACGAAAAATTAAAACGCTAATTGTCTTACCAATccgccct                                                                                                                                                                 >  2:668711‑M2/58‑248 (MQ=255)                                                                                    attaaagtacaaatagcttatcaatgcttatatgcttacggctttatTTATTCAATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAgg                                                                                                                                                       >  2:666827‑M2/48‑251 (MQ=255)                                                                                          gtatacatagcttatccctgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGACAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATc                                                                                                                                                 >  2:457240‑M2/42‑251 (MQ=255)                                                                                              aaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAACATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTgg                                                                                                                                             >  2:260906‑M2/38‑251 (MQ=255)                                                                                                atagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGa                                                                                                                                            >  2:170668‑M2/36‑250 (MQ=255)                                                                                                          catgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAACCCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGCAAACTTTGAAt                                                                                                                                 >  2:401520‑M2/26‑251 (MQ=255)                                                                                                             gcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATaaa                                                                                                                              <  1:503446‑M2/229‑1 (MQ=255)                                                                                                             gcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATaaa                                                                                                                              <  1:41920‑M2/229‑1 (MQ=255)                                                                                                              cttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATaaa                                                                                                                              >  1:350033‑M2/22‑250 (MQ=255)                                                                                                                 atatgcttacgtctttatTTATTCCATATTTCTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAAGGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAATccc                                                                                                                          <  2:354721‑M2/233‑1 (MQ=255)                                                                                                                  tatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTGCCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                                      >  1:229251‑M2/18‑238 (MQ=255)                                                                                                                  tatgcttacggctttatTTATTCCATATTACTTCAGAGCAGGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTGCCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                                      <  2:229251‑M2/221‑1 (MQ=255)                                                                                                                            gctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGTTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACt                                                                                                               <  1:339‑M2/244‑1 (MQ=255)                                                                                                                             ctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTc                                                                                                              >  1:33327‑M2/7‑251 (MQ=255)                                                                                                                             atttattttgtccttattatTTCAGAGTTTGATGTGAAATAGGCGGCTCTAAGTTTAACTAATAGAGGTAGTAATGCACCCTACCGCCCCTGGCTTTACAGATACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCTCTCGCGCTAAACCCGAAAAATTCAAACGTTAATTGTCTTTCCAATCCGCTCTGGAAACAAGGAAATTCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTc                                                                                                              <  2:685704/237‑1 (MQ=255)                                                                                                                               ttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCg                                                                                                             <  1:558294‑M2/246‑1 (MQ=255)                                                                                                                                    tATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCgttgtt                                                                                                        >  1:880791/1‑250 (MQ=255)                                                                                                                                        ccATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCt                                                                                                   >  1:555122/1‑251 (MQ=255)                                                                                                                                          atatTACTTCAGAGCCTGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGCAACAAGTAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTg                                                                                                  >  2:551818/1‑250 (MQ=255)                                                                                                                                                 ttCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATcc                                                                                                                                                                      >  2:56689/1‑175 (MQ=255)                                                                                                                                                 ttCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATcc                                                                                                                                                                      <  1:56689/175‑1 (MQ=255)                                                                                                                                                      agCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACCCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATaa                                                                                     <  2:203774/251‑1 (MQ=255)                                                                                                                                                           gCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAaaca                                                                                 >  2:437403/1‑250 (MQ=255)                                                                                                                                                              gtgAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAaacaac                                                                               <  1:666827/249‑1 (MQ=255)                                                                                                                                                                       gCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCAACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGa                                                                        >  1:350578/1‑247 (MQ=255)                                                                                                                                                                       gCGGCTCTCAGTTTGATTAATAGAGGTATTAAGGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGAt                                                                    <  2:373724/251‑1 (MQ=255)                                                                                                                                                                         ggCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGAt                                                                    >  1:309382/1‑249 (MQ=255)                                                                                                                                                                         ggCTCTCAGTTTAATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGAt                                                                    <  2:309382/249‑1 (MQ=255)                                                                                                                                                                               cggTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGcaa                                                              <  2:22451/247‑1 (MQ=255)                                                                                                                                                                               cAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTCTTCTGATGGTTTATATAAAACAACGGCAGGaatgtcctcaaa                                                             >  2:862291/1‑240 (MQ=255)                                                                                                                                                                                     gATTAATAGTGGTATTAAGGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTGGCAACAAATTa                                                       <  2:628185/250‑1 (MQ=255)                                                                                                                                                                                     gATTAATAGAGGTATCAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTAc                                                      >  1:86473/1‑251 (MQ=255)                                                                                                                                                                                          aTAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtg                                                 <  1:437403/251‑1 (MQ=255)                                                                                                                                                                                             gagGTATTAATGCCCGCTACCGCCCCTGGCGTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGAGTCGCAACAAATTACTTTtgctg                                              <  2:880791/251‑1 (MQ=255)                                                                                                                                                                                                    tAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAAttt                                       >  1:354050/1‑251 (MQ=255)                                                                                                                                                                                                        gCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGAGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAAT‑ACTTTTGCTGCGAATTTTCACt                                  <  1:527929/251‑1 (MQ=255)                                                                                                                                                                                                             cTACCGCCCCTGGCTTTACAGCTCCCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGCAATTTCACTGCCg                              >  2:575665/1‑251 (MQ=255)                                                                                                                                                                                                               aCCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGtt                            <  1:178018/251‑1 (MQ=255)                                                                                                                                                                                                                 cGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGCATTTTCACTGCCGTtat                          <  1:173328/251‑1 (MQ=255)                                                                                                                                                                                                                 cGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTtat                          <  1:346087/251‑1 (MQ=255)                                                                                                                                                                                                                    cccTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAAt                       <  1:260906/251‑1 (MQ=255)                                                                                                                                                                                                                         gCTTTACAGCTACCATAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTAACCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGACACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTACCTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACTAATTCCTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCt                  <  1:659097/251‑1 (MQ=255)                                                                                                                                                                                                                                aGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAAc           >  1:383980/1‑251 (MQ=255)                                                                                                                                                                                                                                    aCCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCg         <  1:45829/249‑1 (MQ=255)                                                                                                                                                                                                                                         agCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATc  >  1:443982/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   GATTCCCGGTTTCCGGCGTACGTTAATTTATCGACATTAACAACACTATCCTGCGTATTATTTATAATGTGACGAACTACAGCTGAACCAATAAATCCTGCGCCACCAGTAACAAGTATTTTCACTTAATTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCATGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCTTATCAACCGTTACATC  >  NC_000913/2112937‑2113419

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