Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,555,655 (T)7→8 coding (15/1698 nt) maeA ← malate dehydrogenase, NAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,6481.T94.1% 167.0 / 2.7 51coding (22/1698 nt)maeAmalate dehydrogenase, NAD‑requiring
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base T (24/24);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CCACTTCCGGCAGTAACCCCAGCAGGTTGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGC  >  NC_000913/1555522‑1555780
                                                                                                                                |                                                                                                                                      
ccACTTCCGGCAGTAACCCCAGCAGGTTGAAGTTAGGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTtttttgtttt                                                                                                                          <  2:301518/140‑1 (MQ=255)
ccACTTCCGGCAGTAACCCCAGCAGGTTGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTtttttgtttt                                                                                                                          <  2:734504/140‑1 (MQ=255)
ccACTTCCGGCAGTAACCCCAGCAGGTTGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTtttttgtttt                                                                                                                          <  1:160027/140‑1 (MQ=255)
           aGTAACCCCAGCAGGTTGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGt                                                                                                               <  2:59040/140‑1 (MQ=255)
              aaCCCCAGCAGGTTGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGtcac                                                                                                            >  2:32093/1‑140 (MQ=255)
                           tGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTg                                                                                               >  1:1883/1‑140 (MQ=255)
                            gAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGa                                                                                              <  1:499852/140‑1 (MQ=255)
                            gAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGa                                                                                              <  2:178261/140‑1 (MQ=255)
                                   cGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATatcca                                                                                       >  2:582022/1‑140 (MQ=255)
                                            tcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTgggg                                                                              <  1:652806/140‑1 (MQ=255)
                                                  cTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCttt                                                                        >  2:35951/1‑140 (MQ=255)
                                                   tGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCtttt                                                                       <  2:683681/140‑1 (MQ=255)
                                                    gAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTa                                                                      <  2:1883/140‑1 (MQ=255)
                                                            tGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCttt                                                              >  1:720908/1‑140 (MQ=255)
                                                                        acaacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCGTCCCTGGGGGGGTTGTATCGTCTTTGCTTTACCGCCa                                                  >  1:228213/1‑140 (MQ=255)
                                                                        acaacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGGTTATATCGTGTTTGCTTTACCGCCa                                                  >  1:299260/1‑140 (MQ=255)
                                                                        acaacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCa                                                  >  2:121816/1‑140 (MQ=255)
                                                                        acaacGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCa                                                  >  2:52407/1‑140 (MQ=255)
                                                                            cGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGGATATCCATTCCTGGGGGCTTTTTATAGTCTAAGCTTTACGGCCAGGGc                                              >  2:355932/1‑140 (MQ=255)
                                                                             gacagttCCAGCAGTACCGGGCCAGCGTAAGGGATTTAATGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCg                                             <  2:301324/135‑1 (MQ=255)
                                                                                  ttCCAGCAGTACAGGTCCAGCGTAACGGATATAAAGCTAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCCCTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        <  1:310158/140‑1 (MQ=255)
                                                                                  ttCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        <  2:452549/140‑1 (MQ=255)
                                                                                  ttCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        <  1:121816/140‑1 (MQ=255)
                                                                                  ttCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        <  2:579298/140‑1 (MQ=255)
                                                                                  ttCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        <  2:272317/140‑1 (MQ=255)
                                                                                     cagcagTACAGGGCCTGCGTAAGGGCTATAAAGCGAAGGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc                                     <  1:355932/140‑1 (MQ=255)
                                                                                     cagcagTACAGGGCCAGCGTAAGGGATATAACGCGCACGCTGTTTTTTCTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc                                     <  2:299260/140‑1 (MQ=255)
                                                                                     cagcagTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc                                     <  1:582022/140‑1 (MQ=255)
                                                                                        cagTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAg                                  >  2:125480/1‑140 (MQ=255)
                                                                                           tACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGtat                               >  1:387384/1‑140 (MQ=255)
                                                                                             cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACTGCCAGGGCGTCGGCCTCAAgctgtc                             >  1:301324/1‑135 (MQ=255)
                                                                                             cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        >  1:593075/1‑129 (MQ=255)
                                                                                             cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        >  1:573031/1‑129 (MQ=255)
                                                                                             cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        <  2:593075/129‑1 (MQ=255)
                                                                                             cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGc                                        <  2:573031/129‑1 (MQ=255)
                                                                                             cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATaa                             >  2:621684/1‑140 (MQ=255)
                                                                                                  ccAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc                                     <  1:561115/127‑1 (MQ=255)
                                                                                                  ccAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc                                     >  2:561115/1‑127 (MQ=255)
                                                                                                     gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc                                     >  2:728472/1‑124 (MQ=255)
                                                                                                     gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTc                                     <  1:728472/124‑1 (MQ=255)
                                                                                                     gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATa                     >  2:520326/1‑140 (MQ=255)
                                                                                                        tAAGGGATATAAAGCGAACGCTCTTTTTTTTGTTTCTGGTTCCATGTCGCTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaa                  <  2:196604/140‑1 (MQ=255)
                                                                                                            ggATATAAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCAAGGGCGTCGGCGTCTAGTATAAAGCAGATaaaaacaa              >  1:196604/1‑140 (MQ=255)
                                                                                                               tataAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaacaaaaa           <  2:27337/140‑1 (MQ=255)
                                                                                                                 taAAGCGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAaaaacaaaaaca         <  1:520326/140‑1 (MQ=255)
                                                                                                                     gCGAACGCTGTTTTTCTTGTTTATGTTTCCATGTCACCCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCAt     <  1:292907/140‑1 (MQ=255)
                                                                                                                      cGAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCAtt    >  1:173603/1‑140 (MQ=255)
                                                                                                                       gAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTg   >  1:151604/1‑140 (MQ=255)
                                                                                                                       gAACGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCAAGGGCGTCGGCCACAAGTATAAAGCCGATAAAAACAAAAACACCATTg   >  2:292907/1‑140 (MQ=255)
                                                                                                                        aaCGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc  <  1:52407/140‑1 (MQ=255)
                                                                                                                        aaCGCTGTTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc  >  1:168541/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                      
CCACTTCCGGCAGTAACCCCAGCAGGTTGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGC  >  NC_000913/1555522‑1555780

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