Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,625,832 T→C Q21Q (CAA→CAG)  guaA ← glutamine‑hydrolyzing GMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,625,8320TC97.9% 173.2 / ‑3.6 47Q21Q (CAA→CAG) guaAglutamine‑hydrolyzing GMP synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (24/22);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TCTTCAGTAGTACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCA  >  NZ_CP009273/2625698‑2625965
                                                                                                                                      |                                                                                                                                     
tcttcAGTAGTACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGa                                                                                                                                    <  2:736864/138‑1 (MQ=255)
   tcAGTAGTACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAgtg                                                                                                                                 >  2:123600/1‑138 (MQ=255)
     agtagtACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTa                                                                                                                               <  2:498905/138‑1 (MQ=255)
       tagtaCTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACt                                                                                                                             >  1:518231/1‑138 (MQ=255)
         gtaCTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTgag                                                                                                                          <  2:291625/139‑1 (MQ=255)
               tCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCg                                                                                                                     <  2:494021/138‑1 (MQ=255)
               tCCGGGCAGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAAGTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCg                                                                                                                     >  1:410907/1‑138 (MQ=255)
                   ggCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTc                                                                                                                <  2:1120328/139‑1 (MQ=255)
                    gccgccGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTc                                                                                                                >  1:793297/1‑138 (MQ=255)
                         cgGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAAt                                                                                                          >  1:609035/1‑139 (MQ=255)
                                  taatGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTc                                                                                                                >  1:1263572/1‑124 (MQ=255)
                                  taatGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTc                                                                                                                <  2:1263572/124‑1 (MQ=255)
                                  taatGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcg                                                                                                 <  1:171622/139‑1 (MQ=255)
                                  taatGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcg                                                                                                 <  2:567167/139‑1 (MQ=255)
                                    atGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgat                                                                                               <  2:1041392/139‑1 (MQ=255)
                                          cTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTc                                                                                                                >  2:809553/1‑116 (MQ=255)
                                          cTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTc                                                                                                                <  1:809553/116‑1 (MQ=255)
                                                        aCGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGt                                                                            >  1:1254210/1‑138 (MQ=255)
                                                           aaTTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCAtt                                                                        <  2:793297/139‑1 (MQ=255)
                                                            aTTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCTGTAAACACCCAGCTCACGCACGCGGCGCGCAAGCAGCTGAGTGTACTGATTGCCGAAGTCCAGAATGAGGATGCTATGCTTATGAATGTTTTCCGTCATTg                                                                       <  2:1047584/139‑1 (MQ=255)
                                                                      cTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcg                                                                                                 >  1:1219856/1‑103 (MQ=255)
                                                                      cTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcg                                                                                                 <  2:1219856/103‑1 (MQ=255)
                                                                            caTCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGc                                                       <  2:187394/139‑1 (MQ=255)
                                                                                   cGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGCGAACCGAAGCCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGa                                                 >  1:491614/1‑138 (MQ=255)
                                                                                   cGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCg                                                                                                                     >  2:822799/1‑70 (MQ=255)
                                                                                   cGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCg                                                                                                                     <  1:822799/70‑1 (MQ=255)
                                                                                   cGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGaa                                                <  2:937735/139‑1 (MQ=255)
                                                                                   cGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGa                                                 >  1:412343/1‑138 (MQ=255)
                                                                                         caGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGa                                           >  1:716373/1‑138 (MQ=255)
                                                                                         caGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGa                                           >  1:34266/1‑138 (MQ=255)
                                                                                                      acacCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGcc                             >  2:516407/1‑139 (MQ=255)
                                                                                                      acacCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGcc                             <  1:1148551/139‑1 (MQ=255)
                                                                                                      acacCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGcc                             <  1:123600/139‑1 (MQ=255)
                                                                                                           cAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACa                        >  1:968440/1‑139 (MQ=255)
                                                                                                           cAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACa                        >  1:610492/1‑139 (MQ=255)
                                                                                                           cAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACa                        >  1:1067490/1‑139 (MQ=255)
                                                                                                             gCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCa                                                                                                              <  2:54853/51‑1 (MQ=255)
                                                                                                             gCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCa                                                                                                              >  1:54853/1‑51 (MQ=255)
                                                                                                             gCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACAt                       >  2:1277511/1‑138 (MQ=255)
                                                                                                             gCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACAt                       >  1:104895/1‑138 (MQ=255)
                                                                                                             gCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACAt                       >  1:82571/1‑138 (MQ=255)
                                                                                                             gCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATg                      <  2:553233/139‑1 (MQ=255)
                                                                                                                  cgcacgcGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTc                 >  1:131567/1‑139 (MQ=255)
                                                                                                                     acgcGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCggg              <  2:1254210/139‑1 (MQ=255)
                                                                                                                        cggcgCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAg                                            <  2:270597/106‑1 (MQ=255)
                                                                                                                        cggcgCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAg                                            >  1:270597/1‑106 (MQ=255)
                                                                                                                                 aCCAGCTGAGTGTACTGAGACCCGAAGTCCAGAATAAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGACTCa  <  2:410907/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TCTTCAGTAGTACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCA  >  NZ_CP009273/2625698‑2625965

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