Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,371,482 (AT)4→3 coding (1267‑1268/1680 nt) ycjM → alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,371,4760A.100.0% 150.2 / NA 36coding (1261/1680 nt)ycjMalpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (20/16);  total (20/16)
*NC_0009131,371,4770T.100.0% 150.2 / NA 36coding (1262/1680 nt)ycjMalpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (20/16);  total (20/16)

GTCCGTATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGATATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAACCCGAGAACTGAACGATGAAG  >  NC_000913/1371351‑1371602
                                                                                                                             ||                                                                                                                             
gTCCGTATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGat                                                                                                                               >  2:471252/1‑127 (MQ=255)
gTCCGTATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGat                                                                                                                               <  1:471252/127‑1 (MQ=255)
  ccGTATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaaag                                                                                                                     >  1:166568/1‑129 (MQ=255)
  ccGTATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaaag                                                                                                                     <  2:166568/135‑7 (MQ=255)
     tATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGa                                                                                                                                <  1:121921/121‑1 (MQ=255)
     tATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGa                                                                                                                                >  2:121921/1‑121 (MQ=255)
         aaataaatGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCt                                                                                                        <  2:492557/139‑1 (MQ=255)
                        atatGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaaag                                                                                                                     >  2:253173/1‑107 (MQ=255)
                        atatGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaaag                                                                                                                     <  1:253173/113‑7 (MQ=255)
                         tatGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGc                                                                                        >  1:1096434/1‑139 (MQ=255)
                              aTGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGgtg                                                                                   >  1:475668/1‑139 (MQ=255)
                                gCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCt                                                                                                        >  2:485295/1‑116 (MQ=255)
                                gCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCt                                                                                                        <  1:485295/116‑1 (MQ=255)
                                 cGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCt                                                                                                        >  1:958882/1‑115 (MQ=255)
                                 cGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCt                                                                                                        <  2:958882/115‑1 (MQ=255)
                                 cGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCg                                                                                >  1:1053294/1‑139 (MQ=255)
                                   ttAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaaag                                                                                                                     >  2:504944/1‑96 (MQ=255)
                                   ttAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaaag                                                                                                                     <  1:504944/102‑7 (MQ=255)
                                   ttAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTgg                                                                                                          <  2:24841/111‑1 (MQ=255)
                                   ttAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTgg                                                                                                          >  1:24841/1‑111 (MQ=255)
                                                  agTAGCGATGAAGGAAGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAAc                                                               >  2:756663/1‑139 (MQ=255)
                                                     aGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGt                                                            >  1:833248/1‑139 (MQ=255)
                                                       cGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatat                                                                                                                           <  1:497936/76‑1 (MQ=255)
                                                       cGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatat                                                                                                                           >  2:497936/1‑76 (MQ=255)
                                                          tGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatat                                                                                                                           <  1:102222/73‑1 (MQ=255)
                                                          tGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatat                                                                                                                           >  2:102222/1‑73 (MQ=255)
                                                           gaagaaCGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGAtt                                                      <  2:616696/139‑1 (MQ=255)
                                                              gaaCGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAAc                                                   >  1:365916/1‑139 (MQ=255)
                                                                 cGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGa                                                                               <  1:878936/108‑1 (MQ=255)
                                                                 cGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGa                                                                               >  2:937844/1‑108 (MQ=255)
                                                                 cGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGa                                                                               >  2:878936/1‑108 (MQ=255)
                                                                 cGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGa                                                                               <  1:937844/108‑1 (MQ=255)
                                                                    tGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatat                                                                                                                           >  1:126142/1‑63 (MQ=255)
                                                                    tGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatat                                                                                                                           <  2:126142/63‑1 (MQ=255)
                                                                                 tCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTa                                <  2:833248/139‑1 (MQ=255)
                                                                                     tGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaa                                                                                                                       <  1:780170/50‑5 (MQ=255)
                                                                                     tGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGatatattcaa                                                                                                                       >  2:780170/1‑46 (MQ=255)
                                                                                     tGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATcaca                                   >  2:289431/1‑132 (MQ=255)
                                                                                     tGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATcaca                                   <  1:289431/132‑1 (MQ=255)
                                                                                      gCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAagag                           >  1:1005325/1‑139 (MQ=255)
                                                                                        ccATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAATAATATCACAGTAAAGAGAt                         <  1:756663/139‑1 (MQ=255)
                                                                                         cATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAAt                                                                                                >  1:521225/1‑67 (MQ=255)
                                                                                         cATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAAt                                                                                                <  2:521225/67‑1 (MQ=255)
                                                                                          aTGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATaa                       >  2:13930/1‑139 (MQ=255)
                                                                                            gCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAata                                        <  1:791382/120‑1 (MQ=255)
                                                                                            gCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAata                                        >  2:791382/1‑120 (MQ=255)
                                                                                            gCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAAcc                     >  1:116275/1‑139 (MQ=255)
                                                                                                            tCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAAc                                                                          >  1:1109693/1‑70 (MQ=255)
                                                                                                            tCCCCGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAAc                                                                          <  2:1109693/70‑1 (MQ=255)
                                                                                                               ccGGTGTGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAACCCGAGAACTGAACGATGAAg  <  2:1096434/139‑1 (MQ=255)
                                                                                                                  gtgtGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAACCCGAGAACTGAACGAt      <  1:404300/132‑1 (MQ=255)
                                                                                                                  gtgtGCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAACCCGAGAACTGAACGAt      >  2:404300/1‑132 (MQ=255)
                                                                                                                    gtgCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAACCCGAGAACTGAACGATg     <  2:439477/131‑1 (MQ=255)
                                                                                                                    gtgCCAGCG‑‑ATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAACCCGAGAACTGAACGATg     >  1:439477/1‑131 (MQ=255)
                                                                                                                             ||                                                                                                                             
GTCCGTATGAAATAAATGTGACCTATATGGATGCGTTAAGCCGCCGTGAGAGTAGCGATGAAGAACGTTGCGCCAGGTTTATCCTTGCCCATGCGATTTTGTTAAGTTTCCCCGGTGTGCCAGCGATATATATTCAAAGTATTCTTGGCTCGCGTAATGATTACGCAGGTGTCGAAAAACTCGGATATAACCGTGCGATTAACCGTAAAAAATATCACAGTAAAGAGATAACCCGAGAACTGAACGATGAAG  >  NC_000913/1371351‑1371602

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: