Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,457,425 T→C V350A (GTT→GCT)  gspD → Type II secretion system protein GspD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,457,4250TC98.0% 175.5 / NA 50V350A (GTT→GCT) gspDType II secretion system protein GspD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (26/23);  minor base A (1/0);  total (27/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

CGATGGATAACGTCGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGAT  >  NC_000913/3457296‑3457561
                                                                                                                                 |                                                                                                                                        
cgATGGATAACGTCGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAatc                                                                                                                                >  2:437608/1‑140 (MQ=255)
  agggATAACGTCGCAGTTACCGCCGCTGACCAGACTATCTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAatcat                                                                                                                              <  1:204760/138‑1 (MQ=255)
       aacggccGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTg                                                                                                                         <  2:462016/134‑1 (MQ=255)
             cGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGTACATTCGCGGTGAACAGGTGCTGGCAGAGGAAATCATCGTTGAAGTTc                                                                                                                   >  2:274693/1‑140 (MQ=255)
             cGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTc                                                                                                                   >  2:494620/1‑140 (MQ=255)
             cGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTc                                                                                                                   >  1:385919/1‑140 (MQ=255)
             cGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGAGCACAGGGGCTGGCTGAGGCAATCATCGTTGAAGTTc                                                                                                                   >  1:229052/1‑140 (MQ=255)
              gCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCAACGTTGAAGTTCa                                                                                                                  >  2:533484/1‑140 (MQ=255)
               ccATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGGACAGGTGATGGATGAGGCAATCATCGTTGAAGTTCAg                                                                                                                 >  1:262254/1‑140 (MQ=255)
                cATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAgg                                                                                                                >  1:99252/1‑140 (MQ=255)
                 aTTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGa                                                                                                               <  1:496839/140‑1 (MQ=255)
                           aTGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGa                                                                                                     <  2:608255/140‑1 (MQ=255)
                              aaCAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTa                                                                                                  >  2:158220/1‑140 (MQ=255)
                                   aCTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCt                                                                                             >  2:606245/1‑140 (MQ=255)
                                   aCTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCt                                                                                             >  1:526615/1‑140 (MQ=255)
                                        ctctctGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCg                                                                                           <  1:500784/137‑1 (MQ=255)
                                        ctctctGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCg                                                                                           >  2:500784/1‑137 (MQ=255)
                                            ctGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCa                                                                                    <  1:384781/140‑1 (MQ=255)
                                             tGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCaa                                                                                   >  2:727367/1‑140 (MQ=255)
                                             tGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCaa                                                                                   >  2:488615/1‑140 (MQ=255)
                                             tGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCAACGTTGAAGTTCAGGATGGAAATGGACTAAACCTCAGCGTGCaa                                                                                   <  1:533484/140‑1 (MQ=255)
                                                  aTTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCACCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGc                                                                              <  1:130311/140‑1 (MQ=255)
                                                   ttACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGCAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCg                                                                             >  2:454959/1‑140 (MQ=255)
                                                            aCCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAAc                                                                    >  2:94193/1‑140 (MQ=255)
                                                              cAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGt                                                                  >  1:577918/1‑140 (MQ=255)
                                                              cAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGt                                                                  <  2:577918/140‑1 (MQ=255)
                                                                   cgtccAGGAAAAACTCGACACGCTAATTGCGCGTCTGGATATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGgcg                                                             <  2:262254/140‑1 (MQ=255)
                                                                    gtccAGGAAAACCTCGCCCCGGTAATTGCGCGTCTGGCCATCCGCCGTGCACAGGTGCTGGCTGAGGAATTCATCGTTAAAGTTCAGGATGGAAATGGCCTAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGca                                                           <  2:99252/140‑1 (MQ=255)
                                                                     tccAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGca                                                           >  1:624073/1‑140 (MQ=255)
                                                                      ccAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGcac                                                          >  2:429472/1‑140 (MQ=255)
                                                                      ccAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGcac                                                          <  1:559832/140‑1 (MQ=255)
                                                                      ccAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGcac                                                          >  2:559832/1‑140 (MQ=255)
                                                                            aaaaCTCGCCACGGTCATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAAt                                                    <  1:274693/140‑1 (MQ=255)
                                                                             aaaCTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAAtt                                                   <  1:488615/140‑1 (MQ=255)
                                                                                 tCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTAcc                                               >  2:644229/1‑140 (MQ=255)
                                                                                    ccACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAAt                                            >  2:380234/1‑140 (MQ=255)
                                                                                    ccACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAAt                                            <  1:494620/140‑1 (MQ=255)
                                                                                              tGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTAc                                  <  2:679825/140‑1 (MQ=255)
                                                                                                      tGGACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGAttttt                          >  1:59286/1‑140 (MQ=255)
                                                                                                       ggACATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTa                         >  2:689619/1‑140 (MQ=255)
                                                                                                         aCATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAAc                       <  1:158220/140‑1 (MQ=255)
                                                                                                         aCATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAAc                       <  1:727367/140‑1 (MQ=255)
                                                                                                          cATTCGCCGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACg                      <  1:167055/140‑1 (MQ=255)
                                                                                                                ccGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTgcgc                >  1:530125/1‑140 (MQ=255)
                                                                                                                 cGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCa               <  2:526615/140‑1 (MQ=255)
                                                                                                                 cGTGCACAGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCa               <  1:689619/140‑1 (MQ=255)
                                                                                                                 cGTGCACAGGTGCTGGCTGAGGAAATCATCGCTGACGTTAAGGATGGACATGGACTAAACCTCGGCGTGCAATGGCCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCa               <  2:229052/140‑1 (MQ=255)
                                                                                                                     cacaGGTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGgt           >  2:543541/1‑140 (MQ=255)
                                                                                                                          gTGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGc      <  2:530125/140‑1 (MQ=255)
                                                                                                                           tGCTGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCt     >  1:82748/1‑140 (MQ=255)
                                                                                                                              tGGCTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGAt  <  1:429472/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                        
CGATGGATAACGTCGCCATTACCGCCGATGAACAGACTAACTCTCTGGTCATTACCGCTGACCAGTCCGTCCAGGAAAAACTCGCCACGGTAATTGCGCGTCTGGACATTCGCCGTGCACAGGTGCTGGTTGAGGCAATCATCGTTGAAGTTCAGGATGGAAATGGACTAAACCTCGGCGTGCAATGGGCGAATAAAAACGTTGGCGCACAGCAATTTACCAATACCGGATTACCGATTTTTAACGCTGCGCAAGGTGTGGCTGAT  >  NC_000913/3457296‑3457561

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: