Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 731,616 C→T H678H (CAC→CAT)  rhsC → rhs element protein RhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913731,6160CT98.0% 129.3 / NA 49H678H (CAC→CAT) rhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/20);  minor base A (1/0);  total (29/20)
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.00e+00

GATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGA  >  NC_000913/731494‑731747
                                                                                                                          |                                                                                                                                   
gATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCaa                                                                                                                    >  2:691567/1‑140 (MQ=17)
gATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCaa                                                                                                                    >  2:645021/1‑140 (MQ=17)
gATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCaa                                                                                                                    >  2:551385/1‑140 (MQ=17)
gATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCaa                                                                                                                    >  1:573030/1‑140 (MQ=17)
gATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGACGCCCACAGCGTAATCAACACGACCTGACCGGCaa                                                                                                                    >  2:307016/1‑140 (MQ=255)
   gctgcCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACt                                                                                                                 >  2:705268/1‑140 (MQ=17)
     tgcCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTa                                                                                                               >  2:563931/1‑140 (MQ=17)
       ccGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATc                                                                                                             >  2:641534/1‑140 (MQ=17)
           aCGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCa                                                                                                         <  2:531121/140‑1 (MQ=17)
               gTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCTCACACAGCGTTAGCATAACGACCTGCCCGGAAAACTTATCCGCAGCGa                                                                                                     >  2:303464/1‑140 (MQ=11)
                        cTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAggg                                                                                            >  1:798009/1‑140 (MQ=17)
                         tGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAgggg                                                                                           >  1:431620/1‑140 (MQ=18)
                              aGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGt                                                                                      <  2:573030/140‑1 (MQ=25)
                               gTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTc                                                                                     <  2:609261/140‑1 (MQ=21)
                                     aCGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAGACTTATCTGCATCGTGGATGACGGGCTGGTCACCCAc                                                                               >  1:286035/1‑140 (MQ=11)
                                      cGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACt                                                                              <  2:431620/140‑1 (MQ=25)
                                       ggCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTg                                                                             >  2:161934/1‑140 (MQ=25)
                                          aGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGc                                                                                          <  2:162517/124‑1 (MQ=21)
                                          aGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGc                                                                                          >  1:162517/1‑124 (MQ=25)
                                          aGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCAcc                                                                                  >  1:34178/1‑132 (MQ=25)
                                          aGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCAcc                                                                                  <  2:34178/132‑1 (MQ=21)
                                           gCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCAc                                                                         <  1:705268/140‑1 (MQ=25)
                                                accaccTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGa                                                                    <  1:799462/140‑1 (MQ=25)
                                                         cGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGATGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGa                                                           <  2:187650/140‑1 (MQ=25)
                                                         cGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGa                                                           >  2:500230/1‑140 (MQ=25)
                                                          gTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGAc                                                          >  1:829899/1‑140 (MQ=25)
                                                           ttACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGAcc                                                         <  1:645021/140‑1 (MQ=21)
                                                           ttACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGAcc                                                         >  1:704611/1‑140 (MQ=25)
                                                           ttACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGAcc                                                         <  1:161934/140‑1 (MQ=25)
                                                            tACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTAACCGACTGGCACTGTGACGAAGCAGACCg                                                        >  1:228390/1‑140 (MQ=255)
                                                               gATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGa                                                           >  1:342311/1‑134 (MQ=25)
                                                               gATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGa                                                           <  2:342311/134‑1 (MQ=25)
                                                               gATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCt                                                     >  1:611920/1‑140 (MQ=25)
                                                                   tACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcg                                                 <  1:154794/140‑1 (MQ=21)
                                                                     cTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGTGGGGCTGGTCACGCACGGGCACTATGACGAGGCAGACCGCCTCAcgca                                               >  2:219651/1‑140 (MQ=11)
                                                                           cGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgcaccgcac                                         >  1:111378/1‑140 (MQ=17)
                                                                                            cGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCg                        >  2:203427/1‑140 (MQ=37)
                                                                                            cGGCTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCg                        >  2:571741/1‑140 (MQ=37)
                                                                                               cTTTGACGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGTACCGTGAATGGCGAAACCGCag                     >  1:316421/1‑140 (MQ=255)
                                                                                                  tGACGGCCGCACACACCGTTATGATTACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGc                  <  1:307016/140‑1 (MQ=255)
                                                                                                    aCGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCgg                <  2:111378/140‑1 (MQ=255)
                                                                                                    aCGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCgg                <  2:614280/140‑1 (MQ=38)
                                                                                                    aCGGCCGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCAACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCgg                >  2:645733/1‑140 (MQ=255)
                                                                                                      ggCCGCACACAGCGTTAACATCTCGACCTGAACGGAAAAGTTATCCGCCGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTg              <  2:228390/140‑1 (MQ=11)
                                                                                                       gCCGCACACAGCGTTATCATCACGACCTGACAGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTgg             <  2:286035/140‑1 (MQ=255)
                                                                                                        ccGCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAAACCGCCTCACGCACCGCACCGTGAAAGGCGATACCGCAGAGCGGTGGc            >  2:354840/1‑140 (MQ=255)
                                                                                                          gCACACAGCGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAg          <  1:563931/140‑1 (MQ=255)
                                                                                                                aGCGTTATCATCACGACCTGACCGGCAAAATTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGAc    <  1:303464/140‑1 (MQ=255)
                                                                                                                  cGTTATCATCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGa  >  2:76380/1‑140 (MQ=37)
                                                                                                                          |                                                                                                                                   
GATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGA  >  NC_000913/731494‑731747

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑