Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,483,440 G→A 18.6% R88H (CGC→CAC)  puuA‑I → glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,483,4400GA18.6% 45.3 / 7.4 27R88H (CGC→CAC) puuA‑Iglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base G (12/10);  new base A (3/2);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

CCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGC  >  NC_002947/2483302‑2483583
                                                                                                                                          |                                                                                                                                               
ccccggACATGCTGGAAGAAGTCGCAGCCGGTTGGCCTCTGAAGCTGCCGCAGAGCTGCGTGCTGGTGGGTACCCAGGGCGGCCTTTTCCGGTTTGGTGACTATTGCTTCGCCGACGGCGCCCCGGATGCTGCCCGGCGCCTGGTGCCgg                                                                                                                                      <  2:294380/146‑1 (MQ=255)
   tGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTACCCGGCACCAGGCGCCgggca                                                                                                                                   <  1:310613/150‑3 (MQ=255)
            tGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCGAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCTACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                          >  1:259692/1‑150 (MQ=255)
            tGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                          <  1:229863/150‑1 (MQ=255)
                 aaaGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTg                                                                                                                     <  2:462984/150‑1 (MQ=255)
                              gTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                          >  1:155383/1‑132 (MQ=255)
                              gTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                          <  2:155383/132‑1 (MQ=255)
                                 cgcCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAg                                                                                                     <  2:545007/150‑1 (MQ=255)
                                    ctctGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCAGGCGCCTGGTGCCGGGTACGCTCAAGCCGGTGCGCCGGGAAAAGGAAGTGa                                                                                                  >  2:158264/1‑150 (MQ=255)
                                           gCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCTTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCa                                                                                           >  1:96229/1‑150 (MQ=255)
                                              gcagCAGAACTGCGTGCTCCTGGGTACCCAGGGCAGCCTTTCCCCGCTTAGTGACGCTCTCGGCGCCGACGGCGACCCGGATGTTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATAGCCCAGAt                                                                                        <  2:285640/147‑1 (MQ=255)
                                                cGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCAGCAACGCGCCCGGCGCGTTGGGCCGGGTACGCTGAAGACGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGc                                                                                      >  1:928426/1‑150 (MQ=255)
                                                        tGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTTGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACg                                                                              >  1:276183/1‑150 (MQ=255)
                                                          cGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTg                                                                                                   <  1:662144/127‑1 (MQ=255)
                                                          cGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTg                                                                                                   >  2:662144/1‑127 (MQ=255)
                                                                              gcggcCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGa                                                        <  2:593914/150‑1 (MQ=255)
                                                                                     ttACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAAtt                                                 <  2:484263/150‑1 (MQ=255)
                                                                                         ccGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAg                                             >  2:1065988/1‑150 (MQ=255)
                                                                                                    cTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCa                                                                    <  2:789872/116‑1 (MQ=255)
                                                                                                    cTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCa                                                                    >  1:789872/1‑116 (MQ=255)
                                                                                                    cTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAgg                                  >  2:723680/1‑149 (MQ=255)
                                                                                                      aTTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGc                                >  1:1004784/1‑150 (MQ=255)
                                                                                                                          ccGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGACGGTGCGCTGGGCAAAGGAAGTGATCGCCCAGATGCTGATGACGTCCGACGGCACTGCAAAGACGATCAAATTCGAGCCTCGTGACGTCCGGGCGCGTGTTCTGGAGCGt            >  2:167887/1‑150 (MQ=255)
                                                                                                                                 gtacccGGCACCAGGCGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCa     >  2:133524/4‑150 (MQ=255)
                                                                                                                                 gtacccGGCACCAGGCGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCa     >  2:1135931/4‑150 (MQ=255)
                                                                                                                                    cccGGCACCAGGCGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGc  >  2:346710/1‑150 (MQ=255)
                                                                                                                                    cccGGCACCAGGCGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGc  <  1:133524/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                               
CCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGC  >  NC_002947/2483302‑2483583

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