Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,224,841 C→T H386Y (CAC→TAC)  pgi → glucosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,224,8410CT90.0% 165.0 / 10.0 61H386Y (CAC→TAC) pgiglucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (33/21);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

TATGCACCGTTTCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTAACTTCTTCGCCCAGACCGAAGC  >  CP009273/4224708‑4224974
                                                                                                                                     |                                                                                                                                     
tatGCACCGTTTCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGAGTACCAGACTGGCCCGATTATCTGGGGTGACCCAGGCACTAACGGTCAGTAcgcg                                                                                                                                  <  1:3718600/139‑1 (MQ=255)
    aaCCGATTCGCGGCGTACTTCCAGCAGGGCAATATGGAGCCCAACGGCACGTATGCTGACCGTGACGGTAACGTTGTGGAGTACCAGACTGGCCCGATTATCTGGGGAGAACCAGGCACTAACGGTCAGACCGCGTTCt                                                                                                                              >  1:3523087/2‑139 (MQ=255)
       cGTTTCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTAcc                                                                                                                           <  1:3202284/139‑1 (MQ=255)
        gTTTCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCa                                                                                                                          >  1:1030491/1‑139 (MQ=255)
        gTTTCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCa                                                                                                                          >  1:1342517/1‑139 (MQ=255)
        gTTTCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGCTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTcctgt                                                                                                                           >  1:4260610/1‑134 (MQ=255)
         tttCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTAcaa                                                                                                                          >  1:1283798/1‑136 (MQ=255)
         tttCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCa                                                                                                                          >  1:3537387/1‑138 (MQ=255)
            cgcgGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTg                                                                                                                      <  1:1351884/139‑1 (MQ=255)
                           gcagGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAAcc                                                                                                       >  1:2696995/1‑139 (MQ=255)
                                          gTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGt                                                                                                                                 >  1:1982685/1‑98 (MQ=255)
                                          gTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGc                                                                                        >  1:1412793/1‑139 (MQ=255)
                                              aaCGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGAtt                                                                                    >  1:1615744/1‑139 (MQ=255)
                                              aaCGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGAtt                                                                                    >  1:2407847/1‑139 (MQ=255)
                                                  gTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCa                                                                                 >  1:1710745/1‑138 (MQ=255)
                                                   tAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTg                                                                                         <  1:3257984/129‑1 (MQ=255)
                                                     aGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGACCCAAAATGGGACCGTGCGATTTCCTc                                                                               >  1:2364627/1‑137 (MQ=255)
                                                     aGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATc                                                                               >  1:491288/1‑137 (MQ=255)
                                                        aTGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCt                                                                            <  1:1366219/137‑1 (MQ=255)
                                                        aTGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCt                                                                            <  1:948417/137‑1 (MQ=255)
                                                           ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTAc                                                                                                                            >  1:751278/1‑86 (MQ=255)
                                                           ttGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCg                                                                         >  1:192664/1‑137 (MQ=255)
                                                            tGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGGCAGTACGCGt                                                                                                                                 >  1:3015145/1‑80 (MQ=255)
                                                             gACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATc                                                                               >  1:1086277/1‑129 (MQ=255)
                                                              aCCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTa                                                                     >  1:3222822/1‑138 (MQ=255)
                                                               ccGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTAt                                                                    <  1:3989169/138‑1 (MQ=255)
                                                               ccGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTAt                                                                    <  1:3455110/138‑1 (MQ=255)
                                                               ccGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTAt                                                                    <  1:1051068/138‑1 (MQ=255)
                                                                 gtaacggtaacgTTGTGGATTACCAGACTGGCCCGATTTCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCAc                                                                 >  1:1932385/1‑138 (MQ=255)
                                                                      ggtaacgTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCAt                                                             <  1:2763288/138‑1 (MQ=255)
                                                                         aacgTTGTGGATTCCCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGCTCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAAcc                                                         <  1:3786411/139‑1 (MQ=255)
                                                                         aacgTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAAcc                                                         <  1:1651390/139‑1 (MQ=255)
                                                                             ttGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGc                                                      >  1:4039290/1‑138 (MQ=255)
                                                                                tgGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAAc                                                          <  1:1749016/131‑1 (MQ=255)
                                                                                tgGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGctct                                                   >  1:590824/1‑138 (MQ=255)
                                                                                tgGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGctc                                                    >  1:3821982/1‑137 (MQ=255)
                                                                                  gATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATc                                                                                                                   <  1:3263907/72‑1 (MQ=255)
                                                                                      aCCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCaaa                                                                                                    <  1:387360/83‑1 (MQ=255)
                                                                                      aCCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGatca                                            <  1:782133/139‑1 (MQ=255)
                                                                                       ccAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGatcat                                           <  1:952781/139‑1 (MQ=255)
                                                                                          gACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGTTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTAATCGCTCCGGCTATCACCCATAACCCGCTCTCTg                                                >  1:3919363/1‑131 (MQ=255)
                                                                                           aCTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTAt                                                                    <  1:1869126/110‑1 (MQ=255)
                                                                                            cTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGaa                                                                                                         >  1:1092212/1‑72 (MQ=255)
                                                                                                cccGATTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGt                                                                                          >  1:2560091/1‑83 (MQ=255)
                                                                                                    aTTATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAAct                                >  1:427246/1‑137 (MQ=255)
                                                                                                     ttATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAAc                                                          >  1:4039265/1‑110 (MQ=255)
                                                                                                     ttATCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAAct                                >  1:4026121/1‑136 (MQ=255)
                                                                                                       aTCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCt                                                                            <  1:676371/90‑1 (MQ=255)
                                                                                                        tCTGGGGTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGatcat                                           <  1:2344998/122‑1 (MQ=255)
                                                                                                              gTGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATAATCACCAGAAACTGCTGTCTAActt                    >  1:146480/1‑139 (MQ=255)
                                                                                                              gTGAACCAGGCACTAACGGTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGAGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTaaat                     >  1:4305002/1‑136 (MQ=255)
                                                                                                               tGAACCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTAAct                     >  1:2425862/1‑137 (MQ=255)
                                                                                                                gAACCAGGCACTAACGGTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAACATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTCActtc                   <  1:3552529/138‑1 (MQ=255)
                                                                                                                 aaCCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCa                                                                  >  1:2848685/1‑90 (MQ=255)
                                                                                                                 aaCCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTAActtct                  >  1:2278600/1‑138 (MQ=255)
                                                                                                                  aCCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAAc                                 >  1:171920/1‑122 (MQ=255)
                                                                                                                  aCCAGGCACTAACGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTAACttcttc                >  1:3585894/1‑139 (MQ=255)
                                                                                                                             aCGGTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTAACTTCTTCGCCCAGACCGa     >  1:1255449/1‑139 (MQ=255)
                                                                                                                               ggTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAAc                                                          <  1:2020247/84‑1 (MQ=255)
                                                                                                                                gTCAGTACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCg                                                                                           <  1:59948/50‑1 (MQ=255)
                                                                                                                                gTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGCACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTAACTTCTTCGCCCAGACCGAAgc  >  1:819173/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TATGCACCGTTTCGCGGCGTACTTCCAGCAGGGCAATATGGAGTCCAACGGTAAGTATGTTGACCGTAACGGTAACGTTGTGGATTACCAGACTGGCCCGATTATCTGGGGTGAACCAGGCACTAACGGTCAGCACGCGTTCTACCAGCTGATCCACCAGGGAACCAAAATGGTACCGTGCGATTTCATCGCTCCGGCTATCACCCATAACCCGCTCTCTGATCATCACCAGAAACTGCTGTCTAACTTCTTCGCCCAGACCGAAGC  >  CP009273/4224708‑4224974

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