Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,962,441 C→T *108* (TAG→TAA ppdC ← putative prepilin peptidase‑dependent protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,4410CT96.7% 180.1 / ‑1.8 60*108* (TAG→TAAppdCputative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (32/26);  total (33/27)
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.70e-01

ACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCG  >  NC_000913/2962312‑2962574
                                                                                                                                 |                                                                                                                                     
aCTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGa                                                                                                                              >  1:1396170/1‑139 (MQ=255)
   ttGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCGACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTc                                                                                                                           >  2:1067341/1‑139 (MQ=255)
    tGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTg                                                                                                                                  >  2:863711/1‑131 (MQ=255)
    tGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTg                                                                                                                                  <  1:863711/131‑1 (MQ=255)
      cACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                        >  2:187286/1‑139 (MQ=255)
              tcatcTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAg                                                                                                                >  1:1166476/1‑139 (MQ=255)
                  ctctGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCgg                                                                                                            <  2:145396/139‑1 (MQ=255)
                        ttCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCAtc                                                                                                      >  2:587409/1‑139 (MQ=255)
                            aaaGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCTTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGc                                                                                                  <  2:1602240/139‑1 (MQ=255)
                                 aTCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCgg                                                                                                            >  1:279717/1‑124 (MQ=255)
                                 aTCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCgg                                                                                                            <  2:279717/124‑1 (MQ=255)
                                            gttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCggg                                                                                  <  2:1166476/139‑1 (MQ=255)
                                            gttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCggg                                                                                  <  2:1396170/139‑1 (MQ=255)
                                               cgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGa                                                                               <  1:587409/139‑1 (MQ=255)
                                                       aCAATAAACTCCATTAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCgg                                                                                                            <  1:45561/102‑1 (MQ=255)
                                                       aCAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCgg                                                                                                            >  2:45561/1‑102 (MQ=255)
                                                              aCTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCt                                                                <  1:187286/139‑1 (MQ=255)
                                                              aCTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCt                                                                >  1:49964/1‑139 (MQ=255)
                                                                tCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGa                                                              >  1:1237853/1‑139 (MQ=255)
                                                                  cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCt                                                                >  1:58092/1‑135 (MQ=255)
                                                                  cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCt                                                                <  2:58092/135‑1 (MQ=255)
                                                                  cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                            >  2:1190936/1‑139 (MQ=255)
                                                                    tCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCt                                                          >  1:17199/1‑139 (MQ=255)
                                                                            tCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATc                                                  >  2:1688452/1‑139 (MQ=255)
                                                                                     tCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         <  2:1237853/139‑1 (MQ=255)
                                                                                     tCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         >  1:1363038/1‑139 (MQ=255)
                                                                                     tCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         >  2:483660/1‑139 (MQ=255)
                                                                                      ccAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGa                                        <  1:1067341/139‑1 (MQ=255)
                                                                                       cAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCa                                                                                                                 >  1:691629/1‑65 (MQ=255)
                                                                                       cAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCa                                                                                                                 <  2:691629/65‑1 (MQ=255)
                                                                                              ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCa                                <  2:1363038/139‑1 (MQ=255)
                                                                                                    tGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGt                          <  2:17199/139‑1 (MQ=255)
                                                                                                      gTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 >  1:558785/1‑82 (MQ=255)
                                                                                                      gTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 <  2:558785/82‑1 (MQ=255)
                                                                                                      gTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 <  2:389395/82‑1 (MQ=255)
                                                                                                      gTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 >  1:389395/1‑82 (MQ=255)
                                                                                                         gACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         <  2:401874/119‑1 (MQ=255)
                                                                                                         gACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         >  1:401874/1‑119 (MQ=255)
                                                                                                             cTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                        <  1:73060/100‑1 (MQ=255)
                                                                                                             cTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                        >  2:73060/1‑100 (MQ=255)
                                                                                                                aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGt              >  1:425099/1‑139 (MQ=255)
                                                                                                                aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGt              >  1:20517/1‑139 (MQ=255)
                                                                                                                     aGCGGCTCCTGCTTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTGCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGGCTGCATTCGGTTGACCTGCCAGTTGGCa         >  2:980633/1‑139 (MQ=255)
                                                                                                                     aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         <  2:1535117/107‑1 (MQ=255)
                                                                                                                     aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCg                                         >  1:1535117/1‑107 (MQ=255)
                                                                                                                     aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGt                          <  1:1381214/122‑1 (MQ=255)
                                                                                                                     aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGt                          >  2:1381214/1‑122 (MQ=255)
                                                                                                                     aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATcc                                                 <  1:548192/98‑1 (MQ=255)
                                                                                                                     aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATcc                                                 >  2:548192/1‑98 (MQ=255)
                                                                                                                      gCGGCTCCTTATTACTGACGATTCGGGCAATGCAgg                                                                                                               >  1:1137993/1‑36 (MQ=255)
                                                                                                                      gCGGCTCCTGATTACTGACGATTCGGGCAATGCAgg                                                                                                               <  2:1137993/36‑1 (MQ=255)
                                                                                                                       cGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAgg       <  2:49964/139‑1 (MQ=255)
                                                                                                                        ggCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcc                                                                                        <  1:1107730/57‑1 (MQ=255)
                                                                                                                        ggCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcc                                                                                        >  2:1107730/1‑57 (MQ=255)
                                                                                                                          cTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                                                                   <  2:1330877/60‑1 (MQ=255)
                                                                                                                          cTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                                                                   >  1:1330877/1‑60 (MQ=255)
                                                                                                                            ccTGATTACTGACGATTCGGGCAATGCAGGCGGGt                                                                                                          >  1:1402595/1‑35 (MQ=255)
                                                                                                                            ccTGATTACTGACGATTCGGGCAATGCAGGCGGGt                                                                                                          <  2:1402595/35‑1 (MQ=255)
                                                                                                                            ccTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg  >  2:305288/1‑139 (MQ=255)
                                                                                                                            ccTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg  >  2:1400320/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                     
ACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCG  >  NC_000913/2962312‑2962574

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: