Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,780,609 A→C G706G (GGT→GGG ypjA ← adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,780,6090AC96.9% 117.7 / ‑4.7 32G706G (GGT→GGGypjAadhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (11/20);  total (11/21)
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.98e-01

CGTCGCCCTGGTGGTTGCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGAACCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTCTCTGTTGACCCACCA  >  NC_000913/2780468‑2780753
                                                                                                                                             |                                                                                                                                                
tgTCGCCCTGGTGGTTGCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCAt                                                                                                                                           >  1:2926948/2‑149 (MQ=255)
 gTCGCCCTGGTGGTTGCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTCCCTGTGCACCCCCATtaa                                                                                                                                        <  1:2549435/151‑1 (MQ=255)
  tCGCCCTGGTGGTTGCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATtaat                                                                                                                                       >  2:2768303/1‑151 (MQ=255)
    gCCCTGGTGGTTGCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATtaata                                                                                                                                      >  2:1424594/1‑150 (MQ=255)
         ggtggtTGCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCCAGTTTTCAGCCAGGCCGACCTCATTAACTACCTGTGACCCCCCATTAATAATGGaa                                                                                                                                <  1:2696872/151‑1 (MQ=255)
                gCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTcc                                                                                                                          <  2:3737622/150‑1 (MQ=255)
                gCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCg                                                                                                                         <  1:1801447/151‑1 (MQ=255)
                      aaaGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTcc                                                                                                                   <  2:1988339/151‑1 (MQ=255)
                         gcgcCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGt                                                                                                                <  2:1498882/151‑1 (MQ=255)
                                   tACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAAt                                                                                                      >  1:2364998/1‑151 (MQ=255)
                                            gTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCAGTCCCGTTTGCCATAATTtgttgt                                                                                               <  1:2566307/149‑1 (MQ=255)
                                             tATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAgg                                                                                            >  1:2333294/1‑151 (MQ=255)
                                                            gCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCg                                                                             <  2:678278/151‑1 (MQ=255)
                                                                ttCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGt                                                                           <  2:2437025/149‑1 (MQ=255)
                                                                    cccGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTattga                                                                      <  1:1920505/150‑1 (MQ=255)
                                                                      cGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTattgattg                                                                   >  1:3130631/1‑151 (MQ=255)
                                                                           aTCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTa                                                              <  1:3555124/151‑1 (MQ=255)
                                                                             cGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTAcc                                                            <  1:691610/151‑1 (MQ=255)
                                                                                      cgccgTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTg                                                   >  2:3310951/1‑151 (MQ=255)
                                                                                         cgTCTTTAAACACCGAGTTTTAAGCCAGGCCGCCATCATTAACTACCTATGACCCCTCATTAATAATGGAACCTTATGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACAATTTCGGCATCGATAGACTTACAATAATTCTGa                                                  >  1:3404587/1‑149 (MQ=255)
                                                                                            cATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACgg                                              <  2:1233482/150‑1 (MQ=255)
                                                                                                  gCACCGAGTTTACAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGt                                       <  1:2422236/151‑1 (MQ=255)
                                                                                                            ttCAGCCAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTg                             >  2:532982/1‑151 (MQ=255)
                                                                                                                  cAGGCCGCCCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGAtgtgt                       >  2:1697620/1‑151 (MQ=255)
                                                                                                                       cgccCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGtttt                  <  1:3470258/151‑1 (MQ=255)
                                                                                                                        gccCTCATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTtc                 <  2:3356360/151‑1 (MQ=255)
                                                                                                                             cATTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTctct              <  2:2174439/149‑1 (MQ=255)
                                                                                                                              aTTAACTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTCTCTGtt           <  2:2454370/151‑1 (MQ=255)
                                                                                                                               ttAACTACCTGTGACCCCCCATTAATAATGGCACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTCTCTGtt           >  1:2822255/1‑150 (MQ=255)
                                                                                                                                  aCTACCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTCTCTGTTGAca       <  1:2028894/151‑2 (MQ=255)
                                                                                                                                     aCCTGTGACCCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTCTCTGTTGACcca     <  2:237570/150‑1 (MQ=255)
                                                                                                                                          tgACCCCCCATTAACAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTCTCTGTTGACccacca  <  1:3527911/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CGTCGCCCTGGTGGTTGCAACCAAAGCGCCCTGGCTACTCTGCTGTATCCCCGTTGCGCTGCCTTTCTCCCGCACATCGAGTGTGCCGCCGTCATTAAGCACCGAGTTTTCAGCCAGGCCGCCCTCATTAACTACCTGTGAACCCCCATTAATAATGGAACCTTCCGCTGTCCCGTTTGCCATAATTTGTTGTAGGCCAGAGACGATATCGGTATTGATTGCCTTACCATAATTCTGAACGGTTTGCGTGCCACCATTGATGTGTGTTTTCTCTGTTGACCCACCA  >  NC_000913/2780468‑2780753

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