Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,688,868	T→C	V98A (GTT→GCT)	manA →	mannose‑6‑phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,688,868	0	T	C	100.0%	192.6 / NA	54	V98A (GTT→GCT)	manA	mannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/25);  total (29/25)

GAGATATCGTTTCACTGCGTGATGTGATTGAGAGTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGTTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGTTTTTGCGCTGACGCCTT  >  NC_000913/1688733‑1688999                                                                                                                                        |                                                                                                                                    gagaTATCGTTTCACTGCGTGATGTGATTGAGAGTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCAt                                                                                                                                 >  2:678748/1‑140 (MQ=255)       tCGTTTCACTGCGTGATGTGATTGAGAGTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCcaaac                                                                                                                           <  2:123640/140‑1 (MQ=255)                tGCGTGATGTGATTGAGAGTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAAt                                                                                                                  >  2:68105/1‑140 (MQ=255)                       tgtgatTGAGAGTGATAAATCGATTCTGCTCGGAGAGGCAGTTGCCAAGCGTTTTGGCGAGCTGCCTTTCATGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGaaa                                                                                                           <  2:446481/140‑1 (MQ=255)                       tgtgatGGAGAGTGATAAATCGACTCTGCTCGGAGAGGCTGTTGCCAAACGCTTTGGGGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGaaa                                                                                                           <  2:305159/140‑1 (MQ=255)                            ttGAGAGTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGt                                                                                                      >  1:318286/1‑140 (MQ=255)                                 agTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGc                                                                                                 >  1:157514/1‑140 (MQ=255)                                 agTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGc                                                                                                 >  2:717555/1‑140 (MQ=255)                                 agTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGc                                                                                                 >  1:47415/1‑140 (MQ=255)                                     aTAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCaaa                                                                                             >  2:796210/1‑140 (MQ=255)                                     aTAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCaaa                                                                                             >  2:425422/1‑140 (MQ=255)                                        aaTCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGGACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGaa                                                                                          <  1:43181/140‑1 (MQ=255)                                        aaTCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGaa                                                                                          <  2:318286/140‑1 (MQ=255)                                        aaTCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGaa                                                                                          <  1:717555/140‑1 (MQ=255)                                        aaTCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGaa                                                                                          <  2:47415/140‑1 (MQ=255)                                         aTCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGaaa                                                                                         <  1:425422/140‑1 (MQ=255)                                               tctGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGaaa                                                                                         >  1:70511/1‑134 (MQ=255)                                               tctGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGaaa                                                                                         <  2:70511/134‑1 (MQ=255)                                               tctGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTACAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGTTCATCCAACCAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCg                                                                                   <  2:856538/140‑1 (MQ=255)                                                  gCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAg                                                                                <  1:678748/140‑1 (MQ=255)                                                    tCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAGTGCCGCAGGt                                                                              >  1:661020/1‑140 (MQ=255)                                                        agagGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATcc                                                                          >  1:81004/1‑140 (MQ=255)                                                            gCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGCATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCgat                                                                      >  2:412011/1‑140 (MQ=255)                                                             ccGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAg                                                                                >  2:69114/1‑129 (MQ=255)                                                             ccGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAg                                                                                <  1:69114/129‑1 (MQ=255)                                                             ccGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATc                                                                           <  1:211515/134‑1 (MQ=255)                                                             ccGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATc                                                                           >  2:211515/1‑134 (MQ=255)                                                             ccGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCgatg                                                                     >  1:713690/1‑140 (MQ=255)                                                                ttGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCa                                                                                 >  2:745335/1‑125 (MQ=255)                                                                ttGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCa                                                                                 <  1:745335/125‑1 (MQ=255)                                                                ttGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCgatggat                                                                  <  1:960222/140‑1 (MQ=255)                                                                    cAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATgccg                                                              <  1:909543/140‑1 (MQ=255)                                                                      aaCGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATgccgcc                                                            <  2:81004/140‑1 (MQ=255)                                                                      aaCGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATgccgcc                                                            >  2:648986/1‑140 (MQ=255)                                                                         gCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAg                                                         <  1:796210/140‑1 (MQ=255)                                                                         gCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAg                                                         <  1:20889/140‑1 (MQ=255)                                                                           tttGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCg                                                       <  2:296656/140‑1 (MQ=255)                                                                           tttGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCg                                                       <  1:68105/140‑1 (MQ=255)                                                                            ttGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGt                                                      <  2:661020/140‑1 (MQ=255)                                                                             tGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTa                                                     >  1:831954/1‑140 (MQ=255)                                                                               gCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAAc                                                   >  2:857392/1‑140 (MQ=255)                                                                                   aCTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATACCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACtata                                               >  2:321776/1‑140 (MQ=255)                                                                                               gTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATaa                                              >  2:190289/1‑129 (MQ=255)                                                                                               gTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATaa                                              <  1:190289/129‑1 (MQ=255)                                                                                                     aGTATTATGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTGACCACAAGc                             >  2:925462/1‑140 (MQ=255)                                                                                                            tGCGCAGCACAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCt                      >  1:207819/1‑140 (MQ=255)                                                                                                                 agcaCAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGttt                 >  2:795371/1‑140 (MQ=255)                                                                                                                 agcaCAGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGttt                 >  2:803481/1‑140 (MQ=255)                                                                                                                   cacaGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGttttt               >  1:850636/1‑140 (MQ=255)                                                                                                                   cacaGCCACTCTCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGttttt               <  1:412011/140‑1 (MQ=255)                                                                                                                           ctctCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGTTTTTGCGCTGAc       >  1:566249/1‑140 (MQ=255)                                                                                                                           ctctCCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGTTTTTGCGCTGAc       >  1:178808/1‑140 (MQ=255)                                                                                                                             ctcCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGTTTTTGCGCTGACGc     >  2:264719/1‑140 (MQ=255)                                                                                                                              tcCATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTGACCACAAGCCGGAGCTGGTTTTTGCGCTGACGcc    <  1:925462/140‑1 (MQ=255)                                                                                                                                cATTCAGGCTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGTTTTTGCGCTGACGCCtt  <  2:557438/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    GAGATATCGTTTCACTGCGTGATGTGATTGAGAGTGATAAATCGACTCTGCTCGGAGAGGCCGTTGCCAAACGCTTTGGCGAACTGCCTTTCCTGTTCAAAGTATTATGCGCAGCACAGCCACTCTCCATTCAGGTTCATCCAAACAAACACAATTCTGAAATCGGTTTTGCCAAAGAAAATGCCGCAGGTATCCCGATGGATGCCGCCGAGCGTAACTATAAAGATCCTAACCACAAGCCGGAGCTGGTTTTTGCGCTGACGCCTT  >  NC_000913/1688733‑1688999

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