Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,097,764 A→G N78N (AAT→AAC csgE ← curli production assembly/transport protein CsgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,097,7640AG96.5% 93.5 / ‑3.2 29N78N (AAT→AACcsgEcurli production assembly/transport protein CsgE
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (7/21);  total (7/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

AATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTAT  >  NZ_CP009273/1097628‑1097897
                                                                                                                                        |                                                                                                                                     
aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGtt                                                                                                                                     <  1:487404/139‑1 (MQ=255)
  tCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGa                                                                                                                                   <  2:631600/139‑1 (MQ=255)
               ataaCGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGAt                                                                                                                      <  1:119494/139‑1 (MQ=255)
                   cGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGtt                                                                                                                                     <  2:788441/120‑1 (MQ=255)
                   cGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGtt                                                                                                                                     >  1:788441/1‑120 (MQ=255)
                   cGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAg                                                                                                                  <  2:23473/139‑1 (MQ=255)
                    gCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGc                                                                                                                 >  2:297098/1‑139 (MQ=255)
                         aTTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGTCCGTTATAGTGATCCAGCTTcc                                                                                                             <  2:438967/138‑1 (MQ=255)
                         aTTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTcc                                                                                                             <  1:268149/138‑1 (MQ=255)
                            tATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCa                                                                                                          <  1:825888/138‑1 (MQ=255)
                                  gCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTg                                                                                                     <  1:297098/137‑1 (MQ=255)
                                        aTTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGc                                                                                             <  2:20733/139‑1 (MQ=255)
                                                   cttcAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAAt                                                                                  >  2:784867/1‑139 (MQ=255)
                                                   cttcAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAAt                                                                                  >  2:195779/1‑139 (MQ=255)
                                                             aaTCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTa                                                                        <  1:46640/139‑1 (MQ=255)
                                                               tCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTAcc                                                                      <  2:387102/139‑1 (MQ=255)
                                                                       aaGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGa                                                                                                                                   <  2:806233/70‑1 (MQ=255)
                                                                       aaGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGa                                                                                                                                   >  1:806233/1‑70 (MQ=255)
                                                                       aaGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGt                                                              <  2:178723/139‑1 (MQ=255)
                                                                         gacgacAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCa                                                            <  1:733866/139‑1 (MQ=255)
                                                                               aGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGt                                                                                                                               <  2:118606/66‑1 (MQ=255)
                                                                               aGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGt                                                                                                                               >  1:118606/1‑66 (MQ=255)
                                                                                        gAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTa                                                <  2:15952/136‑1 (MQ=255)
                                                                                           gTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTa                                          <  1:550584/139‑1 (MQ=255)
                                                                                                   tCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTc                                   <  1:664815/138‑1 (MQ=255)
                                                                                                                       tGGAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCaa                 >  1:7235/1‑136 (MQ=255)
                                                                                                                         gAAAATAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGa             <  1:456978/138‑1 (MQ=255)
                                                                                                                           aaaTAACGTCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGa          <  1:233548/139‑1 (MQ=255)
                                                                                                                                   tCCTGGTTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTAt  <  2:19179/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
AATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTAT  >  NZ_CP009273/1097628‑1097897

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