Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,604,013 C→A Q77H (CAG→CAT ftsY ← signal recognition particle receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,604,0130CA96.2% 84.8 / ‑2.5 26Q77H (CAG→CATftsYsignal recognition particle receptor
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (10/15);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

ACCGCTTCGGCGTTGACGTCTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACCTGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGAACGATTTCTTCTACAAC  >  NC_000913/3603876‑3604144
                                                                                                                                         |                                                                                                                                   
aCCGCTTCGGCGTTGACGTCTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATg                                                                                                                                   >  1:319933/1‑140 (MQ=255)
aCCGCTTCGGCGTTGACGTCTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATg                                                                                                                                   >  2:903595/1‑140 (MQ=255)
    cTTCGGCGTTGACGTTTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACCTGTTCa                                                                                                                               <  2:525280/140‑1 (MQ=255)
      tCGGCGTTGACGTCTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGt                                                                                                                             <  2:559252/140‑1 (MQ=255)
       cGGCGTTGACGTCTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTg                                                                                                                            <  2:924087/140‑1 (MQ=255)
         gCGTTGACGTCTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGAc                                                                                                                          <  1:903595/140‑1 (MQ=255)
                      ccGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGt                                                                                                                             <  1:610827/124‑1 (MQ=255)
                      ccGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGt                                                                                                                             >  2:610827/1‑124 (MQ=255)
                                           gTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCg                                                                                        <  2:319933/140‑1 (MQ=255)
                                               gATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCtt                                                                                    <  2:282043/140‑1 (MQ=255)
                                                aTAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTc                                                                                   >  1:343796/1‑140 (MQ=255)
                                                aTAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTc                                                                                   >  1:918664/1‑140 (MQ=255)
                                                         ggCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATgcgc                                                                          >  1:100839/1‑140 (MQ=255)
                                                                                        gCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCagag                                           <  2:100839/140‑1 (MQ=255)
                                                                                           gtgCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGcc                                        <  2:918664/140‑1 (MQ=255)
                                                                                                  gACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAg                                 <  2:295856/140‑1 (MQ=255)
                                                                                                      tCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGctc                             >  2:656094/1‑140 (MQ=255)
                                                                                                        cGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCaga                                            <  2:175115/123‑1 (MQ=255)
                                                                                                        cGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCaga                                            >  1:175115/1‑123 (MQ=255)
                                                                                                               ggCTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTg                    <  2:952597/140‑1 (MQ=255)
                                                                                                                 cTGCGCTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGaa                  <  2:343796/140‑1 (MQ=255)
                                                                                                                      cTTTTTCACTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGAACGAtt             <  2:389348/140‑1 (MQ=255)
                                                                                                                              cTTTCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGAACGATTTCTTCTAc     >  1:491421/1‑140 (MQ=255)
                                                                                                                                 tCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGAACGATTTCTTCTACa    >  2:1051157/1‑138 (MQ=255)
                                                                                                                                 tCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGAACGATTTCTTCTACa    <  1:1051157/138‑1 (MQ=255)
                                                                                                                                 tCAGCAACATGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGAACGATTTCTTCTACAAc  <  2:491421/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
ACCGCTTCGGCGTTGACGTCTTCCGGCAACGGCAGCTCTTCACGTTCGATAGCCACTGGCTCCGGCGTTTCTTCTACGACCGGTTCCGGCTGTGCAACGACTTCCGCTTCAGGCTGCGCTTTTTCACTTTCAGCAACCTGTTCAGTGACTTCCACAACGTCGGCAGCAAAAGTTTCCGCCTCGGCTTCAGTATGCGCTTGCGGCTGCTCTTCAACGGCTTGTTCAGAGGCCTTCACAGGCTCTTGCGCCTGAACGATTTCTTCTACAAC  >  NC_000913/3603876‑3604144

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