Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 638,614 G→T T14K (ACA→AAA)  thiD ← phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730638,6140GT96.4% 184.4 / ‑2.4 55T14K (ACA→AAA) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (29/24);  total (30/25)
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

CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAAT  >  CP000730/638477‑638750
                                                                                                                                         |                                                                                                                                        
cATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTttgt                                                                                                                                       <  1:112012/141‑1 (MQ=255)
   ccATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAg                                                                                                                                    >  2:89507/1‑141 (MQ=255)
    cATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGa                                                                                                                                   >  1:558094/1‑141 (MQ=255)
     aTTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGaa                                                                                                                                  <  2:652792/141‑1 (MQ=255)
      ttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAc                                                                                                                                 <  2:423917/141‑1 (MQ=255)
      ttGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAc                                                                                                                                 >  1:285111/1‑141 (MQ=255)
           aaTGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCa                                                                                                                            >  2:647860/1‑141 (MQ=255)
             tGGTGTAACATCGTGCGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGAGCATTCCATACGAATCTAATACATGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAAt                                                                                                                          >  1:597741/1‑141 (MQ=255)
              ggTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAAtt                                                                                                                         >  1:354087/1‑141 (MQ=255)
                 gtAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGtt                                                                                                                      <  2:614366/141‑1 (MQ=255)
                  tAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTa                                                                                                                     >  2:272849/1‑141 (MQ=255)
                       tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAct                                                                                                                >  2:92663/1‑141 (MQ=255)
                       tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAct                                                                                                                <  2:482418/141‑1 (MQ=255)
                       tCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAAct                                                                                                                >  1:208791/1‑141 (MQ=255)
                             gcccAGGTATCTTTATCCATAGTAATGTTGCCAGTTATGGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCATTTGTTAAAActttcttt                                                                                                          <  2:345903/139‑1 (MQ=255)
                                   gTATCTGTATCCATAGTAACGATGTCAGTGAAATCCGCCATGCCAGACCAAGCAAAGTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGcc                                                                                                    <  1:592546/141‑1 (MQ=255)
                                         ttATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAg                                                                                              <  1:280736/141‑1 (MQ=255)
                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                                                                             <  1:356498/141‑1 (MQ=255)
                                          tATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGc                                                                                             <  2:53308/141‑1 (MQ=255)
                                               aTAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                        >  2:462689/1‑141 (MQ=255)
                                               aTAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCAc                                                                                        >  1:537018/1‑141 (MQ=255)
                                                    aaCGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTa                                                                                   >  1:382898/1‑141 (MQ=255)
                                                          ggCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACagatt                                                                             <  2:262772/141‑5 (MQ=255)
                                                             aGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCCTTAAAGCCATTGAGCTTCACTCCTACataataat                                                                          >  1:609292/1‑141 (MQ=255)
                                                                tAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGATCTTCACTCCTACATAATAATAtt                                                                       <  1:165340/141‑1 (MQ=255)
                                                                    gCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTAt                                                                   >  2:486571/1‑141 (MQ=255)
                                                                      gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                >  1:697297/1‑78 (MQ=255)
                                                                      gACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAAcc                                                                                                                                <  2:697297/78‑1 (MQ=255)
                                                                          aTGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATtcatca                                                             <  2:597741/141‑1 (MQ=255)
                                                                            gCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcata                                                           >  1:335288/1‑141 (MQ=255)
                                                                             ccATACGTATCTGATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCatcatat                                                          <  2:279404/141‑1 (MQ=255)
                                                                                 aCGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAtt                                                      >  1:340481/1‑141 (MQ=255)
                                                                                  cGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttt                                                     <  1:255115/141‑1 (MQ=255)
                                                                                  cGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCAttt                                                     <  1:92663/141‑1 (MQ=255)
                                                                                       cTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAAc                                                <  1:24919/141‑1 (MQ=255)
                                                                                          aTTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTa                                             >  2:654054/1‑141 (MQ=255)
                                                                                                ggAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGaa                                       >  2:188784/1‑141 (MQ=255)
                                                                                                     gTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAATCTAATTGAAAAata                                  >  1:567415/1‑141 (MQ=255)
                                                                                                         tCAAATCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTaa                              >  1:196682/1‑141 (MQ=255)
                                                                                                             aTCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCAt                          <  1:340228/141‑1 (MQ=255)
                                                                                                             aTCTGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCAt                          >  2:340228/1‑141 (MQ=255)
                                                                                                                tGCTTGCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCa                       >  1:596537/1‑141 (MQ=255)
                                                                                                                     gCATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATAtt                  >  1:322936/1‑141 (MQ=255)
                                                                                                                      cATACCTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCCTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATAttt                 <  2:609292/141‑1 (MQ=255)
                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAg                                                                                                      >  2:509305/1‑51 (MQ=255)
                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAg                                                                                                      <  1:509305/51‑1 (MQ=255)
                                                                                                                           cTGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTtgatga            <  1:299967/141‑1 (MQ=255)
                                                                                                                            tGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat           <  2:339052/141‑1 (MQ=255)
                                                                                                                            tGCGCCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTgatgat           >  2:241149/1‑141 (MQ=255)
                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt        <  2:326994/141‑1 (MQ=255)
                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt        <  1:278609/141‑1 (MQ=255)
                                                                                                                               gcCAGCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGAttgt        >  2:8259/1‑141 (MQ=255)
                                                                                                                                   gCACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaa    >  2:313208/1‑141 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATGGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa   >  2:248324/1‑141 (MQ=255)
                                                                                                                                    cACTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCAGTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGaaa   <  2:235660/141‑1 (MQ=255)
                                                                                                                                     aCTTTTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTCAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAAt  >  1:520307/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CATCCATTGGTAATGGTGTAACATCGTGTGACCATGTATCTTTATCCATAGTAACGATGGCAGTTAAAGCGACCATGCCATACGTATCTAATTCTTGGAACGTTTTCAAATCTGCTTGCATACCTGCGCCAGCACTTGTGTCAGAACCGGCAATTGTTAAAACTTTCTTTAAAGCCATTGAGCTTCACTCCTACATAATAATATTGTATTCATCATATCATTTTTAACCTAATTGAAAAATATTAAGCATTCAATATTTGATGATTGTTGAAAT  >  CP000730/638477‑638750

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