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,4410CT98.0% 155.0 / NA 49*108* (TAG→TAAppdCputative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/23);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCG  >  NC_000913/2962310‑2962574
                                                                                                                                   |                                                                                                                                     
gTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGAc                                                                                                                                >  1:825887/1‑139 (MQ=255)
gTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGAc                                                                                                                                >  1:344223/1‑139 (MQ=255)
   cTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGAt                                                                                                                             <  2:344223/139‑1 (MQ=255)
        cACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCggg                                                                                                                        <  1:1010548/139‑1 (MQ=255)
               atcatcTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCa                                                                                                                 <  2:825887/139‑1 (MQ=255)
                     tctGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCggg                                                                                                           >  2:962058/1‑139 (MQ=255)
                        ggTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCa                                                                                                        >  2:166276/1‑139 (MQ=255)
                         gTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCAt                                                                                                       <  1:962058/139‑1 (MQ=255)
                          ttCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCAtc                                                                                                      <  2:83665/139‑1 (MQ=255)
                                    tCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGAt                                                                                                                             >  2:678309/1‑106 (MQ=255)
                                    tCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGAt                                                                                                                             <  1:678309/106‑1 (MQ=255)
                                          aGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                      <  2:1128712/139‑1 (MQ=255)
                                            ccGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCg                                                                                    >  2:758441/1‑139 (MQ=255)
                                               ttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 >  2:213173/1‑139 (MQ=255)
                                                 cgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGa                                                                               >  1:934801/1‑139 (MQ=255)
                                                    gttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAAc                                                                            <  1:706657/139‑1 (MQ=255)
                                                                aCTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCt                                                                <  1:166276/139‑1 (MQ=255)
                                                                  tCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAgg                                                                                                               <  2:953880/90‑1 (MQ=255)
                                                                  tCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAgg                                                                                                               >  1:953880/1‑90 (MQ=255)
                                                                  tCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGa                                                              >  1:321692/1‑139 (MQ=255)
                                                                    cATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTTGCGTAACGCTGATg                                                            >  2:190801/1‑139 (MQ=255)
                                                                            cTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                      <  2:959483/105‑1 (MQ=255)
                                                                            cTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcccc                                                                                      >  1:959483/1‑105 (MQ=255)
                                                                            cTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGAcaca                                                    >  1:658026/1‑139 (MQ=255)
                                                                                   cACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTTGCGTAACGCTGATGCTGACACATCCCGcc                                             <  1:190801/139‑1 (MQ=255)
                                                                                     cGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTg                                           >  1:1122944/1‑139 (MQ=255)
                                                                                            aCGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                        >  2:244795/1‑119 (MQ=255)
                                                                                            aCGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGa                                                        <  1:244795/119‑1 (MQ=255)
                                                                                               aTTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCggg                                                                                  >  1:338149/1‑90 (MQ=255)
                                                                                               aTTGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCggg                                                                                  <  2:338149/90‑1 (MQ=255)
                                                                                                ttGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCa                                >  1:301887/1‑139 (MQ=255)
                                                                                                 tGGAATGGTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAt                               >  1:239102/1‑139 (MQ=255)
                                                                                                       ggTAGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAAATAGCGTAACACTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGtt                         <  2:698256/139‑1 (MQ=255)
                                                                                                          aGACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGAc                      <  2:1122944/139‑1 (MQ=255)
                                                                                                           gACCCTTAACATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGAcc                     <  2:239102/139‑1 (MQ=255)
                                                                                                                  aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATcc                                                 <  1:477110/104‑1 (MQ=255)
                                                                                                                  aaCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATcc                                                 >  2:477110/1‑104 (MQ=255)
                                                                                                                  aaCATAGCGGCTCCTGATTACTGACGATGCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGt              >  1:702589/1‑139 (MQ=255)
                                                                                                                   aCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCg                                                                                                             >  2:91919/1‑43 (MQ=255)
                                                                                                                   aCATAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCg                                                                                                             <  1:91919/43‑1 (MQ=255)
                                                                                                                     aTAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATg                                                            >  2:771068/1‑90 (MQ=255)
                                                                                                                     aTAGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTCAAACTAGCGTAACGCTTATg                                                            <  1:771068/90‑1 (MQ=255)
                                                                                                                       aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 >  2:623676/1‑67 (MQ=255)
                                                                                                                       aGCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                                                                 <  1:623676/67‑1 (MQ=255)
                                                                                                                        gCGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCAGCCAGTTGGCAg        >  2:1080438/1‑139 (MQ=255)
                                                                                                                         cGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAtt                              <  2:151446/116‑1 (MQ=255)
                                                                                                                         cGGCTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAtt                              >  1:151446/1‑116 (MQ=255)
                                                                                                                            cTCCTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCCGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg    >  1:450531/1‑139 (MQ=255)
                                                                                                                              ccTGATTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCTTGCGACGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg  <  1:758441/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCG  >  NC_000913/2962310‑2962574

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