Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,010,497 T→C 7.4% I41T (ATT→ACT)  prsD_1 → Type I secretion system ATP‑binding proteinPrsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,010,4970TC 7.4% 166.0 / 4.5 54I41T (ATT→ACT) prsD_1Type I secretion system ATP‑binding proteinPrsD
Reads supporting (aligned to +/- strand):  ref base T (20/30);  new base C (2/2);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

CGGGAGCGGGCGTGTTCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATGGGTATGCTCGAAACACTC  >  ACB122/2010365‑2010624
                                                                                                                                    |                                                                                                                               
cgtgAGCGGGCGTGTTCGCATTTTTCGACGCTCCGTTGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTTGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGAtt                                                                                                                                <  3:69078/131‑1 (MQ=255)
cGGGAGCGGGCGTGTTCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGAtt                                                                                                                                <  8:248116/134‑1 (MQ=255)
cGGGAGCGGGCGTGTTCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGAtt                                                                                                                                <  3:74283/134‑1 (MQ=255)
cGGGAGCGGGCGTGTTCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTa                                                                                                                               <  7:38762/135‑1 (MQ=255)
cGGGAGCGGGCGTGTTAGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGAtt                                                                                                                                <  2:134018/134‑1 (MQ=255)
         gCGTGTTCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATcgcg                                                                                                                      <  1:294151/135‑1 (MQ=255)
               tCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACgg                                                                                                                <  7:257771/135‑1 (MQ=255)
               tCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACgg                                                                                                                <  4:369860/135‑1 (MQ=255)
                    tttttCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAg                                                                                                           <  5:49450/135‑1 (MQ=255)
                    tttttCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAg                                                                                                           <  8:267229/135‑1 (MQ=255)
                                              ttATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGaaaa                                                                                 <  8:62437/135‑1 (MQ=255)
                                              ttATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGaaaa                                                                                 <  3:281895/135‑1 (MQ=255)
                                                ataCAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAAc                                                                                <  1:349980/134‑1 (MQ=255)
                                                         tcATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTa                                                                      >  3:185766/1‑135 (MQ=255)
                                                           aTGTTGCTATTCCATCCATGGTTTGGTTGGATTACGTTGGCCTGTGGATCAATTTTTACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAAt                                                                    >  5:127628/1‑135 (MQ=255)
                                                             gTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAAtgt                                                                  <  8:247494/135‑1 (MQ=255)
                                                             gTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAAtgt                                                                  <  1:363792/135‑1 (MQ=255)
                                                             gTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAAtgt                                                                  <  7:120053/135‑1 (MQ=255)
                                                             gTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAAtgt                                                                  <  3:263096/135‑1 (MQ=255)
                                                                   aTTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAAc                                                                                <  2:321108/115‑1 (MQ=255)
                                                                   aTTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAAc                                                                                >  1:321108/1‑115 (MQ=255)
                                                                        atccatGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCt                                                        <  2:347285/134‑1 (MQ=255)
                                                                           catGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGa                                                    >  7:309977/1‑135 (MQ=255)
                                                                           catGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGa                                                    >  2:360527/1‑135 (MQ=255)
                                                                            atGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAg                                                                            >  2:377631/1‑110 (MQ=255)
                                                                            atGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGAc                                                   >  2:353259/1‑135 (MQ=255)
                                                                            atGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGAc                                                   >  8:140256/1‑135 (MQ=255)
                                                                            atGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGAc                                                   >  2:319924/1‑135 (MQ=255)
                                                                            atGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGAc                                                   >  8:288326/1‑135 (MQ=255)
                                                                            atGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGCGTCAGGCTCTAGACATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGAc                                                   >  7:329729/1‑135 (MQ=255)
                                                                            atGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACAGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAg                                                                            <  1:377631/110‑1 (MQ=255)
                                                                                tttggttGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGtt                                                <  4:168066/134‑1 (MQ=255)
                                                                                tttggttGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGtt                                                <  5:83372/134‑1 (MQ=255)
                                                                                      tggATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAAc                                         >  8:265994/1‑135 (MQ=255)
                                                                                      tggATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGACAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAAc                                         >  4:147476/1‑135 (MQ=255)
                                                                                       ggATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACg                                        <  8:229193/135‑1 (MQ=255)
                                                                                       ggATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACg                                        >  7:12514/1‑135 (MQ=255)
                                                                                        gATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACg                                        <  8:12514/134‑1 (MQ=255)
                                                                                          tGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCg                                     <  8:76128/135‑1 (MQ=255)
                                                                                           gACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAAtgt                                                                  >  3:76936/1‑105 (MQ=255)
                                                                                           gACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAAtgt                                                                  <  4:76936/105‑1 (MQ=255)
                                                                                               ttGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTg                                <  1:319924/135‑1 (MQ=255)
                                                                                               ttGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTg                                <  7:307800/135‑1 (MQ=255)
                                                                                               ttGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTg                                <  7:140256/135‑1 (MQ=255)
                                                                                                  gCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCGATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGAt                              <  4:145541/134‑1 (MQ=255)
                                                                                                          aTCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCAt                     >  1:26481/1‑135 (MQ=255)
                                                                                                          aTCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCAt                     >  6:127684/1‑135 (MQ=255)
                                                                                                          aTCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCAt                     >  6:138487/1‑135 (MQ=255)
                                                                                                          aTCAATTTTGACCGTCTTTGCAGTGACTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTg                                >  7:34000/1‑124 (MQ=255)
                                                                                                          aTCAATTTTGACCGTCTTTGCAGTGACTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTg                                <  8:28074/124‑1 (MQ=255)
                                                                                                          aTCAATTTTGACCGTCTTTGCAGTGACTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTg                                <  8:34000/124‑1 (MQ=255)
                                                                                                          aTCAATTTTGACCGTCTTTGCAGTGACTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTg                                >  7:28074/1‑124 (MQ=255)
                                                                                                           tCAATTTTAACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATg                    <  7:24112/135‑1 (MQ=255)
                                                                                                                     ccGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATGGGTATGCTCg          <  1:353259/135‑1 (MQ=255)
                                                                                                                       gTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATGGGTATGCTCGaa        <  4:185766/135‑1 (MQ=255)
                                                                                                                           ttGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATGGGTATGCTCGAAacac    >  1:334765/1‑135 (MQ=255)
                                                                                                                            tGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATGGGTATGCTCGAAacac    >  5:268462/1‑134 (MQ=255)
                                                                                                                              cAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATGGGTATGCTCGAAACACTc  <  7:267081/134‑1 (MQ=255)
                                                                                                                                    |                                                                                                                               
CGGGAGCGGGCGTGTTCGCATTTTTCGACGCTCCGTGGTTTCCCATTTATACAGTCGTCATGTTCCTATTCCATCCATGGTTTGGTTGGATGACGTTGGCCTGTGGATCAATTTTGACCGTCTTTGCAGTGATTAATCATCGCGTGACGGGTCAGGCTCTAGCCATAGCGAACAAAGAAAACGTAGCGTCTAATGTGATCACGTCTAAGACGTTGCGTAACGCCGAAGTGATCGAGTCCATGGGTATGCTCGAAACACTC  >  ACB122/2010365‑2010624

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑