Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,457,801 G→A G475G (GGG→GGA)  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,8010GA97.7% 147.5 / NA 44G475G (GGG→GGA) gspDType II secretion system protein GspD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/20);  minor base T (0/1);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGGGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGTTGACTCTTCCTCTAACTCG  >  NC_000913/3457664‑3457936
                                                                                                                                         |                                                                                                                                       
cTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGgaga                                                                                                                                       >  2:437004/1‑140 (MQ=255)
  acTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATa                                                                                                                                     <  2:136499/139‑1 (MQ=255)
   cTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCAGGCCACAGAACACTTCAGGAGATaa                                                                                                                                    >  1:219248/1‑140 (MQ=255)
    ttACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAAc                                                                                                                                   <  1:35225/140‑1 (MQ=255)
    ttACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCACGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAAc                                                                                                                                   <  2:469834/140‑1 (MQ=255)
        cgcgCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCt                                                                                                                               >  1:150741/1‑140 (MQ=255)
                                  aCATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTg                                                                                                     >  1:2171/1‑140 (MQ=255)
                                  aCATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTg                                                                                                     >  2:346163/1‑140 (MQ=255)
                                       cTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGAc                                                                                                <  2:150741/140‑1 (MQ=255)
                                       cTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGAc                                                                                                <  1:519173/140‑1 (MQ=255)
                                         cGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACaa                                                                                              <  2:101125/140‑1 (MQ=255)
                                             aCCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACt                                                                                          <  1:104628/140‑1 (MQ=255)
                                             aCCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACt                                                                                          <  2:75540/140‑1 (MQ=255)
                                                   aGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCCAAACGGTGGGGACAAAACTCAAAGt                                                                                    <  2:452568/140‑1 (MQ=255)
                                                        cGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTc                                                                               >  2:15330/1‑140 (MQ=255)
                                                          tAACGCTGGATAATAAACTCGCTTCCTTCAACGTGTGGCGGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCg                                                                             <  2:219248/140‑1 (MQ=255)
                                                             cGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAg                                                                          <  1:359564/140‑1 (MQ=255)
                                                                  gATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCaa                                                                     >  1:374874/1‑140 (MQ=255)
                                                                  gATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCaa                                                                     >  2:29583/1‑140 (MQ=255)
                                                                     aataaACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGa                                                                  <  2:17232/140‑1 (MQ=255)
                                                                       taaACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAg                                                                >  2:111715/1‑140 (MQ=255)
                                                                        aaaCTCGCGTCCTTCAACGTGGTGCGGGAGGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAgg                                                               <  1:261218/140‑1 (MQ=255)
                                                                              gcgTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGAcgc                                                         >  1:162255/1‑140 (MQ=255)
                                                                                   cTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGgtgt                                                    >  1:379232/1‑140 (MQ=255)
                                                                                     tCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAGAGTTACTACGCAGGACAATGAAGGCGACGCGGTGTTg                                                  >  1:211434/1‑140 (MQ=255)
                                                                                       aaCGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCt                                                <  2:7885/140‑1 (MQ=255)
                                                                                            ggggCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGaaa                                           >  2:664077/1‑140 (MQ=255)
                                                                                               gCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATag                                        <  1:346163/140‑1 (MQ=255)
                                                                                                    aTGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGAGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAg                                   >  1:73653/1‑140 (MQ=255)
                                                                                                     tgtgCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTGCTCCGCAGCTCAATGACGGCGACGCGGTGTTGCTCGAAATAGAGCAgg                                  >  2:216565/1‑140 (MQ=255)
                                                                                                      gtgCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGa                                 <  1:459171/140‑1 (MQ=255)
                                                                                                        gCCGGTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAg                               >  1:648668/1‑140 (MQ=255)
                                                                                                           ggTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGtct                            <  2:2171/140‑1 (MQ=255)
                                                                                                            gTGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAATGCGACGCGGTGGTGCTCGAAACAGAGCAGGAAGtctc                           >  1:240721/1‑140 (MQ=255)
                                                                                                             tGCTATCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTcc                          <  1:437004/140‑1 (MQ=255)
                                                                                                                 aTCCGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCg                      >  1:443768/1‑140 (MQ=255)
                                                                                                                   ccGGGTCACAGACCACTTCAGGTGATAACGTGTTTGATATCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGtt                    <  2:232518/140‑1 (MQ=255)
                                                                                                                   ccGGGTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGtt                    >  2:131472/1‑140 (MQ=255)
                                                                                                                       gTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGTAATGGAACAGGAATTCTCCAGCGTTGAct                >  1:232518/1‑140 (MQ=255)
                                                                                                                       gTCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGTTGAct                >  2:621187/1‑140 (MQ=255)
                                                                                                                        tCACAGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGTTGActc               <  1:15330/140‑1 (MQ=255)
                                                                                                                         cacaGACCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGTTGActct              <  1:158967/140‑1 (MQ=255)
                                                                                                                              aCCACTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGTTGACTCTTCctc         >  2:660508/1‑140 (MQ=255)
                                                                                                                                  cTTCAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGTTGACTCTTCCTCTAAc     >  2:349735/1‑140 (MQ=255)
                                                                                                                                     cAGGAGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGCGCAGGAAATCTGCAGCGTTGAATCTTCCTCTAACTCg  >  1:202272/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCCTTCAACGTGGGGCAGGATGTGCCGGTGCTATCCGGGTCACAGACCACTTCAGGGGATAACGTCTTTAATACCGTCGAACGCAAAACGGTGGGGACAAAACTCAAAGTTACTCCGCAGGTCAATGAAGGCGACGCGGTGTTGCTCGAAATAGAGCAGGAAGTCTCCAGCGTTGACTCTTCCTCTAACTCG  >  NC_000913/3457664‑3457936

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: ‑