Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,554,731 C→T 100% M313I (ATG→ATA maeA ← malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,554,7310CT93.1% 220.7 / 5.0 76M313I (ATG→ATAmaeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (43/26);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGT  >  NC_000913/1554635‑1554827
                                                                                                |                                                                                                
cATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGattatt                                                                                               >  1:1864304/1‑100 (MQ=255)
cATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGattatt                                                                                               <  2:1743397/100‑1 (MQ=255)
cATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGATTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGattatt                                                                                               <  2:1743455/100‑1 (MQ=255)
 aTCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTAttt                                                                                              >  1:1115534/1‑100 (MQ=255)
 aTCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTAttt                                                                                              >  1:298494/1‑100 (MQ=255)
  tCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTc                                                                                             >  1:2046880/1‑100 (MQ=255)
  tCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTc                                                                                             <  1:7966/100‑1 (MQ=255)
   cTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGCTTATTTc                                                                                             >  2:2311584/1‑99 (MQ=255)
   cTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCg                                                                                            >  1:1349241/1‑100 (MQ=255)
             agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATCTCGGCAATGCCGc                                                                                  >  1:630333/1‑100 (MQ=255)
             agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGc                                                                                  >  1:65402/1‑100 (MQ=255)
              gCAAGCCAAAGCGGTCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCa                                                                                 >  2:1309511/1‑100 (MQ=255)
              gCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCa                                                                                 >  2:1024150/1‑100 (MQ=255)
              gCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCa                                                                                 >  2:599103/1‑100 (MQ=255)
              gCAAGCCAAAGCGATCGACCATAAAGACTTGCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCa                                                                                 <  1:1915869/100‑1 (MQ=255)
                 aGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATcc                                                                              >  2:1443003/1‑100 (MQ=255)
                 aGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATcc                                                                              >  1:1822200/1‑100 (MQ=255)
                  gCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATccc                                                                             >  1:890020/1‑100 (MQ=255)
                     aaaGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCt                                                                          >  1:393885/1‑100 (MQ=255)
                      aaGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTg                                                                         >  1:2040881/1‑100 (MQ=255)
                        gCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGaa                                                                       <  2:1864290/100‑1 (MQ=255)
                         ngATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAAc                                                                      >  1:452821/2‑100 (MQ=255)
                         cGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAAc                                                                      <  2:867586/100‑1 (MQ=255)
                         cGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAAc                                                                      >  1:453998/1‑100 (MQ=255)
                          gATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGCTGATTATTTCGGCAATGCCGCATCCCGCTGAAcc                                                                     >  2:2125776/1‑100 (MQ=255)
                             cGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTgc                                                                  >  2:237344/1‑100 (MQ=255)
                              gACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTgcg                                                                 <  1:2268416/100‑1 (MQ=255)
                                ccATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGcc                                                               <  2:55987/100‑1 (MQ=255)
                                 cATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCa                                                              <  2:1338171/100‑1 (MQ=255)
                                  aTAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCaa                                                             >  1:48811/1‑100 (MQ=255)
                                   tAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAg                                                            >  1:1532625/1‑100 (MQ=255)
                                      aGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGaa                                                         >  1:2023691/1‑100 (MQ=255)
                                      aGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGaa                                                         <  2:733393/100‑1 (MQ=255)
                                       gACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAg                                                        <  2:1729059/100‑1 (MQ=255)
                                        aCTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGa                                                       >  2:1617965/1‑100 (MQ=255)
                                         cTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAc                                                      >  2:1939075/1‑100 (MQ=255)
                                         cTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAc                                                      <  1:793298/100‑1 (MQ=255)
                                         cTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAc                                                      >  1:1089641/1‑100 (MQ=255)
                                         cTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAc                                                      >  2:779300/1‑100 (MQ=255)
                                              tGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttt                                                 >  1:619238/1‑100 (MQ=255)
                                               gccgagccgcTTCCTCGCTTAATCCTTCGCGCTGGGACTGGGAGCTGCTTATTTCGGCCATGCCGCACCCCGCGGCATCTGCTCCCAGGaaacc                                                      >  2:1413627/1‑91 (MQ=255)
                                                ccgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttttt                                               <  1:47633/100‑1 (MQ=255)
                                                 cgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAtttttt                                              >  1:1203909/1‑100 (MQ=255)
                                                   cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                            <  1:1562877/100‑1 (MQ=255)
                                                   cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                            <  2:403426/100‑1 (MQ=255)
                                                     ccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtctc                                          <  1:1945254/100‑1 (MQ=255)
                                                        cTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCt                                       <  1:154262/100‑1 (MQ=255)
                                                          tCCTCTCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTa                                     <  2:601920/100‑1 (MQ=255)
                                                           ccTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAAAATTTTTTTCCCGCTTa                                     >  1:438281/1‑99 (MQ=255)
                                                            cTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAAc                                   >  2:618528/1‑100 (MQ=255)
                                                            cTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAAc                                   >  1:1199708/1‑100 (MQ=255)
                                                              cGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGGCGATTTTTTTCTCGCTTAACTg                                 <  2:1381453/100‑1 (MQ=255)
                                                              cGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTg                                 >  1:1289943/1‑100 (MQ=255)
                                                              cGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTg                                 >  1:1928175/1‑100 (MQ=255)
                                                              cGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTg                                 <  2:1699243/100‑1 (MQ=255)
                                                              cGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTg                                 >  1:1462429/1‑100 (MQ=255)
                                                                  tAATCCTTCGCGCTGGGTCTTGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGacc                              <  2:1589082/99‑1 (MQ=255)
                                                                  tAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGacca                             >  1:731200/1‑100 (MQ=255)
                                                                      ccTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGaccacc                           <  2:373644/98‑1 (MQ=255)
                                                                       cTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                        >  1:64550/1‑100 (MQ=255)
                                                                       cTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                        <  1:618530/100‑1 (MQ=255)
                                                                       cTTCGCGCTGGGTCTGGGAGATGATAACTTCGGCAATGCCGCATCCGACTGAACCTGCACCTCGGAAGACGGCTTTTTTCTCGcct                                      >  2:721849/1‑84 (MQ=255)
                                                                        ttCGCTCTGGAACTAGGAGACGACAATTTCGGCAATGCCGCATCCCGCTGAACCTGCCCCAAGAAACACGATTTTTTTCTCGCTCAACTGACCACCTGcc                       >  2:244300/1‑100 (MQ=255)
                                                                            cgcTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcg                      >  2:950344/1‑97 (MQ=255)
                                                                              cTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCgg                 <  1:2191202/100‑1 (MQ=255)
                                                                              cTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCgg                 <  2:65401/100‑1 (MQ=255)
                                                                               tGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGc                <  1:628637/100‑1 (MQ=255)
                                                                                  gTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtg             <  1:1482086/100‑1 (MQ=255)
                                                                                  gTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtg             >  2:756357/1‑100 (MQ=255)
                                                                                   tCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATATTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgc            <  2:1463134/100‑1 (MQ=255)
                                                                                    cTGGGAGCTGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTCTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgct           <  2:1275839/100‑1 (MQ=255)
                                                                                    cTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgct           <  2:708268/100‑1 (MQ=255)
                                                                                      gggAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgctgc         >  2:152091/1‑100 (MQ=255)
                                                                                        gagaTGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGa       >  2:1860382/1‑100 (MQ=255)
                                                                                         agaTGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGAt      <  2:1893603/100‑1 (MQ=255)
                                                                                            tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATc     >  2:1653471/1‑98 (MQ=255)
                                                                                            tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAg   >  2:1777873/1‑100 (MQ=255)
                                                                                            tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAg   >  1:1229358/1‑100 (MQ=255)
                                                                                             gatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAgt  >  1:1877735/1‑100 (MQ=255)
                                                                                                |                                                                                                
CATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGT  >  NC_000913/1554635‑1554827

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: