Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,483,443 T→A 17.9% L89Q (CTG→CAG)  puuA‑I → glutamate‑putrescine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,483,4430TA17.9% 58.2 / 9.7 28L89Q (CTG→CAG) puuA‑Iglutamate‑putrescine ligase
Reads supporting (aligned to +/- strand):  ref base T (13/10);  new base A (3/2);  total (16/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 = 1.53e-01

CCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCATT  >  NC_002947/2483302‑2483591
                                                                                                                                             |                                                                                                                                                    
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                                                                                                                                  <  2:155383/132‑1 (MQ=255)
                              gTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGc                                                                                                                                  >  1:155383/1‑132 (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                                                                            >  1:789872/1‑116 (MQ=255)
                                                                                                    cTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCa                                                                            <  2:789872/116‑1 (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)
                                                                                                                                            cTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCAtt  >  1:552228/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCATGGACATGCTGGAAAAAGTCGCAGCCGGTGCGCCTCTGAAGCTGCCGCAGAACTGCGTGCTGCTGGGTACCCAGGGCGGCCTTTACCCGATTGGTGACTATTGCTTCGCCGACGGCGACCCGGATGCTGCCCGGCGCCTGGTGCCGGGTACGCTGAAGCCGGTGCGCTGGGAAAAGGAAGTGATCGCCCAGATGCTGATCACGTCCGACGGCACTGCAAAGCCGATCGAATTCGAGCCTCGTGAAGTGCTGGCGCGTGTTCTTCAGCGTCTGGCCAAGCGTGGCATT  >  NC_002947/2483302‑2483591

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: