Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,193,996	T→G	D70A (GAT→GCT)	thiE ←	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,193,996	0	T	G	100.0%	191.5 / NA	59	D70A (GAT→GCT)	thiE	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (26/33);  total (26/33)

GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCT  >  NC_000913/4193751‑4194240                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     gCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGtt                                                                                                                                                                                                                                                  <  1:104983/250‑1 (MQ=255) gCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTg                                                                                                                                                                                                                                                 <  1:171453/251‑1 (MQ=255)  ccGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAATGCATCTGTTTGTTTTGCGTCGGGAACACATGTCCCAGCGT‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGa                                                                                                                                                                                                                                                <  2:280472/251‑1 (MQ=255)    gTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGAt                                                                                                                                                                                                                                               >  1:356494/1‑250 (MQ=255)      gCCAGCTGTTCCAGACCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATaa                                                                                                                                                                                                                                             >  1:127248/1‑250 (MQ=255)             gTTCCAGCCCCTGCGGTGCAGAAGGCAGCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcgc                                                                                                                                                                                                                                     <  1:236133/251‑1 (MQ=255)                    ccccTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTtat                                                                                                                                                                                                                              >  1:334537/1‑251 (MQ=255)                         gCGGTGCAGAAGGCATCTGTTTGGTTGGCGTCGGGAACACTTGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCTAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcg                                                                                                                                                                                                                          <  1:14856/250‑1 (MQ=255)                          cGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGCTTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGACCCCAATGGACGCAATACGCCTGATGCTTGATCGCCAGCCGCCAGTACGCGTTGATAAACAATCGCGCGTtatcacggc                                                                                                                                                                                                                        >  1:264824/1‑245 (MQ=255)                            gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGcggcgg                                                                                                                                                                                                                      >  2:106732/1‑251 (MQ=255)                               cAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGccc                                                                                                                                                                                                                   <  1:34781/251‑1 (MQ=255)                                  aaGGCATCTGTGTGGTTTGCGCCGGGAACACATTTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACTTCGATTTCCATATCGTCATGTGTCGAAACTCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCGGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAgc                                                                                                                                                                                                                <  1:70682/251‑1 (MQ=255)                                   agtcATCTGTTTGGTTTGCGCCGGGACCCCATGTCCCA‑CGCCGATATAAGAGGCGCGTGCTCCCGGCGCTACGCCGATTGCCATTTCGGCATGTGTCGAAACGCCCGGCCGCAGGCCTGCCGCGCGGATGGAATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATCGACTCCATACGTCTGATTCTTGAGCGCCAGCCGCCAGTAAGCGTTGATAACCAATCGCGCGTTGTAGCGGCGGACCAgcg                                                                                                                                                                                                               <  2:208669/248‑1 (MQ=255)                                   aGGCAGCTGTGTGGTTTGCGTCGGGAACTCATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATACCGTCATGTGTCGAAACGCCCCGCCTCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAgcg                                                                                                                                                                                                               <  2:369457/251‑1 (MQ=255)                                     gCATCTGTTGGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGGCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCa                                                                                                                                                                                                             <  1:75930/251‑1 (MQ=255)                                          tGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgc                                                                                                                                                                                                        >  1:306847/1‑251 (MQ=255)                                           gtttggtttgCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgcc                                                                                                                                                                                                       >  1:338889/1‑251 (MQ=255)                                                         gggAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCCCGCCATCGGc                                                                                                                                                                                         >  2:348185/1‑251 (MQ=255)                                                           gAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCGGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAACCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGTTGCTTGACCTCCAGCGGCCAGGAAGCGTTGATAAACAATCGCTCGTTATAGCGGCGGCCAAGCGCAATTGCCGCCACGACATCGGCtt                                                                                                                                                                                       <  2:194612/251‑1 (MQ=255)                                                                caTGCCCCAGCGC‑GATATAAGGGGGGGGTGCTCCCAGCGCGGCGCCGATTTCCATATCGTCATGGGCCGAAACGCCCAGCCGCAGGCCGGCCGCGGGGATGGCATTGAGATCGGGGGCTTGCAAAT‑CTCCCTGCCCAAAAGGGACGCCATACGCCGGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAcc                                                                                                                                                                                  <  2:289445/251‑1 (MQ=255)                                                                   ttCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACctc                                                                                                                                                                                <  1:249368/249‑1 (MQ=255)                                                                      cccagctggatatagtgggggCGTGCGGCCAGCGCGACGTCGCTTCCCATACCGTAATGGGCCGAAACGCCCGGCCGACGGCCTGCCGCGCGGAGGGCATTGAGATCGGTGGTTTGCAAAT‑CTTCCTGCCCCAAAGGGACGCCATACGCCTGTAGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGTCACGACATCGGCTTCCACCTCTTc                                                                                                                                                                             <  2:338124/234‑1 (MQ=255)                                                                        aGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCAATTTCCATATCGTCATGTGCCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATc                                                                                                                                                                          <  1:326645/251‑1 (MQ=255)                                                                                   aaGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGCCATCGGCTTCCACCCCTTCATCTCGCCGATCt                                                                                                                                                                >  2:13839/1‑251 (MQ=255)                                                                                       ggggcggtctgcCGGCGCGCCGTCGTTTTCCATATCGTCATGTGCCGAAACGCCCAGCCGCAGGCCTGCCGGGCGGATGGCATTGAGATCGGTGGCTTCCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCTTGAGGTTTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGa                                                                                                                                                            <  2:378047/243‑1 (MQ=255)                                                                                       ggggCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATGGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGa                                                                                                                                                            <  1:434232/251‑1 (MQ=255)                                                                                                  ccAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTgg                                                                                                                                                 <  1:360201/251‑1 (MQ=255)                                                                                                    aGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCAGACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGgag                                                                                                                                               <  1:346375/251‑1 (MQ=255)                                                                                                        cgACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGCTCTTTGATGCGTAGCTGGAGAgt                                                                                                                                            >  2:341882/1‑250 (MQ=255)                                                                                                                  ttCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATTGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTg                                                                                                                                 >  1:151847/1‑251 (MQ=255)                                                                                                                    ccATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTCCGTACGCCTGc                                                                                                                                >  2:88487/1‑250 (MQ=255)                                                                                                                           gTCATGTGTCGAACCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATCGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGCGAGTACGTTCGCCTGCATcccac                                                                                                                         >  2:27007/1‑247 (MQ=255)                                                                                                                                gtgtCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGt                                                                                                                   >  1:209163/1‑251 (MQ=255)                                                                                                                                       aacccccAGCCCCGGGCCCGCCACGCGGAGGGCATTAAGATCGGGGGCTGACAAAT‑CTTCCTGCCCATAATGGCCGCCATACGCCTGATGCTTGATCGCCACCCGCCGGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTTCGCCTGCTTCCAACAGACGGTCGATcc                                                                                                            <  2:387140/247‑1 (MQ=255)                                                                                                                                         cGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAc                                                                                                          <  1:333537/251‑1 (MQ=255)                                                                                                                                          gCCCAGCCGCAGGCCAGCCGCGCGGATGGCATTGTGATAGGTGGCTTGCAAAT‑CTTCCTGCCCCAACTGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCCGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCAATTGCCGCCACGACACCGGCTTCCACCTCTGCCTCGCGCCGATCATTCATGCGTAGCTGGAGAGTACGTCCGCCTGCCTCCAACAGACGATCGATCCACt                                                                                                         >  2:279081/1‑251 (MQ=255)                                                                                                                                               gctgCAGCCCTGCCGCGCGGATGGCATTGAAAGCGGGGGATTGCAGAC‑TTTCCGGCCCCAAAGGTACGCCAACCGCCGGATGTTTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTGCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACg                                                                                                    <  2:130752/248‑1 (MQ=255)                                                                                                                                                   cAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGt                                                                                                >  2:59218/1‑251 (MQ=255)                                                                                                                                                     ggCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTcc                                                                                              <  1:106732/251‑1 (MQ=255)                                                                                                                                                            cgcgcgGATGGCATTGAGCTCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGAACCACTGTACGCTGTCCACCACCg                                                                                       >  2:372487/1‑251 (MQ=255)                                                                                                                                                                   aTGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATGCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAg                                                                                <  1:202934/251‑1 (MQ=255)                                                                                                                                                                        aTTGAGATCGGTGGCTTGCACAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGCTCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCCACAGACGTTCGCTCCACTGGACGCTGTCCACCACCGGGGCCAGTCCTg                                                                           >  2:29417/1‑251 (MQ=255)                                                                                                                                                                                 ggTGGCT‑GCAAATGCTTCCTGCCCCAAATGGTCGCCATACGCCGTATGCTTGATCCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTTGAGCGGCGGCCTAGCGCAATGGCCGCCAGGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTATCTTGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAg                                                                  <  1:35084/251‑1 (MQ=255)                                                                                                                                                                                    ggCTTGCAAAT‑CTTCCTGCCCCACATGGAAGCCATCCGCCTGATGCTTGCTCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACCTCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATACGTCGCTGGAGAGTAAGTACGCCTGCCTCCCACAGACGTTCGATCCCCTGTACGATGTCCACCACCGGGTACAGTCCTGAACgcacaggg                                                                >  2:399212/1‑243 (MQ=255)                                                                                                                                                                                      cTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTAc                                                              >  1:142054/1‑250 (MQ=255)                                                                                                                                                                                        tGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACag                                                            <  1:392235/250‑1 (MQ=255)                                                                                                                                                                                        tGCAAAT‑CGTCCTGCCCCAAATGGACGCCATACGCAGGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACag                                                            <  1:192756/250‑1 (MQ=255)                                                                                                                                                                                          gcaatcctcccTGCCCCAATGGGGCACCATACGCGTGTTGCTGGATCGCCAGCCGCCAGTAAGCGTGGATAAACAAGCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCGGGAGATTGCGTACGCCTGCATCCAACAGACGTTTGATCCACTGAACGCTGTCCACCACCGGGTACGGTCCTGTACGAAAATGTACagg                                                           <  2:267652/241‑1 (MQ=255)                                                                                                                                                                                                cgtCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGaaaa                                                    <  1:222968/249‑1 (MQ=255)                                                                                                                                                                                                    cTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATc                                                  >  2:291011/1‑249 (MQ=255)                                                                                                                                                                                                      gCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCTCCACGACATCGGCTTCCGCCTCGTCATCGCTCCGATCTGTGATGCGTAGCTGGAGAGTGCGTGCGCCTGCATCCAACAGACGGTCGATCCACTGTTCTCTGTCCACCGCCGGGTATAGTCCTGAACGAATAGGTACAGTAGGAAACTCAgg                                               <  2:218359/250‑1 (MQ=255)                                                                                                                                                                                                          cAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGAt                                          <  1:341882/251‑1 (MQ=255)                                                                                                                                                                                                           aaaTGGACGCCATACGCCTGATGCTTGCTCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATa                                         <  1:291011/251‑1 (MQ=255)                                                                                                                                                                                                              tGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGGCCCCAACCGGGTACAGTCCTGAACGAAAAGGTACAGGCGGAAAATCAGGCTGATACAt                                      >  1:178081/1‑251 (MQ=255)                                                                                                                                                                                                                                     gATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGTCGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGAACAGGAGGACAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGt                >  2:324136/1‑250 (MQ=255)                                                                                                                                                                                                                                     gATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGt                >  1:53186/1‑250 (MQ=255)                                                                                                                                                                                                                                           cAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGAACAGTCCTGAACGAAAAGGTACAGGCGGAAAATCAAGCTGATACATCACGCTTCCTCCTTACGACGGAAGATTtc         >  2:69864/1‑251 (MQ=255)                                                                                                                                                                                                                                            aGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGCACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTtct        >  2:27476/1‑251 (MQ=255)                                                                                                                                                                                                                                                  cAGTAAGCGTTGAAAAACATGCGCGCGTTATCGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTACGctt  <  2:24394/251‑3 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGC‑GATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAT‑CTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCT  >  NC_000913/4193751‑4194240

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