Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,174,381 T→G C270G (TGC→GGC)  yegU → ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,174,3810TG97.1% 111.5 / ‑3.8 34C270G (TGC→GGC) yegUADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (15/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

GTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGG  >  NZ_CP009273/2174300‑2174463
                                                                                 |                                                                                  
gTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTc                                                                            >  2:617231/1‑90 (MQ=255)
gTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTc                                                                            >  1:573467/1‑90 (MQ=255)
                   cTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcg                                                         >  2:277369/1‑90 (MQ=255)
                     aTTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcggc                                                       <  1:278621/90‑1 (MQ=255)
                      ttGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGgcggcg                                                      <  1:771620/90‑1 (MQ=255)
                      ttGAGTCCGTTCCGTGCGCCATTGCGCTAGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGgcggcg                                                      <  2:512942/90‑1 (MQ=255)
                        gAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGac                                                    >  1:733655/1‑90 (MQ=255)
                          gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                  >  2:872254/1‑90 (MQ=255)
                          gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                  <  2:807118/90‑1 (MQ=255)
                            ccgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAg                                                <  1:196929/90‑1 (MQ=255)
                               ttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGaca                                             >  1:270741/1‑90 (MQ=255)
                                   gtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCtg                                                                         >  1:218826/1‑58 (MQ=255)
                                   gtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCtg                                                                         <  2:218826/58‑1 (MQ=255)
                                    tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                        <  2:77602/90‑1 (MQ=255)
                                    tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                        <  1:516172/90‑1 (MQ=255)
                                    tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                        <  1:774860/90‑1 (MQ=255)
                                                cTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCg                            >  2:201900/1‑90 (MQ=255)
                                                  ggTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGAc                          >  2:526919/1‑90 (MQ=255)
                                                            gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgc                >  2:978913/1‑90 (MQ=255)
                                                            gCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgc                >  2:994660/1‑90 (MQ=255)
                                                              acaAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgcgg              <  1:164504/90‑1 (MQ=255)
                                                                aaaCCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTgcggcg            <  1:238410/90‑1 (MQ=255)
                                                                 aaCCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                        >  2:468235/1‑61 (MQ=255)
                                                                 aaCCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                        <  1:468235/61‑1 (MQ=255)
                                                                     gACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTg       >  2:236425/1‑90 (MQ=255)
                                                                     gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                        <  2:585459/57‑1 (MQ=255)
                                                                     gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                        >  1:585459/1‑57 (MQ=255)
                                                                     gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTg       >  1:553364/1‑90 (MQ=255)
                                                                        ccGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAt    <  2:553364/90‑1 (MQ=255)
                                                                        ccGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAt    <  2:982315/90‑1 (MQ=255)
                                                                         cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg   <  2:3704/90‑1 (MQ=255)
                                                                         cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg   <  2:245862/90‑1 (MQ=255)
                                                                         cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg   <  2:636660/90‑1 (MQ=255)
                                                                          gAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATgg  >  1:518802/1‑90 (MQ=255)
                                                                                 |                                                                                  
GTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGG  >  NZ_CP009273/2174300‑2174463

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