Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 183,409 A→G F577L (TTT→CTT)  glnD ← bifunctional uridylyltransferase/uridylyl‑removing protein GlnD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273183,4090AG97.5% 137.0 / ‑4.0 40F577L (TTT→CTT) glnDbifunctional uridylyltransferase/uridylyl‑removing protein GlnD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (19/20);  total (19/21)
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.89e-01

CAAAGTAGAGCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGTTCAGCCCGTGGAGTTCGGCAAAATGCA  >  NZ_CP009273/183280‑183540
                                                                                                                                 |                                                                                                                                   
cAAAGTAGAGCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTtgat                                                                                                                            <  2:990745/139‑1 (MQ=255)
       gagCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTcc                                                                                                                     >  1:349761/1‑139 (MQ=255)
          cTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCggg                                                                                                                  >  2:393355/1‑139 (MQ=255)
             aCGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTcc                                                                                                               >  1:660724/1‑139 (MQ=255)
             aCGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTc                                                                                                                <  1:638921/138‑1 (MQ=255)
             aCGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTc                                                                                                                <  2:129558/138‑1 (MQ=255)
               gCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTcc                                                                                                                     >  1:527612/1‑131 (MQ=255)
               gCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTcc                                                                                                                     <  2:527612/131‑1 (MQ=255)
               gCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCt                                                                                                              >  2:107788/1‑138 (MQ=255)
                      aCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATcg                                                                                                      <  1:166529/139‑1 (MQ=255)
                      aCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATcg                                                                                                      <  2:530939/139‑1 (MQ=255)
                           gCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCgcggcg                                                                                                 >  2:140852/1‑139 (MQ=255)
                               ccAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCCTGATGACTTCCGGGTCCTGAATATCGCGGCGTTg                                                                                              <  1:449650/138‑1 (MQ=255)
                                                tttCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCaa                                                                            <  2:172154/139‑1 (MQ=255)
                                                          ggCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAAc                                                                           <  2:519582/130‑1 (MQ=255)
                                                          ggCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAAc                                                                           >  1:519582/1‑130 (MQ=255)
                                                           gcgcAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGc                                                                 >  2:429195/1‑139 (MQ=255)
                                                                       gCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAgg                                                     >  2:1071498/1‑139 (MQ=255)
                                                                           caGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGAc                                                 <  1:670297/139‑1 (MQ=255)
                                                                             gTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGTTGCTGGCGAACCAGCCAGGCGAc                                                 >  1:707661/1‑137 (MQ=255)
                                                                                aGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGc                                             <  2:457688/138‑1 (MQ=255)
                                                                                      aCCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGAc                                                 >  2:721695/1‑128 (MQ=255)
                                                                                      aCCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGAc                                                 <  1:721695/128‑1 (MQ=255)
                                                                                       ccAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCg                                                                                   >  2:476762/1‑93 (MQ=255)
                                                                                       ccAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCg                                                                                   <  1:476762/93‑1 (MQ=255)
                                                                                          gATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTcgcg                                  <  1:955054/139‑1 (MQ=255)
                                                                                           aTAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTcc                                                                                                                     >  1:749297/1‑55 (MQ=255)
                                                                                           aTAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTcc                                                                                                                     <  2:749297/55‑1 (MQ=255)
                                                                                            tAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGt                                 >  1:1275337/1‑138 (MQ=255)
                                                                                             aGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTcgcg                                  >  2:1222777/1‑136 (MQ=255)
                                                                                             aGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTcgcg                                  <  1:1222777/136‑1 (MQ=255)
                                                                                             aGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGt                                 >  2:306064/1‑137 (MQ=255)
                                                                                               cgcAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGt                             <  1:1137744/139‑1 (MQ=255)
                                                                                               cgcAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGt                             >  1:711539/1‑139 (MQ=255)
                                                                                               cgcAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGt                             <  2:9459/139‑1 (MQ=255)
                                                                                                 cAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGtt                            >  2:771387/1‑138 (MQ=255)
                                                                                                 cAGACGATTTTCCGTTTGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGTTc                           <  2:1003992/139‑1 (MQ=255)
                                                                                                                 tGCACTTCTTCGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGTTCAGCCCGTGGAGTTCgg           >  2:318676/1‑139 (MQ=255)
                                                                                                                           cGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGTTCAGCCCGTGGAGTTCGGCAAAATGCa  <  2:385393/138‑1 (MQ=255)
                                                                                                                           cGGCAAGCTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGTTCAGCCCGTGGAGTTCGGCAAAATGCa  <  2:711539/138‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                   
CAAAGTAGAGCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCCTGAATATCGCGGCGTTGGGCGGTCACCGACATCAACAGGTGCTGGCGAACCAGCCAGGCGACCAGCTGTGTTTCGCGTGAGTTCAGCCCGTGGAGTTCGGCAAAATGCA  >  NZ_CP009273/183280‑183540

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