Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,780,007 T→C E907G (GAA→GGA)  ypjA ← adhesin‑like autotransporter YpjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,780,0070TC94.1% 44.8 / ‑4.2 17E907G (GAA→GGA) ypjAadhesin‑like autotransporter YpjA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (10/6);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTTGTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCA  >  NC_000913/2779774‑2780253
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
ccGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGGACCTGCTAATACCAATAAATTACCGCCATTTCCCCGCAACATATTGGTc                                                                                                                                                                                                                                       >  1:218250/1‑251 (MQ=255)
ccGATGCATACGCCAGCGAACCTGCGTAGACGCTTGCAGTACAGCCAACGACCTGGAGAACATCTGCCCGGGACAGAGCCTGAACCTCATCATTTGATCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGAGGCGGAGGTCTGAACCAGGTTTTGAATTGACCCACCCTCGCCAACCGTGGTGTCGCGGGCTTCGGTACCTGCTAATAACATTAACTTACCGACAATTCCAAGCACCATCTTGGTc                                                                                                                                                                                                                                       >  2:305533/1‑251 (MQ=255)
            ccAGCGTTCCTGCGGAGACGATTGCCGTCCCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACGCATCTTTTGCCCGCCCAAGGTTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCTTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCg                                                                                                                                                                                                                           <  1:75587/251‑1 (MQ=255)
                                             cAGCGTCCGGGAGATCTTCGGCCCGGGCTAGGGCCTTAACCTCATCTTTGGACCGCCCAAGGGTCTTTTGCCCGCCAGAGTTAACCTTTGCGGCGGAGTCCTGCCCCAGGTTTTGCATTGCCCCCCCCTTGCCAACCGTGGAGTCGCGAGCTTCGGCACCTGCTAATACCAATAACTTACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCg                                                                                                                                                                                          <  2:61173/251‑1 (MQ=255)
                                                            cttctGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTCTATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCACCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATaacg                                                                                                                                                                           >  2:135492/1‑248 (MQ=255)
                                                             ttctGCCCGGGCCAGCGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGGAATGGAAAAAGTAACGTAATGGTGAGTACCGCTGATTTCAATAcagt                                                                                                                                                                          >  1:81221/1‑248 (MQ=255)
                                                                            atccTGAAACTCATCTTTTGACCGCCCAGGGGTATATGCCCCGCCAGATTTATCTGTTGTGGCGTAGTCTTGACCCAGGTTTTGCTTTGCCCCCCCCTTGCCAACCGTGGAGTCGCCACTTTCGGTACCTTCTAATACCAACAAATTACCGCCATTTCCCAGCAACATATTGTTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGACTACCGCTGATTTCAATACCGTTAGCCTTGCTCGTCt                                                                                                                                                           <  1:162338/249‑1 (MQ=255)
                                                                                                          ggTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCCTTTCCCAGCCACATATTGGTCGCTAACTTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGGGAGCGGCCCCGGTGTTCTGGTCGattt                                                                                                                             >  2:273492/1‑249 (MQ=255)
                                                                                                                   cccGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCtt                                                                                                                    >  1:144807/1‑251 (MQ=255)
                                                                                                                                      gtgGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCAACCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGt                                                                                                 >  2:218118/1‑251 (MQ=255)
                                                                                                                                                          ggttttcaatGCCCCACCCTTGCCAACCGTGGAATCGCGAGTTTCGGTACCTGCTAATACAAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGATGAGTACCGCTGATTTCAATACCGTTAGCCGTGCGCGTCTGGAGAGCGGCACCGCTGTTCTTGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGa                                                                             <  1:152285/242‑1 (MQ=255)
                                                                                                                                                                                                                      caataccttaacaccATTTCCCAGCAACATATTGGGCGCTAAATTGCCAGAAATGGAAAAAGTACAGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTTGTTTTTATTGCCTCTCCGTTAGCTGACACTGTTTGTATTCCGCCGTCGTTAAttg                 >  2:102352/14‑251 (MQ=255)
                                                                                                                                                                                                                              tACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCAtt                                                                               >  2:146705/1‑181 (MQ=255)
                                                                                                                                                                                                                              tACCGCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCAtt                                                                               <  1:146705/181‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCATTTCCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGAACCGTACTGGTGAGTACCGCTGATTTCAATAACGTTAGCCGTGCTCGTCTGGAGCGCGGCACCGTTGTTCTGGACGATATCTGCCGCTTTGCCCTTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTTGTTTTTCTTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCg     >  2:120394/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gCCATTTCCCAGCAACATATTGGTCGCTAAAGTGTCGGAAATGGAAAAAGTACCGTACTGGTGAGTGCCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTTGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGTCTTCATTGATCTTTGTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATGGTTGGCTCATTCg     <  2:214545/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTTCCCATCAACATATTGGTCGCTCAATAGCCGGAAATGGAAAATGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGGTAGCCGTGCTCGGCTGGAGATCGTCGCCGCTGGTCTGGACGATATCTTCCGCTTTCCCATTATCGTTTACTGGCAGCGTACCGCCTTTATTGATCTTTGTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTTTATGCCGCCGTCGTTAATTTTTGTCTCATTCGCca  <  1:305533/251‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
CCGATGCATCCGCCAGCGTACCTGCGTAGACGATTGCCGTACCGCCAGCGACCTGGAGATCTTCTGCCCGGGCCAGAGCCTGAAACTCATCTTTTGACCGCCCAAGGGTATATTGCCCGCCAGAGTTAACCTTTGTGGCGGAGTCCTGACCCAGGTTTTGCATTGCCCCACCCTTGCCAACCGTGGAGTCGCGAGCTTCGGTACCTGCTAATACCAATAAATTACCGCCATTTTCCAGCAACATATTGGTCGCTAAATTGCCGGAAATGGAAAAAGTACCGTACTGGTGAGTACCGCTGATTTCAATACCGTTAGCCGTGCTCGTCTGGAGAGCGGCACCGCTGTTCTGGACGATATCTGTCGCTTTGCCATTATCGTTGACTGTCAGCGTACCGCCTTCATTGATCTTTGTTTTTATTGCCTCTCCGTTAGCTGAAACTGTTTGTATTCCGCCGTCGTTAATTGTTGTCTCATTCGCCA  >  NC_000913/2779774‑2780253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑