Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	958,173	A→G	D174G (GAC→GGC)	serC →	phosphoserine/phosphohydroxythreonine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	958,173	0	A	G	100.0%	208.3 / NA	59	D174G (GAC→GGC)	serC	phosphoserine/phosphohydroxythreonine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/28);  total (31/28)

ACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTAAATGGCT  >  NC_000913/957926‑958419                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       accaccGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAAGGAAACCATCGATGGTATCGCCATCGACGAAACGCCGGACTTCGTCGCAGATGTGGTGGTCGCCGCTGGCtt                                                                                                                                                                                                                                                     <  2:40932/251‑1 (MQ=255)      cGCAGATTATGTCCCTGCAGGTTACTGGGCGGCAAGCGACTTCACAGAAGCGCAAAAACACTGCACGCCTAATGTCCTTGCCGCCACAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGACTGCCAACTCTATGATAATGCTGCTTATCTGCATTATTGCCCGACTGACACAATCGATGGTCTCGCCATCGCCGAAACGACAGACTTCGGCGAAGATGTGGTGGTCGCCGCTTGCTtctcat                                                                                                                                                                                                                                                >  2:219374/1‑249 (MQ=255)        cAGATTATGTTGATGCCGGCTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTAGATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGTCTTCGGCGCAGATGTGGTGGTCGACGCTGGCTTCTCTTCa                                                                                                                                                                                                                                              <  2:37403/251‑1 (MQ=255)         aGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCCTTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCa                                                                                                                                                                                                                                              >  2:150225/1‑250 (MQ=255)              aTGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCAtt                                                                                                                                                                                                                                        <  1:18201/251‑1 (MQ=255)                gTTGATGCCGATTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTc                                                                                                                                                                                                                                       <  1:260273/250‑1 (MQ=255)                             aCTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCGTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCGTATATGCATTATTGCCCGAATGAAACCATCGATGGGATCGCCATCGACGAAACGCCAGGCTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGAtt                                                                                                                                                                                                                         <  1:469631/251‑1 (MQ=255)                                gggCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGAc                                                                                                                                                                                                                      <  1:292789/251‑1 (MQ=255)                                       aaGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGcc                                                                                                                                                                                                               >  2:381610/1‑251 (MQ=255)                                          tGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGtt                                                                                                                                                                                                            >  1:329474/1‑251 (MQ=255)                                               ttAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGAc                                                                                                                                                                                                                      >  2:341115/1‑236 (MQ=255)                                               ttAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGAc                                                                                                                                                                                                                      <  1:341115/236‑1 (MQ=255)                                                tAAAGAAGCGAAAAAATACTGCACGCCTCATGTCTTTGACGCCACAGTGCCTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGTTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGCTGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTTTGgtg                                                                                                                                                                                                      >  2:461284/1‑251 (MQ=255)                                                        cGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTAc                                                                                                                                                                                              >  1:412448/1‑251 (MQ=255)                                                        cGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTAc                                                                                                                                                                                              <  2:157659/251‑1 (MQ=255)                                                              aaTACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGgc                                                                                                                                                                                        >  1:178750/1‑251 (MQ=255)                                                              aaTACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGgc                                                                                                                                                                                        >  1:163736/1‑251 (MQ=255)                                                                    gCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTGAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAATCGCCAGACTTCGGCGCAGATGTGGTGTTCGCCGCTGGCTTCTCGTCAACCATTCTGTCCCTTCCTATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCa                                                                                                                                                                                   <  1:426272/250‑1 (MQ=255)                                                                      aCGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGa                                                                                                                                                                                 >  1:33562/1‑250 (MQ=255)                                                                          cTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAat                                                                                                                                                                            >  1:323957/1‑251 (MQ=255)                                                                            aaTGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAatat                                                                                                                                                                          >  1:127311/1‑251 (MQ=255)                                                                            aaTGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAatat                                                                                                                                                                          >  2:477158/1‑251 (MQ=255)                                                                                        gCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTgg                                                                                                                                                              >  2:121672/1‑251 (MQ=255)                                                                                                    gTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTGCCCGTTCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGTCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATc                                                                                                                                                    <  1:216718/249‑1 (MQ=255)                                                                                                    gTTGATGGGCTGCGCTCGGTGAAGCCAAGGAGCGAATGGCAACTGTCTGATAATGCAGATTATATGCATTAATGCCCGAATGAACCCATCGACGGTATCGCCATCGACGAAACGACAGCCTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACCCTGGCGCGCAGAAAAATATCGGCCCGGCTGGccggtaaatcgt                                                                                                                                                  >  2:312143/1‑241 (MQ=255)                                                                                                      tGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCTGCTGGCCTGACAATCGTc                                                                                                                                                 <  1:383515/250‑1 (MQ=255)                                                                                                               gcgcgcGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGAc                                                                                                                                                        >  2:166191/1‑234 (MQ=255)                                                                                                               gcgcgcGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGAc                                                                                                                                                        <  1:166191/234‑1 (MQ=255)                                                                                                                           gCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTgg                                                                                                                           <  2:101063/251‑1 (MQ=255)                                                                                                                            ccAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTggg                                                                                                                          <  1:477158/251‑1 (MQ=255)                                                                                                                                      aaTGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAat                                                                                                                >  1:248307/1‑251 (MQ=255)                                                                                                                                            aaCTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCCTTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATcgcg                                                                                                          >  2:230504/1‑251 (MQ=255)                                                                                                                                                    gATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTc                                                                                                  >  1:97605/1‑251 (MQ=255)                                                                                                                                                     aTAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCGTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCTGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCg                                                                                                 <  1:230504/251‑1 (MQ=255)                                                                                                                                                        aTGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTCTGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGAtt                                                                                              >  1:444825/1‑251 (MQ=255)                                                                                                                                                          gctgctTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCAGCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCt                                                                                            <  2:329474/251‑1 (MQ=255)                                                                                                                                                           ctgctTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTg                                                                                           >  1:240998/1‑251 (MQ=255)                                                                                                                                                             gctTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTAGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCTTGTCCGTCGATTCTg                                                                                           <  1:375351/249‑1 (MQ=255)                                                                                                                                                                 atatGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGattat                                                                                     <  1:25841/251‑1 (MQ=255)                                                                                                                                                                 atatGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAAGATCGGCCCGGCTGGCCTGACAATCGGCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGatta                                                                                      <  1:461284/250‑1 (MQ=255)                                                                                                                                                                   aggCATTATTGCCCGAATGAAACCATCGAGGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGTTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTc                                                                                   <  2:76460/249‑1 (MQ=255)                                                                                                                                                                         tattGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCTGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATcc                                                                              <  1:101540/250‑1 (MQ=255)                                                                                                                                                                         tattGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCt                                                                             >  1:99833/1‑251 (MQ=255)                                                                                                                                                                         tattGCCCGAAGGAAACCAACGGTGGAATCGACACCGACGTACCGCCAGCCTTCGGCGCGGTTGGGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCTTCCGAGTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCATGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATGCCATcc                                                                              <  2:272920/250‑1 (MQ=255)                                                                                                                                                                          attGCCCGAATGAAACCCTCGATGGTATCGCCATCGAAGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTc                                                                            >  1:193583/1‑251 (MQ=255)                                                                                                                                                                          attGCCCGAAGGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACTTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTc                                                                            <  2:323957/251‑1 (MQ=255)                                                                                                                                                                                 gAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGAt                                                                      >  1:220159/1‑250 (MQ=255)                                                                                                                                                                                  aaTGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTACGATTGACGTCAGCCGTTATGGTGAAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTAGTGAAGGTTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTCTTCCATCCTCAACGATa                                                                     >  1:469718/1‑250 (MQ=255)                                                                                                                                                                                             cGATGGTATCGCCATCGACGAAACGCCATACTTCGGCGCAGATGTGGTGGGCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCTCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATg                                                          <  1:463578/250‑1 (MQ=255)                                                                                                                                                                                             cGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATg                                                          <  1:334043/250‑1 (MQ=255)                                                                                                                                                                                                    aTCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAacac                                                  >  1:13024/1‑251 (MQ=255)                                                                                                                                                                                                       gCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACAcgcc                                               <  1:121672/251‑1 (MQ=255)                                                                                                                                                                                                         cATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACAcgccgc                                             >  1:417651/1‑251 (MQ=255)                                                                                                                                                                                                                  aaCGCCAGACTTCGGCGCAGATGTGGCGGTCGCCGCTGGCTTCTCTTCAACCCTTCTTTCCCGTCCGGTTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCAGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGACGCCGACCTTTg                                    >  1:117390/1‑251 (MQ=255)                                                                                                                                                                                                                    cGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGGGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGcc                                  >  1:392463/1‑251 (MQ=255)                                                                                                                                                                                                                                     gATGTGGTGGTCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtct                 <  2:97605/251‑1 (MQ=255)                                                                                                                                                                                                                                           gtggtCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGAc                                                                                                                                                        <  2:264411/110‑1 (MQ=255)                                                                                                                                                                                                                                           gtggtCGCCGCTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGAc                                                                                                                                                        >  1:264411/1‑110 (MQ=255)                                                                                                                                                                                                                                                    gcTGGCTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTAAATGGCt  >  2:446421/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ACCACCGCAGATTATGTTGATGCCGGTTACTGGGCGGCAAGTGCCATTAAAGAAGCGAAAAAATACTGCACGCCTAATGTCTTTGACGCCAAAGTGACTGTTGATGGTCTGCGCGCGGTTAAGCCAATGCGTGAATGGCAACTCTCTGATAATGCTGCTTATATGCATTATTGCCCGAATGAAACCATCGATGGTATCGCCATCGACGAAACGCCAGACTTCGGCGCAGATGTGGTGGTCGCCGCTGACTTCTCTTCAACCATTCTTTCCCGTCCGATTGACGTCAGCCGTTATGGTGTAATTTACGCTGGCGCGCAGAAAAATATCGGCCCGGCTGGCCTGACAATCGTCATCGTTCGTGAAGATTTGCTGGGCAAAGCGAATATCGCGTGTCCGTCGATTCTGGATTATTCCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTAAATGGCT  >  NC_000913/957926‑958419

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