Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,196,972 A→C H1072P (CAC→CCC)  BW25113_RS11045 → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,196,9720AC96.4% 88.3 / ‑3.1 28H1072P (CAC→CCC) BW25113_RS11045DUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (13/14);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GGGGTTTATCCGGCAAATGAACCGCGTTGTTTACTCAAAGCCTTTCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCACCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGTGCTGCTCACACCTGACGAGTCAACCAGTAA  >  NZ_CP009273/2196845‑2197101
                                                                                                                               |                                                                                                                                 
ggggTTTATCCGGCAAATGAACCGCGTTGTTTACTCAAAGCCTTTCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGaaa                                                                                                                        >  1:88894/1‑139 (MQ=255)
            gCAAATGAACCGCGTTGTTTACTCAAAGCCTTTCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGaaa                                                                                                            >  1:113685/1‑139 (MQ=255)
                             tttACTCAAAGCCTTTCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGca                                                                                           <  2:795513/139‑1 (MQ=255)
                             tttACTCAAAGCCTTTCGTGTGGCCGCAGAGGGGGAGATCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGAACCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGca                                                                                           >  2:498932/1‑139 (MQ=255)
                               tACTCAAAGCCTTTCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGcac                                                                                          <  2:528139/138‑1 (MQ=255)
                                       gCCTTTCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTg                                                                                 <  2:302947/139‑1 (MQ=255)
                                            tCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCTGAt                                                                            >  2:158855/1‑139 (MQ=255)
                                                   gCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGaaa                                                                     >  2:743598/1‑139 (MQ=255)
                                                                    cTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGc                                                    >  2:726772/1‑139 (MQ=255)
                                                                      gCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGc                                                    >  2:386738/1‑137 (MQ=255)
                                                                       cAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGca                                                                                           >  1:582054/1‑97 (MQ=255)
                                                                       cAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGca                                                                                           <  2:582054/97‑1 (MQ=255)
                                                                       cAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAg                                                 >  2:360825/1‑139 (MQ=255)
                                                                          tGCGCAAGATGAGCCAATTGGCCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATg                                                                                                          <  2:85028/79‑1 (MQ=255)
                                                                          tGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATg                                                                                                          >  1:85028/1‑79 (MQ=255)
                                                                            cgcAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTc                                            >  2:151127/1‑139 (MQ=255)
                                                                            cgcAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTc                                            >  2:626525/1‑139 (MQ=255)
                                                                             gcAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTcg                                           <  1:626525/139‑1 (MQ=255)
                                                                                     aGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCa                                    <  2:46025/138‑1 (MQ=255)
                                                                                     aGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCa                                    <  2:549075/138‑1 (MQ=255)
                                                                                     aGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCa                                    <  1:657089/138‑1 (MQ=255)
                                                                                     aGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCAc                                   >  2:222463/1‑139 (MQ=255)
                                                                                         aaTTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCACCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGtg                               <  2:472492/139‑1 (MQ=255)
                                                                                              aCCTGCCTGCGGACGCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGTgctgct                          <  1:769890/139‑1 (MQ=255)
                                                                                                            gCTCTGATTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGTGCTGCTCACACCTGACGAGt            <  1:386738/139‑1 (MQ=255)
                                                                                                                  aTTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGTGCTGCTCACACCTGACGAGTCAAc        <  1:743598/137‑1 (MQ=255)
                                                                                                                  aTTGGCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGTGCTGCTCACACCTGACGAGTCAAc        <  1:360825/137‑1 (MQ=255)
                                                                                                                      gCATTGCCCCCCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGTGCTGCTCACACCTGACGAGTCAACCAGTaa  <  2:72971/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                 
GGGGTTTATCCGGCAAATGAACCGCGTTGTTTACTCAAAGCCTTTCGTGTGGCCGCAGAGGGGGAGTTCTGCAATGCGCAAGATGAGCCAATTGACCTGCCTGCGGACGCTCTGATTGGCATTGCCCACCCGTTAGAAATGACAGCAGAAATGCGCAGTGAATTTGCACAGCTTTTTGCCGATTACGAAATTATGCCGCCTTTTCGCCAGTTGTCGCGCCGCACGGTGCTGCTCACACCTGACGAGTCAACCAGTAA  >  NZ_CP009273/2196845‑2197101

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