Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,718,745 T→G T121P (ACG→CCG)  kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,718,7450TG98.3% 204.9 / ‑3.5 58T121P (ACG→CCG) kgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (29/28);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGG  >  NZ_CP009273/2718613‑2718878
                                                                                                                                    |                                                                                                                                     
gAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCta                                                                                                                                  >  1:1265964/1‑138 (MQ=255)
 aaaTGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCtat                                                                                                                                 <  1:867648/138‑1 (MQ=255)
 aaaTGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCtat                                                                                                                                 <  1:586646/138‑1 (MQ=255)
 aaaTGATGCGTAAAAACCTTTGCGCCCTTCAACGGAAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCtata                                                                                                                                <  2:134397/139‑1 (MQ=255)
       tgCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCa                                                                                                                          <  2:55339/139‑1 (MQ=255)
       tgCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCa                                                                                                                          <  2:250196/139‑1 (MQ=255)
         cGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCAt                                                                                                                         <  1:891882/138‑1 (MQ=255)
         cGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAACCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCAt                                                                                                                         >  2:891882/1‑138 (MQ=255)
           tAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAAc                                                                                                                      >  1:1268373/1‑139 (MQ=255)
                    ttGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCt                                                                                                                                   >  2:315785/1‑117 (MQ=255)
                    ttGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCt                                                                                                                                   <  1:315785/117‑1 (MQ=255)
                            ttCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGAt                                                                                                      >  2:227682/1‑138 (MQ=255)
                                 cGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                                                                                <  1:873704/139‑1 (MQ=255)
                                   gCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCg                                                                                              <  2:600374/139‑1 (MQ=255)
                                      aCTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATc                                                                                           >  1:852424/1‑139 (MQ=255)
                                          cACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGa                                                                                        <  1:1160668/138‑1 (MQ=255)
                                           aCTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGaa                                                                                       <  2:1268373/138‑1 (MQ=255)
                                           aCTCATATAGGTGGCGCTGGTGACATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGaa                                                                                       <  1:358161/138‑1 (MQ=255)
                                                 ataGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                                                                                <  2:840806/123‑1 (MQ=255)
                                                 ataGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                                                                                >  1:840806/1‑123 (MQ=255)
                                                  taGGTGGCGCTGGGGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACAcatc                                                                                <  2:259369/138‑1 (MQ=255)
                                                     gTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATac                                                                            <  1:299942/139‑1 (MQ=255)
                                                       ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCa                                                                                                                          <  2:825641/91‑1 (MQ=255)
                                                       ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCa                                                                                                                          >  1:825641/1‑91 (MQ=255)
                                                       ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAac                                                                                     <  2:859013/128‑1 (MQ=255)
                                                       ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAac                                                                                     >  1:859013/1‑128 (MQ=255)
                                                       ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATaca                                                                           >  2:281742/1‑138 (MQ=255)
                                                       ggCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATaca                                                                           >  2:50334/1‑138 (MQ=255)
                                                        gcgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATacac                                                                          <  1:269360/138‑1 (MQ=255)
                                                         cgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCtata                                                                                                                                >  1:864187/1‑83 (MQ=255)
                                                         cgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCtata                                                                                                                                <  2:864187/83‑1 (MQ=255)
                                                         cgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGaa                                                                                       <  1:275275/124‑1 (MQ=255)
                                                         cgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGaa                                                                                       >  2:275275/1‑124 (MQ=255)
                                                         cgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATacacac                                                                        >  2:76823/1‑139 (MQ=255)
                                                         cgcTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATacacac                                                                        >  1:431593/1‑139 (MQ=255)
                                                             ggTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGa                                                                     >  1:346858/1‑138 (MQ=255)
                                                              gTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATac                                                                            >  1:824275/1‑130 (MQ=255)
                                                              gTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATac                                                                            <  2:824275/130‑1 (MQ=255)
                                                                          ccGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAtt                                                       >  1:69407/1‑139 (MQ=255)
                                                                           cgacAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCACAAGGACATAAAGTTTCATAACATGGGAGGCAGGCGATAACAAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttt                                                      >  2:967290/4‑139 (MQ=255)
                                                                           cGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCtata                                                                                                                                <  1:1063513/65‑1 (MQ=255)
                                                                           cGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCtata                                                                                                                                >  2:1063513/1‑65 (MQ=255)
                                                                           cGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttt                                                      >  2:1166678/1‑139 (MQ=255)
                                                                           cGCCAACAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGAttt                                                      >  1:1187120/1‑139 (MQ=255)
                                                                                aCAGATAATCCCTGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGc                                                 >  2:151280/1‑139 (MQ=255)
                                                                                     tAATCCCTGAAATAAACGAGCGATGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACTCATCATACACTCCGATAACAGCATCGATTTTTTGCGTCCa                                            <  1:967290/139‑1 (MQ=255)
                                                                                            tGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                                                                                >  2:154057/1‑80 (MQ=255)
                                                                                            tGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAg                                                                                                <  1:154057/80‑1 (MQ=255)
                                                                                            tGAAATAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTAt                                     >  2:207176/1‑139 (MQ=255)
                                                                                               aaTAAACGAGCGAGGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCgg                                  <  1:281742/139‑1 (MQ=255)
                                                                                                        gcgagGAGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGc                         >  2:370866/1‑139 (MQ=255)
                                                                                                           aggaGCAATAATGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCa                       >  1:523773/4‑138 (MQ=255)
                                                                                                                  ataatGCCGGAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCa               <  2:69407/139‑1 (MQ=255)
                                                                                                                           gAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACAcatca                                                                               >  2:470274/1‑66 (MQ=255)
                                                                                                                           gAGCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACAcatca                                                                               <  1:470274/66‑1 (MQ=255)
                                                                                                                             gCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTg                                                  >  2:29103/1‑93 (MQ=255)
                                                                                                                             gCCCACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTg                                                  <  1:29103/93‑1 (MQ=255)
                                                                                                                               ccACGGACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTggg  <  1:1116191/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
GAAATGATGCGTAAAAACCTTTGCGCCCTTCAACGGCAACTTCACTCATATAGGTGGCGCTGGTGCCATATTCTCCGCCAACAGATAATCCCTGAAATAAACGAGCGAGAAGCAATAATGCCGGAGCCCACGTACCTATAGTTTCATAACCTGGGAGGCAGGCGATAACCAGCGATCCGAAACACATCATACACACCGATAACAGCATCGATTTTTTGCGACCATGTTTATCGGCTATGCGGCCAAATAGCCAACCGCCTATTGGG  >  NZ_CP009273/2718613‑2718878

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