Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,003,978 A→G E140G (GAG→GGG)  yigL → sugar/pyridoxal phosphate phosphatase YigL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,003,9780AG96.5% 99.1 / ‑6.4 29E140G (GAG→GGG) yigLsugar/pyridoxal phosphate phosphatase YigL
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/12);  minor base T (0/1);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAACGACAATCCGGACATCATTACTAACGTTTATCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGAGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTG  >  NZ_CP009273/4003843‑4004111
                                                                                                                                       |                                                                                                                                     
tCAACGACAATCCGGACATCATTACTAACGTTTATCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGcc                                                                                                                                    >  2:96223/1‑139 (MQ=255)
          tCCGGACATCATTACTAACGTTTATCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGCGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGt                                                                                                                           >  1:454043/1‑138 (MQ=255)
          tCCGGACATCATTACTAACGTTTATCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGt                                                                                                                           >  2:523414/1‑138 (MQ=255)
                           aCGTTTATCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCAc                                                                                                         >  2:112448/1‑139 (MQ=255)
                              tttATCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTg                                                                                                      <  2:296222/139‑1 (MQ=255)
                                  tCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGa                                                                                                   <  1:465124/138‑1 (MQ=255)
                                                    tATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAActgct                                                                                 <  1:96223/138‑1 (MQ=255)
                                                            gccatcgccCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGATg                                                                          >  2:36500/1‑137 (MQ=255)
                                                                         agagATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAAc                                                           <  1:347582/139‑1 (MQ=255)
                                                                             aTGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGc                                                                            >  2:88062/1‑118 (MQ=255)
                                                                             aTGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGc                                                                            <  1:88062/118‑1 (MQ=255)
                                                                              tGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGTGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCg                                                      <  1:420605/139‑1 (MQ=255)
                                                                                  cTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTgg                                                  >  2:569077/1‑139 (MQ=255)
                                                                                                     ttCCAATATGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTc                                <  1:434960/138‑1 (MQ=255)
                                                                                                           tatGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTtct                         <  1:523414/139‑1 (MQ=255)
                                                                                                             tGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACg                                                          <  1:218026/104‑1 (MQ=255)
                                                                                                             tGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACg                                                          >  2:218026/1‑104 (MQ=255)
                                                                                                             tGCGCTGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTtctc                        >  1:406939/1‑138 (MQ=255)
                                                                                                                  tGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCt                   >  1:415548/1‑138 (MQ=255)
                                                                                                                  tGTATGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCt                   >  1:384965/1‑138 (MQ=255)
                                                                                                                      tGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAAcc              <  1:290201/139‑1 (MQ=255)
                                                                                                                      tGAGCCTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAAc               <  2:406939/138‑1 (MQ=255)
                                                                                                                          ccTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCACCGTCAGTTTCTCTACCTTAACCTGTc          >  2:94968/1‑139 (MQ=255)
                                                                                                                          ccTGGATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTc          >  2:456748/1‑139 (MQ=255)
                                                                                                                              gATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTgg       >  1:216922/1‑138 (MQ=255)
                                                                                                                              gATTACTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTgg       >  1:35909/1‑138 (MQ=255)
                                                                                                                                   cTGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTg  <  1:79772/138‑1 (MQ=255)
                                                                                                                                    tGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGAt                                            >  2:205016/1‑95 (MQ=255)
                                                                                                                                    tGGGGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGAt                                            <  1:205016/95‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
TCAACGACAATCCGGACATCATTACTAACGTTTATCGCGACGACGAATGGTTTATGAATCGCCATCGCCCGGAAGAGATGCGCTTTTTTAAAGAAGCGGTGTTCCAATATGCGCTGTATGAGCCTGGATTACTGGAGCCGGAAGGCGTCAGCAAAGTGTTCTTCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTG  >  NZ_CP009273/4003843‑4004111

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: