Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,538 A→T 13.3% intergenic (+114/+28) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,5380AT13.3% 71.9 / 5.7 30intergenic (+114/+28)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (11/15);  new base T (2/2);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

GGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTA  >  NC_002947/4338395‑4338676
                                                                                                                                                 |                                                                                                                                          
ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                      <  2:125610/150‑1 (MQ=255)
 gCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCcgcg                                                                                                                                     >  1:506638/1‑150 (MQ=255)
   aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGaa                                                                                                                                   >  2:347788/1‑150 (MQ=255)
    aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                  <  2:506638/150‑1 (MQ=255)
     cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTTTACCCGCGAAAg                                                                                                                                 >  2:427700/1‑150 (MQ=255)
     cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                 <  1:347788/150‑1 (MQ=255)
        aCTTAAAGGGCATGTGCGGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTTAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                                                                                                              <  2:415454/150‑1 (MQ=255)
          ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGCCg                                                                                                                            <  1:427700/150‑1 (MQ=255)
                                cGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCaa                                                                                                      <  1:334565/150‑1 (MQ=255)
                                       cAATAGTCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAAGCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                               <  1:7407/150‑1 (MQ=255)
                                        aaTAGACCCGCCTGGGCGGGTCTTTTGCGTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTAGCGGGTAAACCCGCGAAAGGGCCGGCAAATGCTGCTCAGGGCTCAAACCCTTcc                                                                                              >  1:395589/1‑150 (MQ=255)
                                                 gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGCATCTGCCTCGATCTGCTGTTTGCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATc                                                                                     <  1:391510/150‑1 (MQ=255)
                                                                  gCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAc                                                                    >  2:350593/1‑150 (MQ=255)
                                                                   cTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAACCTTCCACAATGATCACCTCAGCCTTCGCCAcc                                                                   >  1:452935/1‑150 (MQ=255)
                                                                    ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCa                                                                                                              <  2:69452/106‑1 (MQ=255)
                                                                    ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCa                                                                                                              >  1:69452/1‑106 (MQ=255)
                                                                         tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTgg                                                             <  1:350593/150‑1 (MQ=255)
                                                                              ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCCCAAACCCTTCCACAATGAGCACCTCAGCCTTGGCCACCCCTGGGCGCTg                                                        >  1:380162/1‑150 (MQ=255)
                                                                                   cAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                                                                                                >  1:355448/1‑105 (MQ=255)
                                                                                   cAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                                                                                                <  2:355448/105‑1 (MQ=255)
                                                                                       ccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGc                                               <  1:11977/150‑1 (MQ=255)
                                                                                                  aaTGAATCTGCCAGCGGTGCGGCAGGT‑‑CACGTCATTCGGCGTGTAAACCTGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACa                                    <  1:494196/150‑1 (MQ=255)
                                                                                                           gCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAGTCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGc                           >  1:392496/1‑150 (MQ=255)
                                                                                                            ccAGCGGTGCGGCAGGC‑‑CCGGCCCTTTCGCGCGTAAACCCGCGCGAGGAGCAGGTAAGCCTGCTCAGGGCTCAACCCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg                          <  2:395589/150‑1 (MQ=255)
                                                                                                              aGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCcaca                                                                                           >  2:15786/1‑83 (MQ=255)
                                                                                                              aGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCcaca                                                                                           <  1:15786/83‑1 (MQ=255)
                                                                                                                         agat‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACg             <  1:266778/147‑1 (MQ=255)
                                                                                                                                      tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAt                                                                                      >  2:15665/1‑66 (MQ=37)
                                                                                                                                      tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAt                                                                                      <  1:15665/66‑1 (MQ=37)
                                                                                                                                      tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa  >  2:496872/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
GGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGT‑‑CCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTA  >  NC_002947/4338395‑4338676

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