Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,004,879 A→G 7.5% L259S (TTG→TCG)  pdxB ← D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,004,8790AG 7.5% 151.9 / 3.7 53L259S (TTG→TCG) pdxBD‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base G (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CATTCTTATTGAGAAAATTAACTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCATCTTTCATT  >  CP000730/2004745‑2005007
                                                                                                                                      |                                                                                                                                
cATTCTTATTGAGAAAATTAACTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGtt                                                                                                                            >  2:375708/1‑141 (MQ=255)
  ttCTTATTGAGAAAATTAACTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCa                                                                                                                          <  2:278847/141‑1 (MQ=255)
      tATTGAGAAAATTAACTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCACTTCATATAATTCATGATTAGGTTTCAAAGGTTCAtttt                                                                                                                      >  1:84823/1‑141 (MQ=255)
         tGAGAAAATTAACTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCaa                                                                                                                   >  1:339785/1‑141 (MQ=255)
         tGAGAAAATTAACTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATCCATGATTAGGTTTCAAAGGTTCATTTTCaa                                                                                                                   <  2:706414/141‑1 (MQ=255)
                    aCTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTaaa                                                                                                        >  1:475864/1‑141 (MQ=255)
                     cTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAat                                                                                                       <  1:349644/141‑1 (MQ=255)
                      tAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAata                                                                                                      >  2:543439/1‑141 (MQ=255)
                           tGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCAt                                                                                                 <  2:634375/141‑1 (MQ=255)
                           tGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCAt                                                                                                 <  2:379395/141‑1 (MQ=255)
                                 tAAAAATATCTAATATGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGaa                                                                                           <  1:242232/141‑1 (MQ=255)
                                 tAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGaa                                                                                           >  1:200064/1‑141 (MQ=255)
                                   aaaaTATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATa                                                                                         <  1:688274/141‑1 (MQ=255)
                                     aaTATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAac                                                                                       <  1:682475/141‑1 (MQ=255)
                                       tttttaataaGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAactt                                                                                     >  2:75990/5‑141 (MQ=255)
                                       tatCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAactt                                                                                     >  1:66053/1‑141 (MQ=255)
                                        atCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAacttt                                                                                    <  2:187966/141‑1 (MQ=255)
                                         tCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCGAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAacttta                                                                                   >  1:13210/1‑141 (MQ=255)
                                         tCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCGAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAacttta                                                                                   >  1:13386/1‑141 (MQ=255)
                                                  tCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACCCCTAAATATGCATGTCGAATAACTTTACTTTTTAAt                                                                          <  2:626361/141‑1 (MQ=255)
                                                  tCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAAt                                                                          <  1:25467/141‑1 (MQ=255)
                                                   cATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATa                                                                         >  2:189670/1‑141 (MQ=255)
                                                     tACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACt                                                                       <  2:200064/141‑1 (MQ=255)
                                                             ttCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTaa                                                               >  2:72580/1‑141 (MQ=255)
                                                              tCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAag                                                              >  2:125142/1‑141 (MQ=255)
                                                               cATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGCCGAATAACTTTACTTTTTAATACTTCTATTAaga                                                             <  2:66371/141‑1 (MQ=255)
                                                                 tAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAgc                                                           >  2:660347/1‑141 (MQ=255)
                                                                  aaTCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAgcg                                                          <  1:567618/141‑1 (MQ=255)
                                                                    tCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGct                                                        <  2:535788/141‑1 (MQ=255)
                                                                      aTTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGcttc                                                      >  2:721608/1‑141 (MQ=255)
                                                                      aTTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTACGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGcttc                                                      >  2:313743/1‑141 (MQ=255)
                                                                        tACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGcttctt                                                    >  2:233292/1‑141 (MQ=255)
                                                                         aCCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCttt                                                   <  2:496318/141‑1 (MQ=255)
                                                                             gTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACt                                               >  1:540098/1‑141 (MQ=255)
                                                                                 tatGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGc                                           <  2:170191/141‑1 (MQ=255)
                                                                                  atGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCt                                          <  1:543439/141‑1 (MQ=255)
                                                                                    gcgcTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTAc                                        <  1:290115/141‑1 (MQ=255)
                                                                                      gcTGTNATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCt                                      <  1:261670/141‑1 (MQ=255)
                                                                                            ataGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCt                                >  1:165606/1‑141 (MQ=255)
                                                                                                ttACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATATAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGtaaaaaaa                            <  2:341090/141‑8 (MQ=255)
                                                                                                        tCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATaaaa                    >  2:392858/1‑141 (MQ=255)
                                                                                                         ccAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATaaaaa                   >  1:764207/1‑141 (MQ=255)
                                                                                                              tCATATANTTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCt              <  2:588564/141‑1 (MQ=255)
                                                                                                              tCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCt              <  1:59032/141‑1 (MQ=255)
                                                                                                                 tataATTCATGATTAGGTTTCGAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCa           <  2:13386/141‑1 (MQ=255)
                                                                                                                 tataATTCATGATTAGGTTTCGAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCa           <  2:13210/141‑1 (MQ=255)
                                                                                                                 tataATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCa           >  2:87948/1‑141 (MQ=255)
                                                                                                                  ataATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTTTaatcagtccatttt          >  2:228097/1‑129 (MQ=255)
                                                                                                                  ataATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCAt          >  2:467201/1‑141 (MQ=255)
                                                                                                                    aaTTCATGATTAGGTTTCAAAGGTTCATTTTCAAAGACATCTTAATATGCATGTCGAATAACTTTACATTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCATCt        >  1:676403/1‑141 (MQ=255)
                                                                                                                     aTTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCATCtt       <  2:17184/141‑1 (MQ=255)
                                                                                                                         aTGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCATCTTTCAt   >  2:496230/1‑141 (MQ=255)
                                                                                                                          tGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCATCTTTCAtt  >  1:145167/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                
CATTCTTATTGAGAAAATTAACTAGATTGTTTTTAAAAATATCTAATAAGTCATACTTTGCTTCATAATCATTACCAGTTATATGCGCTGTTATAGTTACATTATCCAATTCATATAATTCATGATTAGGTTTCAAAGGTTCATTTTCAAACACATCTAAATATGCATGTCGAATAACTTTACTTTTTAATACTTCTATTAAGAGCGCTTCTTTAACTATGCTACCTCGTCCTATATTTATAAAAAGTGCTTCATCTTTCATT  >  CP000730/2004745‑2005007

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