Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,820,976 C→T A69T (GCG→ACG)  astE ← succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,820,9760CT92.6% 176.5 / 3.4 54A69T (GCG→ACG) astEsuccinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (19/31);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

GCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACACCATCGCCCAGCCAGCGCCAGCGAA  >  NZ_CP009273/1820847‑1821110
                                                                                                                                 |                                                                                                                                      
gCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCaa                                                                                                                               <  2:513010/139‑1 (MQ=255)
 cTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAAc                                                                                                                              >  1:1143477/1‑139 (MQ=255)
       aCGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTcc                                                                                                                        <  2:804348/139‑1 (MQ=255)
       aCGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTGGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTcc                                                                                                                        <  2:1043393/139‑1 (MQ=255)
         gACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAg                                                                                                                      >  2:752287/1‑139 (MQ=255)
          aCCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGc                                                                                                                     <  1:78131/139‑1 (MQ=255)
          aCCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGc                                                                                                                     <  1:524301/139‑1 (MQ=255)
          aCCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGc                                                                                                                     <  2:760311/139‑1 (MQ=255)
              ccAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCAtct                                                                                                                 <  2:1271872/139‑1 (MQ=255)
              ccAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCAtct                                                                                                                 >  1:837897/1‑139 (MQ=255)
                 aaCATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCca                                                                                                              <  1:53762/139‑1 (MQ=255)
                  aCATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCcac                                                                                                             <  2:85584/139‑1 (MQ=255)
                        gATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGc                                                                                                       <  1:51266/139‑1 (MQ=255)
                        gATTCATATCGCTATGGCAAAAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGc                                                                                                       <  2:199541/139‑1 (MQ=255)
                         aTTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGcc                                                                                                      <  2:346032/139‑1 (MQ=255)
                              tatCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGtct                                                                                                  >  2:391333/1‑138 (MQ=255)
                                  gCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTa                                                                                             <  1:197838/139‑1 (MQ=255)
                                  gCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTa                                                                                             <  1:148758/139‑1 (MQ=255)
                                    tATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTAcc                                                                                           >  1:475423/1‑139 (MQ=255)
                                           aTAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTAtt                                                                                    >  1:537016/1‑139 (MQ=255)
                                            tAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTc                                                                                   >  1:487424/1‑139 (MQ=255)
                                               cGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCg                                                                                <  1:1237297/139‑1 (MQ=255)
                                               cGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCg                                                                                <  2:487424/139‑1 (MQ=255)
                                                   tCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTg                                                                             >  2:298900/1‑138 (MQ=255)
                                                     ccTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGaa                                                                           <  1:752287/138‑1 (MQ=255)
                                                                     gaggaTTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCAt                                                                                                                   <  2:465716/82‑1 (MQ=255)
                                                                     gaggaTTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCAt                                                                                                                   >  1:465716/1‑82 (MQ=255)
                                                                     gaggaTTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTg                                                          <  2:837897/139‑1 (MQ=255)
                                                                      aggaTTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGc                                                         <  2:537016/139‑1 (MQ=255)
                                                                       ggaTTCCCGAGGATCCCCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCg                                                        <  1:981752/139‑1 (MQ=255)
                                                                        gaTTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCgg                                                       <  2:475423/139‑1 (MQ=255)
                                                                          ttCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGc                                                                                                       <  1:93721/89‑1 (MQ=255)
                                                                          ttCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGc                                                                                                       >  2:93721/1‑89 (MQ=255)
                                                                            cccGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTgg                                                    <  1:298900/138‑1 (MQ=255)
                                                                               gAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGt                                                 >  2:574313/1‑138 (MQ=255)
                                                                                       ccAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGt                                        <  1:680931/139‑1 (MQ=255)
                                                                                             aCCGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTc                                  >  2:885869/1‑139 (MQ=255)
                                                                                               cGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCa                                                                                                <  2:307766/75‑1 (MQ=255)
                                                                                               cGCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCa                                                                                                >  1:307766/1‑75 (MQ=255)
                                                                                                gCCAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAg                               >  2:900697/1‑139 (MQ=255)
                                                                                                  cAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAAc                                                                                                                              >  2:1175856/1‑42 (MQ=255)
                                                                                                  cAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAAc                                                                                                                              <  1:1175856/42‑1 (MQ=255)
                                                                                                  cAACGTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGca                             <  2:944545/139‑1 (MQ=255)
                                                                                                      gTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACAcc                         <  2:1143477/139‑1 (MQ=255)
                                                                                                      gTAACGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACAcc                         <  2:101213/139‑1 (MQ=255)
                                                                                                        aaCGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTc                                       <  2:334198/123‑1 (MQ=255)
                                                                                                        aaCGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTc                                       >  1:334198/1‑123 (MQ=255)
                                                                                                         aCGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCAtc                                                                                                                  <  2:607841/47‑1 (MQ=255)
                                                                                                         aCGGGATCTCGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCAtc                                                                                                                  >  1:607841/1‑47 (MQ=255)
                                                                                                               tctcGCCGTGAGATATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACACCATCGCccag                >  2:494967/1‑139 (MQ=255)
                                                                                                                         agaTATCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACACCATCGCCCAgccagcgccag      >  1:1263501/1‑139 (MQ=255)
                                                                                                                         agaTATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACACCATCGCCCAgccagcgccag      >  1:834567/1‑139 (MQ=255)
                                                                                                                           atatCGTGCCAAGCAACGCGTCCAGCATCTACACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGCGGCGTTAATGGCGTCAGTTCCAGCACACCATCGCCCAgccagcgccagc     >  2:288422/1‑138 (MQ=255)
                                                                                                                              tCGTGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACACCATCGCCCAGCCAGCGCCAGCGaa  <  1:880632/138‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                      
GCTGCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGCCAAGCAACGCGTCCAGCATCTCCACAGGTGCCGTCTCATTACCGTGTATTCCCGCTGAAATCACCAGTGCGCCTTGCGGTGGCGTTAATGGCGTCAGTTCCAGCACACCATCGCCCAGCCAGCGCCAGCGAA  >  NZ_CP009273/1820847‑1821110

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