Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,048,572 G→A G40G (GGC→GGT gfcC ← capsule biosynthesis GfcC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,5720GA98.0% 159.4 / NA 49G40G (GGC→GGTgfcCcapsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/36);  minor base C (1/0);  total (13/36)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGT  >  NC_000913/1048436‑1048694
                                                                                                                                        |                                                                                                                          
caagcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCa                                                                                                                         <  2:337189/140‑1 (MQ=255)
      gtgCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGaa                                                                                                                   >  1:418782/1‑140 (MQ=255)
       tgCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGaaa                                                                                                                  <  2:600771/140‑1 (MQ=255)
       tgCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGaaa                                                                                                                  <  2:120515/140‑1 (MQ=255)
          cATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAgtg                                                                                                               >  1:665224/1‑140 (MQ=255)
             gACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCt                                                                                                            >  1:412145/1‑140 (MQ=255)
                aTGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCtgtt                                                                                                         <  1:521790/140‑1 (MQ=255)
                 tGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCtgttg                                                                                                        >  2:97926/1‑140 (MQ=255)
                 tGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCtgttg                                                                                                        >  2:132510/1‑140 (MQ=255)
                      gTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGc                                                                                                   <  2:418782/140‑1 (MQ=255)
                      gTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGc                                                                                                   <  1:124999/140‑1 (MQ=255)
                          tCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAgg                                                                                               <  1:28414/140‑1 (MQ=255)
                          tCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAgg                                                                                               <  2:14117/140‑1 (MQ=255)
                          tCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAgg                                                                                               <  1:660847/140‑1 (MQ=255)
                          tCTTTAAGCGCTTTAGCCTTCGCCGCAGAGTCGGTAAGCCAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAgg                                                                                               <  2:205254/140‑1 (MQ=255)
                            tttAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCa                                                                                             >  1:643135/1‑140 (MQ=255)
                            tttAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACAGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCa                                                                                             >  2:551708/1‑140 (MQ=255)
                              tAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGa                                                                                           >  1:621915/1‑140 (MQ=255)
                                   gcTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACAGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAAt                                                                                      <  1:551708/140‑1 (MQ=255)
                                    cTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATa                                                                                     >  1:4856/1‑140 (MQ=255)
                                    cTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTATCCACGGGCCCAGCGGGAAGTGTCTGTTGTGCGCCAGGCAGATAAATa                                                                                     >  2:189647/1‑140 (MQ=255)
                                     tttAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAg                                                                                    <  1:593366/140‑1 (MQ=255)
                                               gccgcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTcc                                                                          <  2:131227/140‑1 (MQ=255)
                                                  gcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTg                                                                       <  2:371767/140‑1 (MQ=255)
                                                    agagTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTgcg                                                                     >  1:351737/1‑140 (MQ=255)
                                                      agTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                                                   <  1:632803/140‑1 (MQ=255)
                                                      agTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                                                   <  2:412145/140‑1 (MQ=255)
                                                      agTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                                                   <  2:643135/140‑1 (MQ=255)
                                                      agTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                                                   <  2:665221/140‑1 (MQ=255)
                                                      agTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCACGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                                                   <  1:294618/140‑1 (MQ=255)
                                                      agTCGGTCAGCAAGGCACCAGGACACCAGAGTCGACCGCGCAGTTGCGCTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                                                   <  2:179792/140‑1 (MQ=255)
                                                            tCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCa                                                             <  1:114043/140‑1 (MQ=255)
                                                              aGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATg                                                           <  2:601270/140‑1 (MQ=255)
                                                                     cACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTc                                                    <  1:512603/140‑1 (MQ=255)
                                                                     cACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTc                                                    <  2:471280/140‑1 (MQ=255)
                                                                      aCCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCa                                                   <  2:563015/140‑1 (MQ=255)
                                                                      aCCAGGCCACCAGAGGCCATCGTGCAGTTGCAGTTCTGTTACCAACTGAACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCa                                                   <  2:200765/140‑1 (MQ=255)
                                                                        cAGGCCACCAGAGGCGATCGCGTAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCAtt                                                 <  1:267898/140‑1 (MQ=255)
                                                                          ggCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTAc                                               <  1:133721/140‑1 (MQ=255)
                                                                                     gCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGc                                    >  2:113695/1‑140 (MQ=255)
                                                                                        aTCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTgg                                 <  1:132510/140‑1 (MQ=255)
                                                                                           gcgcAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCaa                              <  2:432893/140‑1 (MQ=255)
                                                                                                  tGCGGTTGTGTCACCTACTGGACAATATTCTCCACGGTACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTAATTACGTAAAGTACGCTGGCAATGAAATa                       <  1:189647/140‑1 (MQ=255)
                                                                                                    cGGTTGTGTCACCAACTAGACAGCATTTTCCTCGGAACCTACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACg                     <  2:301095/140‑1 (MQ=255)
                                                                                                     ggTTGTGTCACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGa                    <  1:97926/140‑1 (MQ=255)
                                                                                                       ttGTGTCACCAACTGGACAACATTTCCCAAGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACt                  <  1:269170/140‑1 (MQ=255)
                                                                                                             cACCAACTGGACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAAtt            <  1:113695/140‑1 (MQ=255)
                                                                                                                      gACAACATTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGCCATGGGGCGTCGTTACGTAATGTATGCTGGCAATGAAGTACGGCTGTAATTTATTCATgg   >  2:175547/1‑140 (MQ=255)
                                                                                                                       acaacaTTTTCCACGGGACCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGt  <  2:665224/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                          
CAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGT  >  NC_000913/1048436‑1048694

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