Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,102,221 T→C E107E (GAA→GAG csgE ← curlin secretion specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,102,2210TC98.2% 195.3 / NA 55E107E (GAA→GAGcsgEcurlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (27/27);  minor base A (0/1);  total (27/28)
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

GAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCC  >  NC_000913/1102087‑1102336
                                                                                                                                      |                                                                                                                   
gAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                 <  1:307843/139‑1 (MQ=255)
 aaaTAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCctcttc                                                                                                                >  1:172971/1‑139 (MQ=255)
          tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                       <  2:11375/139‑1 (MQ=255)
          tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                       <  2:960361/139‑1 (MQ=255)
          tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                       <  2:172971/139‑1 (MQ=255)
                      ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTGTTCAGTTTGAATCAGTGCAAAgac                                                                                           >  1:798857/1‑139 (MQ=255)
                           tGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAg                                                                                      >  1:37881/1‑139 (MQ=255)
                                    cGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGGCGACAGTTTTCTCGa                                                                             <  2:955708/139‑1 (MQ=255)
                                     gCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGaa                                                                            >  2:493028/1‑139 (MQ=255)
                                       aTTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                       <  2:577129/110‑1 (MQ=255)
                                       aTTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                       >  1:577129/1‑110 (MQ=255)
                                        ttGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCg                                                                              >  1:299912/1‑134 (MQ=255)
                                        ttGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCg                                                                              <  2:299912/134‑1 (MQ=255)
                                              ccTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                       <  1:234666/103‑1 (MQ=255)
                                              ccTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                       >  2:234666/1‑103 (MQ=255)
                                                 ccGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCa                                                                                                     >  2:693780/1‑102 (MQ=255)
                                                 ccGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCa                                                                                                     <  1:693780/102‑1 (MQ=255)
                                                 ccGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCaa                                                                >  2:526772/1‑139 (MQ=255)
                                                  cGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAAc                                                               <  2:37881/139‑1 (MQ=255)
                                                            tttatttaGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                             <  1:1398137/115‑1 (MQ=255)
                                                            tttatttaGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                             >  2:1398137/1‑115 (MQ=255)
                                                                          atcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATa                                       >  2:753655/1‑139 (MQ=255)
                                                                               gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGtct                                                                        <  2:835410/101‑1 (MQ=255)
                                                                               gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGtct                                                                        >  1:835410/1‑101 (MQ=255)
                                                                                     aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            <  1:526772/139‑1 (MQ=255)
                                                                                      aTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGt                                                                                                             >  2:165687/1‑57 (MQ=255)
                                                                                      aTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCATCTTCAGt                                                                                                             <  1:165687/57‑1 (MQ=255)
                                                                                       tCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGa                          <  2:798857/139‑1 (MQ=255)
                                                                                        cGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGAc                         <  1:454113/139‑1 (MQ=255)
                                                                                          cccGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGt                                                 >  2:424766/1‑113 (MQ=255)
                                                                                          cccGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGt                                                 <  1:424766/113‑1 (MQ=255)
                                                                                                    ataaCGCCTGATTTATCTGGCGACGAATTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGAt             <  1:493028/139‑1 (MQ=255)
                                                                                                        cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCaaa                                                                                              <  2:535434/54‑1 (MQ=255)
                                                                                                        cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCaaa                                                                                              >  1:535434/1‑54 (MQ=255)
                                                                                                        cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                          <  2:395771/74‑1 (MQ=255)
                                                                                                        cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                          >  1:395771/1‑74 (MQ=255)
                                                                                                        cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAg         >  2:240526/1‑139 (MQ=255)
                                                                                                         gCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGc        >  1:304511/1‑139 (MQ=255)
                                                                                                         gCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGc        >  1:963966/1‑139 (MQ=255)
                                                                                                          ccTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGAt             <  2:128319/133‑1 (MQ=255)
                                                                                                          ccTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGAt             >  1:128319/1‑133 (MQ=255)
                                                                                                           cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            >  2:583045/1‑117 (MQ=255)
                                                                                                           cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            <  1:583045/117‑1 (MQ=255)
                                                                                                             gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc    <  1:240526/139‑1 (MQ=255)
                                                                                                             gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc    <  2:304511/139‑1 (MQ=255)
                                                                                                             gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc    >  1:4431/1‑139 (MQ=255)
                                                                                                               tttATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGCAGTCTCTTTTCAACGGAAATAAAAAAGTCTGTAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcccc  >  2:931674/1‑139 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                             <  2:469139/55‑1 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                             >  1:469139/1‑55 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            <  2:419568/104‑1 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            <  2:1334331/104‑1 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            <  2:63972/104‑1 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            >  1:63972/1‑104 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            >  1:419568/1‑104 (MQ=255)
                                                                                                                        cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                            >  1:1334331/1‑104 (MQ=255)
                                                                                                                             aTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAg         <  2:533374/118‑1 (MQ=255)
                                                                                                                             aTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAg         >  1:533374/1‑118 (MQ=255)
                                                                                                                                      |                                                                                                                   
GAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCC  >  NC_000913/1102087‑1102336

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