Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,738,635	(A)6→5	coding (1139/2031 nt)	malS →	alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,738,635	0	A	.	100.0%	151.9 / NA	34	coding (1139/2031 nt)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/15);  total (19/15)

CGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCTATAAAAACAAAATGGATACCCACGCCAAAGCCATTGACGGCTATACGCCGCGCGATTACTTAACCCACTGGT‑TAAGTCAGTGGGTCCGCGACTATGGG  >  NC_000913/3738403‑3738855                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              cgcgGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATc                                                                                                                                                                                                            >  1:169966/1‑251 (MQ=255) cgcgGTATTCGTATTCTCGTTGATGTCGTGGTGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATc                                                                                                                                                                                                            <  1:250738/251‑1 (MQ=255)     gtattcgtattcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATCTCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCg                                                                                                                                                                                                        >  1:126362/1‑251 (MQ=255)     gtattcgtattcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCg                                                                                                                                                                                                        <  2:169966/251‑1 (MQ=255)             attcTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCTGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGa                                                                                                                                                                                                 <  2:262002/250‑1 (MQ=255)                  ctgtGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGATCGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGGTTTAACGATTACATGAAGTGCAGCGACGAAACAGGCTGGGATAAATGGTTGGGAAAAA‑CTGGAGCAGAACGGATATCGGCGATTACGACAAt                                                                                                                                                                                             <  2:60066/246‑1 (MQ=255)                  ctTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAAt                                                                                                                                                                                             <  2:46540/249‑1 (MQ=255)                      gATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGAtt                                                                                                                                                                                                     <  1:107782/237‑1 (MQ=255)                      gATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGAtt                                                                                                                                                                                                     >  2:107782/1‑237 (MQ=255)                         gTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTTCCGCCGGTCAAACCTGGCATAGCTTTAACGATTACATTAATTGCAGCGACAAAACAGGCTGGGATAAATGGTTGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTc                                                                                                                                                                                    <  1:307960/251‑1 (MQ=255)                            gTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGACGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTAAGACAATACTGGATTcgac                                                                                                                                                                                 >  2:264906/1‑251 (MQ=255)                              gatgaACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAAGTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATAGGCGCTTAAGGCAATCCTGGATTcgaagg                                                                                                                                                                               >  2:145322/1‑248 (MQ=255)                                                ggCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGc                                                                                                                                                             <  2:126362/250‑1 (MQ=255)                                                              atatGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGata                                                                                                                                               <  1:145322/251‑1 (MQ=255)                                                                    aGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAACTGAAATAAAAGCCGGGAGAACGCTGGAGCGACTGGACACCTGCCGCCAGGCAAACCTGGCATAGCCTTCACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATCAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGCAAATCCAGGCTTCGACGATCACACTATGACGCCAGCCTTTTTGCCGGATATCaaa                                                                                                                                          >  1:48378/1‑250 (MQ=255)                                                                        gTATCAGTTTGGCGCGTTATATCTTTATGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTAGCTAGCCTTTTTGCCGGATATCAAAACCg                                                                                                                                      >  2:213661/1‑250 (MQ=255)                                                                                  ggCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATATCTTTAACGATTACATTAATTTCAGCTACAAAACAGGCTGGGATAAATGGTTGGGAAAAA‑CTGGATCAGATCGGATATCGGCGATTACGACAATTCTGGAGTCGACGATCTCACGAGGTCGCTAGCTTGTTTGCCGGATATCAAAACCGAATCAACTAcc                                                                                                                           <  1:13354/251‑1 (MQ=255)                                                                                          atatCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTc                                                                                                                                                                                    <  1:118039/186‑1 (MQ=255)                                                                                          atatCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTc                                                                                                                                                                                    >  2:118039/1‑186 (MQ=255)                                                                                              cTTTCTGGTGACGCAGAGAAAAAATCGCTGGGTGAACGCCGCAGCGACTGGACACCTGCCGCCGGGCAAACCTGGCATAGCTTTCACGATTACATTACTTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACACTCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTg                                                                                                               >  2:307113/1‑251 (MQ=255)                                                                                                     gTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGTGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGc                                                                                                                                                                                                         >  1:133168/1‑154 (MQ=255)                                                                                                     gTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTGTAACGATTACATTAATTTCAGCGACAAAACATGCTGGGATAAAGTGTGGGGAAAAA‑CTTGATCAGAACGGATATCGGc                                                                                                                                                                                                         <  2:133168/154‑1 (MQ=255)                                                                                                               gAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCTATaaaa                                                                                               >  1:183344/1‑250 (MQ=255)                                                                                                                          tGGGTGAACGCTGGAGAGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCTATAAAAACAAAATGGAt                                                                                    >  2:31117/1‑250 (MQ=255)                                                                                                                                                                       ggCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGATGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTTCCGGTGTTCTATAAAAACAAAATGGCTACCCACGCCAAAGCCATTGACGGCTATACGCCGCGCGATTACTTa                                       >  1:125497/1‑250 (MQ=255)                                                                                                                                                                            aGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCTATAAAAACAAAATGGATACCCACGCCAAAGCCATTGACGGCTATACGCCGCGCGATTACTTAACCCAc                                 >  1:73492/1‑251 (MQ=255)                                                                                                                                                                                    cGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCTATAAAAACAAAATGGATACCCACGCCAAAGCCATTGACGGCTATACGCCGCGCGATTACTTAACCCACTGGT‑Ta                          >  1:53624/1‑249 (MQ=255)                                                                                                                                                                                       ttACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAAcc                                                                                                                                       <  2:243307/138‑1 (MQ=255)                                                                                                                                                                                       ttACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAAcc                                                                                                                                       >  1:243307/1‑138 (MQ=255)                                                                                                                                                                                               aTTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAA‑CTGGATCAGAAAGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCTATAAAAACAAAATGGATACCCACGCCAAAGCCATTGACGGCTATACGCCGCGCGATTACTTAACCCACTGGT‑TAAGTCAGTGGGTcc             >  1:175634/1‑251 (MQ=255)                                                                                                                                                                                                     ccGACAAAACAGGCTGTGATAAATGGTTGGGAAAAA‑CTGGATCAGAACGGAGATCGGCGATTACGACAAGCCTGGATTCGACGATCTCACTATGTCGCTAGCCGTTTTGCCGGATATCAAAACCTAAGCAACTTCCGCTTCTGGGCTGCTGGTGTTCTATAAAAACAAAATGGATTCCCTCGCCATAGCCAGTGACGGCTATACGCCGCGCGAGTTCGTAGCCCTC‑GGTGTAAGTCAGTGGGTCCGCGACt       <  2:131422/250‑1 (MQ=255)                                                                                                                                                                                                           ataCAGGCTGGGATAAATGGAGGGGAAAAA‑CTGGCTCAGCACGGATATCGGCGATTACGACCATACAGGATTCGACGATATCACTTTGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTATGCCGGTGTTCTATAAAAACAAAATGGATACCCACGCCCAAGCCCTTGACGGCTCTCCGCCGCGCGATTACTTACCCCACTGGG‑TAAGTCCGTGGGTCCGCGACTATggg  >  2:51017/3‑250 (MQ=255)                                                                                                                                                                                                                         aaaTGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCt                                                                                                     <  2:105685/138‑1 (MQ=255)                                                                                                                                                                                                                         aaaTGGTGGGGAAAAA‑CTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCt                                                                                                     >  1:105685/1‑138 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              CGCGGTATTCGTATTCTCTTTGATGTCGTGATGAACCACACCGGCTATGCCACGCTGGCGGATATGCAGGAGTATCAGTTTGGCGCGTTATATCTTTCTGGTGACGAAGTGAAAAAATCGCTGGGTGAACGCTGGAGCGACTGGAAACCTGCCGCCGGGCAAACCTGGCATAGCTTTAACGATTACATTAATTTCAGCGACAAAACAGGCTGGGATAAATGGTGGGGAAAAAACTGGATCAGAACGGATATCGGCGATTACGACAATCCTGGATTCGACGATCTCACTATGTCGCTAGCCTTTTTGCCGGATATCAAAACCGAATCAACTACCGCTTCTGGTCTGCCGGTGTTCTATAAAAACAAAATGGATACCCACGCCAAAGCCATTGACGGCTATACGCCGCGCGATTACTTAACCCACTGGT‑TAAGTCAGTGGGTCCGCGACTATGGG  >  NC_000913/3738403‑3738855

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