Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,215,879 G→A 9.8% G31R (GGG→AGG)  ycgZ → RcsB connector protein for regulation of biofilm and acid‑resistance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,215,8790GA 9.8% 103.6 / 4.1 41G31R (GGG→AGG) ycgZRcsB connector protein for regulation of biofilm and acid‑resistance
Reads supporting (aligned to +/- strand):  ref base G (23/14);  new base A (2/2);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.81e-01

AAGTTTGCTATTGTTACTTCACTTAACATTGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACATTATAATGCACTAATTGGGTTGCTTTTTGAATAATG  >  NC_000913/1215732‑1216027
                                                                                                                                                   |                                                                                                                                                    
aaGTTTGCTATTTTTACTTCACTTAACATTGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCgg                                                                                                                                                     <  1:962400/149‑1 (MQ=255)
           tGTTCCTTCACTTAACATGGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAACAGCTAACTTGCCTATTCAGCAGGAAACTCTCGGGGTCATtgtta                                                                                                                                          <  2:634689/149‑3 (MQ=255)
           tGTTACTTCACTTAACATTGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATtgtgg                                                                                                                                          <  2:948842/149‑2 (MQ=255)
           tGTTAATTGCCTAAAAATTGATTAACATTTTTAACAGGCACGTAGCATGCATGAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGa                                                                                                                                          <  2:234513/149‑1 (MQ=255)
            gTTACTTCACTTAACATTGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGAc                                                                                                                                         >  1:1248965/1‑149 (MQ=255)
                  tCACTTAACATTGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTACCTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGAt                                                                                                                                   >  2:235995/1‑149 (MQ=255)
                   cACTTAACATTGACTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCCGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGAtt                                                                                                                                  <  1:1248799/149‑1 (MQ=255)
                           aTTGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGa                                                                                                                          >  1:858663/1‑149 (MQ=255)
                              gATTCACATTTTTAACAGATGCGTAGCATGCATCAATATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTCGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATgg                                                                                                                       <  2:727153/149‑1 (MQ=255)
                                  aaCATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGCTTTTGAAAGATGGACGt                                                                                                                   >  1:986040/1‑149 (MQ=255)
                                     aTTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTGACTTTGCAAAAGCTAACTTGCATACTCAGCAGGCAACTCGCGGGGAAATAGTGACTGAGATTTTGAAAGATGGACGTAAt                                                                                                                >  1:715073/1‑149 (MQ=255)
                                      tttttAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCAGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATc                                                                                                               >  1:125769/1‑149 (MQ=255)
                                      tttttAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCAGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATc                                                                                                               >  1:1136469/1‑149 (MQ=255)
                                         ttAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGCCTGAGATTTTGAAAGATGAACGTAATCTGa                                                                                                            >  2:863755/1‑149 (MQ=255)
                                               gagGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGAGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGCAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGaa                                                                                                      >  1:375728/1‑149 (MQ=255)
                                                     tAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGc                                                                                                >  1:206309/1‑149 (MQ=255)
                                                     tAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAAGCAGAGTCGAAAGTCGc                                                                                                >  1:457135/1‑149 (MQ=255)
                                                     tAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACCCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGCCGTAATCTGGGTCGAAAGTCGc                                                                                                >  1:839771/1‑149 (MQ=255)
                                                      aGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCt                                                                                               <  2:875042/149‑1 (MQ=255)
                                                        catgcatCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCttt                                                                                             >  2:953118/1‑149 (MQ=255)
                                                             atCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGGCGAAAGTCGCTTTGTGcc                                                                                        >  1:1055507/1‑149 (MQ=255)
                                                                              ttAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGCCt                                                                       >  1:418735/1‑149 (MQ=255)
                                                                                      tGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGCACAcg                                                               >  2:774608/1‑149 (MQ=255)
                                                                                      tGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAaaacc                                                               >  2:774662/1‑148 (MQ=255)
                                                                                      tGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACAcg                                                               >  2:40000/1‑149 (MQ=255)
                                                                                             gCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCAGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTgg                                                        <  2:1136469/149‑1 (MQ=255)
                                                                                             gCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCAGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTgg                                                        <  2:125769/149‑1 (MQ=255)
                                                                                                aTCACTCGTTACTTTGCAAAAGCTAACTTGCATAATCAACAGGAAACTCTCGGGGAAATAATGACTGAGAATTTGAAAGATGGACGGAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTggggc                                                     >  2:1097097/1‑148 (MQ=255)
                                                                                                   ccacGTTACGTTGCAAAAGCTCACTTCCATACTCCGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGCACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAaga                                                  <  2:611389/146‑1 (MQ=255)
                                                                                                     tCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAgg                                                >  1:1085325/1‑149 (MQ=255)
                                                                                                               gCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGTGACTGGGGAAGAGGAACAGAAACa                                      >  1:877650/1‑149 (MQ=255)
                                                                                                               gCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGGGGAACAGAAACa                                      >  2:547176/1‑149 (MQ=255)
                                                                                                                cAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGGGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGGACAGAAACag                                     >  1:122321/1‑148 (MQ=255)
                                                                                                                 aaaaGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACAtt                                    <  1:1263025/149‑1 (MQ=255)
                                                                                                                     gCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACATTATaa                                <  1:1160128/149‑1 (MQ=255)
                                                                                                                       tAACTGGTATTCTTAGCAGGAAACTCTCGGGGAAATTGTGACCGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACATTATAATg                              <  1:235995/149‑1 (MQ=255)
                                                                                                                         aCTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACATTATAATGCa                            <  1:422564/149‑1 (MQ=255)
                                                                                                                         aCTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACAGGGGAAGAGGAACAGAAACATTATCATGCa                            >  1:106958/1‑149 (MQ=255)
                                                                                                                                            aaCTCTCGGGGAAATTGTTACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACATTATAATGCACTAATTGGGTTGCTTTTTg         <  2:740767/149‑1 (MQ=255)
                                                                                                                                             agtgtcGGGGAATTTGTGCCTGGGATTTTGAAAGATGGACGTAATCAGCATATAGACGCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGATTGGGGAAGAGGAACAGAAACATTATAATGCACTAATTGGGTTGCTTTTTGa        <  1:1213248/145‑1 (MQ=255)
                                                                                                                                               tctcGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGAGACTGGGGAAGAGGAACAGAAACATTATAATGCACTAATTGGGTTGCTTTTTGaat      >  2:163239/1‑149 (MQ=255)
                                                                                                                                                ctcGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACATTATAATGCACTAATTGGGTTGCTTTTTGaata     <  1:239332/149‑1 (MQ=255)
                                                                                                                                                   ggggAAATTGTGACTGAGATTTTGAAAGACGGACGTNACCTGAGTCGAAAGTCGCTTTGTGCCAAGCTTCTTTGCCGACTGGAACTCGCGACTGGGGAAGAGGAACAGAATCATTATAATGCACTAATTGGGTTGCTTTTTGAACAATg  <  2:106958/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGTTTGCTATTGTTACTTCACTTAACATTGATTAACATTTTTAACAGAGGCGTAGCATGCATCAAAATTCAGTGACTTTAGATTCTGCGGGAGCAATCACTCGTTACTTTGCAAAAGCTAACTTGCATACTCAGCAGGAAACTCTCGGGGAAATTGTGACTGAGATTTTGAAAGATGGACGTAATCTGAGTCGAAAGTCGCTTTGTGCCAAACTTCTTTGCCGACTGGAACACGCGACTGGGGAAGAGGAACAGAAACATTATAATGCACTAATTGGGTTGCTTTTTGAATAATG  >  NC_000913/1215732‑1216027

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