Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,376,277 T→C intergenic (+41/+2) dcuD → / ← sspB putative transporter DcuD/ClpXP protease specificity‑enhancing factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,376,2770TC95.5% 68.3 / ‑3.8 22intergenic (+41/+2)dcuD/sspBputative transporter DcuD/ClpXP protease specificity‑enhancing factor
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (10/11);  total (10/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 = 9.33e-01

ACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGAT  >  NC_000913/3376171‑3376411
                                                                                                          |                                                                                                                                      
aCAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGaccac                                                                                                       <  1:323590/140‑1 (MQ=255)
 cAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGaccacc                                                                                                      >  1:368986/1‑140 (MQ=255)
      aTCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCgtg                                                                                                 >  2:325499/1‑140 (MQ=255)
      aTCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCgtg                                                                                                 >  1:277887/1‑140 (MQ=255)
          ccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGtggtgg                                                                                             <  2:368986/140‑1 (MQ=255)
                 aGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGa                                                                                      <  2:306705/140‑1 (MQ=255)
                       tttCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCa                                                                                <  2:277887/140‑1 (MQ=255)
                         tCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATc                                                                              <  2:17619/140‑1 (MQ=255)
                                         tGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtca                                                              <  1:241266/140‑1 (MQ=255)
                                         tGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtca                                                              <  1:144065/140‑1 (MQ=255)
                                           ccACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtcatc                                                            >  2:176128/1‑140 (MQ=255)
                                            cacaATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGtcatca                                                           <  1:8215/140‑1 (MQ=255)
                                              caATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCatcatc                                                         <  1:176128/140‑1 (MQ=255)
                                                                   tCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                    >  2:201240/1‑140 (MQ=255)
                                                                   tCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                    >  1:388990/1‑140 (MQ=255)
                                                                      ttAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                 >  1:130913/1‑140 (MQ=255)
                                                                                                     ttttGCATTACTTCACAACGCGTAATGCCGGTCGACCGCCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGAt  <  1:8110/140‑1 (MQ=255)
                                                                                                     ttttGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCg              <  2:17131/128‑1 (MQ=255)
                                                                                                     ttttGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCg              >  1:327129/1‑128 (MQ=255)
                                                                                                     ttttGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCg              <  2:327129/128‑1 (MQ=255)
                                                                                                     ttttGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCg              >  1:17131/1‑128 (MQ=255)
                                                                                                     ttttGCATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGAt  >  2:8110/1‑140 (MQ=255)
                                                                                                          |                                                                                                                                      
ACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGAT  >  NC_000913/3376171‑3376411

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