Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,214,964 A→G I229V (ATT→GTT)  metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,214,9640AG95.8% 235.8 / ‑1.9 72I229V (ATT→GTT) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (35/34);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCG  >  NC_000913/4214867‑4215059
                                                                                                 |                                                                                               
cTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGtt                                                                                               <  1:578785/100‑1 (MQ=255)
    ccGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCt                                                                                           >  1:1295219/1‑100 (MQ=255)
     cGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCtt                                                                                          <  2:888732/100‑1 (MQ=255)
     cGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCtt                                                                                          >  1:999531/1‑100 (MQ=255)
      ggCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGAGCGTTgaatt                                                                                          >  2:1387334/1‑95 (MQ=255)
       gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtg                                                                                        >  2:293633/1‑100 (MQ=255)
       gCAGAGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtg                                                                                        <  2:718531/100‑1 (MQ=255)
        cAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtgt                                                                                       >  2:250918/1‑100 (MQ=255)
         aGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtgtg                                                                                      <  1:246167/100‑1 (MQ=255)
               gATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGc                                                                                >  1:131707/1‑100 (MQ=255)
               gATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGc                                                                                <  1:1777340/100‑1 (MQ=255)
                     tGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATccc                                                                          >  2:889898/1‑100 (MQ=255)
                     tGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATccc                                                                          <  2:588885/100‑1 (MQ=255)
                       aTTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGa                                                                        <  2:791444/100‑1 (MQ=255)
                        ttACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGATGGGCCATCCCGa                                                                        <  1:371049/99‑1 (MQ=255)
                         tACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAat                                                                      <  1:1065140/100‑1 (MQ=255)
                         tACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAat                                                                      <  2:1509519/100‑1 (MQ=255)
                          acacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAata                                                                     <  2:753256/100‑1 (MQ=255)
                           cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAatat                                                                    >  1:72386/1‑100 (MQ=255)
                             ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatga                                                                  >  2:463545/1‑100 (MQ=255)
                              cGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatgat                                                                 >  2:1175802/1‑100 (MQ=255)
                              cGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatgat                                                                 >  2:748215/1‑100 (MQ=255)
                               gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatgatg                                                                <  2:1162364/100‑1 (MQ=255)
                               gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatgatg                                                                <  1:1499485/100‑1 (MQ=255)
                                aTCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATgc                                                               >  1:1067609/1‑100 (MQ=255)
                                  cTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATgcgc                                                             >  2:1444488/1‑100 (MQ=255)
                                  cTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATgcgc                                                             <  2:131707/100‑1 (MQ=255)
                                   tGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCa                                                            <  2:29586/100‑1 (MQ=255)
                                   tGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCa                                                            <  2:1295230/100‑1 (MQ=255)
                                    ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaa                                                           >  1:211122/1‑100 (MQ=255)
                                      aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                         >  1:793839/1‑100 (MQ=255)
                                      aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                         >  2:1361150/1‑100 (MQ=255)
                                            tGGCAGAGACGGAAGAAGTGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGgc                                                   <  1:1387324/100‑1 (MQ=255)
                                                agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                               >  2:1185484/1‑100 (MQ=255)
                                                   gaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGACTATGATGCGCAAACGCTGGCGCAGGaa                                            >  2:722756/1‑100 (MQ=255)
                                                        aagaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAttttt                                       <  2:1654391/100‑1 (MQ=255)
                                                         agaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTc                                      <  1:169629/100‑1 (MQ=255)
                                                           aagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcg                                    >  1:10101/1‑100 (MQ=255)
                                                            agGGTATGCATATCTGTTTGCCAGTACAGATATGCACATTTCCTTTGTGGCGGTCCATCCCGACTATGATACGCAAGAGCTGGCCAAGGAATTTATcct                                    <  1:484021/99‑2 (MQ=255)
                                                            agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcgc                                   >  2:1337629/1‑100 (MQ=255)
                                                            agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGGTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcgc                                   <  1:1570658/100‑1 (MQ=255)
                                                               ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                <  1:384643/100‑1 (MQ=255)
                                                                 aTGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtgt                              <  1:492749/100‑1 (MQ=255)
                                                                    cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                           >  2:245946/1‑100 (MQ=255)
                                                                    cATATCTGTTTGCCAGTAAAGATAAACGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCAATGTGGa                           >  1:605966/1‑100 (MQ=255)
                                                                      tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                         <  2:922342/100‑1 (MQ=255)
                                                                      tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                         >  2:1413796/1‑100 (MQ=255)
                                                                      tatCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                         >  1:145615/1‑100 (MQ=255)
                                                                        tCTGTTTGCCAGTAAAGATAAGCGCGTTTCCTTTGTTACTGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGcc                       <  2:210735/100‑1 (MQ=255)
                                                                          tGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCgg                     <  2:1067619/100‑1 (MQ=255)
                                                                          tGTTTGCCAGTAAAGATAAGAGCGTTGCCTTTGTGACGGACCATGCCGACAATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGAAAGCAgg                     >  1:151033/1‑100 (MQ=255)
                                                                           gTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                    <  2:1554955/100‑1 (MQ=255)
                                                                           gTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                    >  1:308790/1‑100 (MQ=255)
                                                                            tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                   <  2:211123/100‑1 (MQ=255)
                                                                            cttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                   <  1:405825/99‑1 (MQ=255)
                                                                              tGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAcag                 >  2:405827/1‑98 (MQ=255)
                                                                              tGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTa                 <  2:879624/100‑1 (MQ=255)
                                                                                 cAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAc              >  2:99025/1‑100 (MQ=255)
                                                                                  aGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAcc             >  2:257432/1‑100 (MQ=255)
                                                                                  aGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAc              >  2:1083634/1‑99 (MQ=255)
                                                                                   gTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAccc            >  1:1636809/1‑100 (MQ=255)
                                                                                    tAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTTTGTTTTTTGGAAGCCGGAGTAGAccc            >  1:359500/1‑99 (MQ=255)
                                                                                      aaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa         <  2:1503229/100‑1 (MQ=255)
                                                                                       aGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt        >  1:1229519/1‑100 (MQ=255)
                                                                                       aGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt        <  1:77742/100‑1 (MQ=255)
                                                                                        gATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATg       >  1:1599619/1‑100 (MQ=255)
                                                                                        gATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATg       >  2:1706063/1‑100 (MQ=255)
                                                                                        gATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGACGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATg       <  2:793844/100‑1 (MQ=255)
                                                                                         aTAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGt      <  1:1359932/100‑1 (MQ=255)
                                                                                          tAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTa     >  1:158360/1‑100 (MQ=255)
                                                                                            aGCGCGTTGCCTTTGTGACGGGCCATCCCGAATACGATCCGCAACCGCTGGCGCAGGAATTTTTCCCCAATGTAGAAGCCGGACTAGACCCGGATGTAcc   <  1:81564/100‑1 (MQ=255)
                                                                                             gcgcGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCg  <  2:1048710/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCG  >  NC_000913/4214867‑4215059

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