Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,100,789 G→T 15.5% intergenic (+2/+62) ycdZ → / ← csgG DUF1097 family inner membrane protein/curli production assembly/transport outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,100,7890GT15.5% 137.0 / 14.2 45intergenic (+2/+62)ycdZ/csgGDUF1097 family inner membrane protein/curli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base T (4/3);  total (24/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.92e-01

CTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACC  >  NC_000913/1100655‑1100937
                                                                                                                                      |                                                                                                                                                    
cTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCttt                                                                                                                                        <  3:369417/149‑1 (MQ=255)
 tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCt                                                                                                                                          <  4:1059507/146‑1 (MQ=255)
 tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTc                                                                                                                                       >  3:68453/1‑149 (MQ=255)
 tGGAAACTGGTGTTACCTTCGCTGGAGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCt                                                                                                                                          >  3:1059507/1‑146 (MQ=255)
               aCCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACa                                                                                                                              >  1:1144356/1‑144 (MQ=255)
               aCCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACa                                                                                                                              <  2:1144356/144‑1 (MQ=255)
                  ttCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAGTg                                                                                                                      >  3:1393531/1‑149 (MQ=255)
                             ttGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATaa                                                                                                           >  2:883020/1‑149 (MQ=255)
                                ggCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAGTGCCGGATGATAATTc                                                                                                        <  4:1393531/149‑1 (MQ=255)
                                             ttACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgat                                                                                                                >  3:543715/1‑128 (MQ=255)
                                             ttACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgat                                                                                                                <  4:543715/128‑1 (MQ=255)
                                               aCGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGga                                                                                                                 <  3:771856/125‑1 (MQ=255)
                                               aCGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGga                                                                                                                 >  4:771856/1‑125 (MQ=255)
                                                 gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCttt                                                                                                 >  4:35296/1‑139 (MQ=255)
                                                 gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCttt                                                                                                 <  3:35296/139‑1 (MQ=255)
                                                           aCAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGcc                                                                                                                    >  2:331073/1‑110 (MQ=255)
                                                           aCAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGcc                                                                                                                    <  1:331073/110‑1 (MQ=255)
                                                                    tgtgGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGa                                                                    <  1:883020/149‑1 (MQ=255)
                                                                                   gTAGTGCAAAGACGGCGCACCTTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATc                                                                                           <  3:1356710/111‑1 (MQ=255)
                                                                                   gTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATc                                                                                           >  4:1356710/1‑111 (MQ=255)
                                                                                       tGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCttt                                                                                                                                        >  4:1266169/1‑62 (MQ=255)
                                                                                       tGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCttt                                                                                                                                        <  3:1266169/62‑1 (MQ=255)
                                                                                                          gAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCttt                                                                                                 <  1:647189/82‑1 (MQ=255)
                                                                                                          gAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCttt                                                                                                 >  2:647189/1‑82 (MQ=255)
                                                                                                             cAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTaaaa                                                                                                                        <  4:530265/56‑1 (MQ=255)
                                                                                                             cAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTaaaa                                                                                                                        >  3:530265/1‑56 (MQ=255)
                                                                                                              aGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGc                          >  4:581184/1‑149 (MQ=255)
                                                                                                                  aaTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCtttt                                                                                                >  3:331054/1‑75 (MQ=255)
                                                                                                                  aaTCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCtttt                                                                                                <  4:331054/75‑1 (MQ=255)
                                                                                                                      aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGaa                                             <  1:109565/122‑1 (MQ=255)
                                                                                                                      aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGaa                                             >  2:109565/1‑122 (MQ=255)
                                                                                                                      aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTg                                   <  3:1051806/132‑1 (MQ=255)
                                                                                                                      aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTg                                   >  4:1051806/1‑132 (MQ=255)
                                                                                                                       aaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCg                                                                              <  3:561466/88‑1 (MQ=255)
                                                                                                                       aaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCg                                                                              >  4:561466/1‑88 (MQ=255)
                                                                                                                           tAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCaaa             >  3:466822/1‑149 (MQ=255)
                                                                                                                           gaaaGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTc                                                                                <  4:1228510/81‑1 (MQ=38)
                                                                                                                           gaaaGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTc                                                                                >  3:1228510/2‑82 (MQ=38)
                                                                                                                            aaaGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCa                                                   >  4:420854/1‑110 (MQ=255)
                                                                                                                            aaaGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCa                                                   <  3:420854/110‑1 (MQ=255)
                                                                                                                              aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCt                         <  1:525144/134‑1 (MQ=255)
                                                                                                                              aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCAGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCt                         >  2:525144/1‑134 (MQ=255)
                                                                                                                               gCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc         >  2:561331/1‑149 (MQ=255)
                                                                                                                                    aGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTc                                                    >  3:872875/1‑101 (MQ=255)
                                                                                                                                    aGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTc                                                    <  4:872875/101‑1 (MQ=255)
                                                                                                                                      gggCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGAcc  >  4:487748/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                                    
CTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACC  >  NC_000913/1100655‑1100937

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