Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,544,276 T→C 8.5% Y168Y (TAT→TAC hlgC → gamma hemolysin component C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,544,2760TC 8.5% 143.3 / 4.1 59Y168Y (TAT→TAChlgCgamma hemolysin component C
Reads supporting (aligned to +/- strand):  ref base T (33/18);  new base C (2/3);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AATTAATTATTTACCTAAAAATAAAATTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCT  >  CP000730/2544138‑2544414
                                                                                                                                          |                                                                                                                                           
aattaattATTTACCTAAAAATAAAATTGTATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCtata                                                                                                                                            <  1:233801/140‑1 (MQ=255)
 attaattaTTTACCTAAAAATAAAATTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATac                                                                                                                                           <  1:688202/140‑1 (MQ=255)
             ccTAAAAATAAAATTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTACCTATACACAACAAAATTa                                                                                                                               >  2:138296/1‑140 (MQ=255)
               tAAAAATAAAATTGCACCTTCAAACGTGAGTCAGAGATTAGGATACAAGATCGGTGGTAATTTCCAATTAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACTAAAtgttg                                                                                                                             <  2:89681/140‑5 (MQ=255)
               tAAAAATAAAATTGAATCTACAAACGTGAGTTAGATATTAGGATACAATATAGGTGGTAATTTCCAATCAGCCCCATCACTAGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAATTTATg                                                                                                                             <  2:738691/140‑1 (MQ=255)
                aaaaaTAAAATTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGt                                                                                                                            >  2:103423/1‑140 (MQ=255)
                 aaaataaaatTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCCATACACAACAAAATTaagga                                                                                                                           >  2:264425/1‑136 (MQ=255)
                      aaaatTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGa                                                                                                                      >  1:500535/1‑140 (MQ=255)
                      aaaatTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGGGAGAATTTCCAATCAGCCCCAACACTCGGAGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACCACAAAATTATGTAAGTGa                                                                                                                      >  1:400649/1‑140 (MQ=255)
                       aaatTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGaa                                                                                                                     <  1:757541/140‑1 (MQ=255)
                                 tACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTACACACAACAAAATTATGTAAGTGAAGTAGaacaac                                                                                                           >  2:367675/1‑140 (MQ=255)
                                 tACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTACACACAACAAAATTATGTAAGTGAAGTAGaacaac                                                                                                           >  2:325095/1‑140 (MQ=255)
                                   cAAACGTGAGTCAGACATTAGGATACACTATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAAcaacaa                                                                                                         <  1:172988/140‑1 (MQ=255)
                                    aaaCGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATACTAAATCGATTAGCTATACACAAAAAAATTATGTAAGTGAAGTAGAACaaaaaa                                                                                                        >  2:395057/1‑140 (MQ=255)
                                      aCGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAAc                                                                                                      <  1:80104/140‑1 (MQ=255)
                                         tGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACACGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCa                                                                                                   >  1:507019/1‑140 (MQ=255)
                                            gTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTAAAGTCGAACAACAAAACTCaaaa                                                                                                >  2:581570/1‑140 (MQ=255)
                                                   aTTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTACACACAACAAAATTATGTACGAGAAGAAGAACAACAAAACTCAAAAAGTGttt                                                                                         <  1:367675/140‑1 (MQ=255)
                                                   aTTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTACACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGttt                                                                                         <  1:325095/140‑1 (MQ=255)
                                                    ttataaTACAAAATCGGTGGTAATCTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCAATTAGCCATACACAACAAAATTAGGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGtttt                                                                                        <  1:64590/135‑1 (MQ=255)
                                                    ttAGGATACGATATCGGTAGTAATTTCTAATCAGCCCCATCACTCTGTTGTAATGGATCATTTAACTATTCTAAAACCATCAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGtttt                                                                                        <  1:441087/140‑1 (MQ=255)
                                                     tAGGATACAAAATCGGTGGTAATTTCCAATCAACCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGAATAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTa                                                                                       <  1:124554/140‑1 (MQ=255)
                                                           aCAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACGATGCTAAATCGATTAGCTATACACAACAAAAAGATGTAATTTAAGAATAAAAACAAAACTCagaaaatcatgtattcgac                                                                                 >  1:343443/1‑126 (MQ=255)
                                                                atCGGTGGTAGTTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTaa                                                                           >  2:514077/1‑140 (MQ=255)
                                                                atCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTGTCTGGGGCGTTaa                                                                           >  2:744559/1‑140 (MQ=255)
                                                                atCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTATAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTaa                                                                           >  2:66587/1‑140 (MQ=255)
                                                                 tCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTAAGTAAGTGAAGTAGAACAACAACACTCAAAAAGTGTTTTATGGGACGTTaaa                                                                          >  2:31044/1‑140 (MQ=255)
                                                                          aTTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattca                                                                 >  2:451060/1‑140 (MQ=255)
                                                                          aTTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAAAAACAAAACTCAAAAAGGGTTTTATGGGGCGTTAAAGCGAattaa                                                                 >  2:20823/1‑138 (MQ=255)
                                                                            ttCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACCAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattcatt                                                               >  2:296749/1‑140 (MQ=255)
                                                                            ttCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATGAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGCTTTATGGGGCGTTACAGCGAattcatt                                                               >  2:442444/1‑140 (MQ=255)
                                                                             tCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattcattc                                                              >  2:715276/1‑140 (MQ=255)
                                                                             tCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAattcattc                                                              >  2:729093/1‑140 (MQ=255)
                                                                                  tCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAc                                                         >  2:752121/1‑140 (MQ=255)
                                                                                  tCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAc                                                         >  1:194888/1‑140 (MQ=255)
                                                                                  tCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAAAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAc                                                         >  1:194937/1‑140 (MQ=255)
                                                                                    aGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAATTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTg                                                       >  2:169903/1‑140 (MQ=255)
                                                                                    aGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTg                                                       <  2:663571/140‑1 (MQ=255)
                                                                                       cccATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAAt                                                    >  2:362631/1‑140 (MQ=255)
                                                                                         cATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCa                                                  >  1:759357/1‑140 (MQ=255)
                                                                                            tacccGGTAGTAATGGAACAATTAACTATTCAAAATCAATCACCAAAACAAAACCTAAATATATAAAACAAATAAAAAAATAAAAATAAAAAAGTATTTTATGGCGCGTTAAAGCTAATTCATTCGCCACTGAATCAGGt                                               <  2:283898/136‑1 (MQ=255)
                                                                                                 ggtggtAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAATGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCaaaa                                          >  1:102073/1‑140 (MQ=255)
                                                                                                    ggtAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAATTCCTTTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATc                                       <  1:590032/140‑1 (MQ=255)
                                                                                                        aTGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGcc                                   >  1:96729/1‑140 (MQ=255)
                                                                                                        aTGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGcc                                   <  1:567271/140‑1 (MQ=255)
                                                                                                          ggATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGAGAATTAATACGCCCCTGAATCAGAACAAAAagcagcata                                 >  2:716535/1‑137 (MQ=255)
                                                                                                                  tAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCg                         >  1:454680/1‑140 (MQ=255)
                                                                                                                   aaCTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGa                        <  1:598544/140‑1 (MQ=255)
                                                                                                                   aaCTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGAATAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGa                        >  1:589925/1‑140 (MQ=255)
                                                                                                                       aatCTAAATCGATTAACAAAACACAACAAAATTATGTAAGTGAAGTAGAAAAAAAAAAAAAAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGattta                    <  2:446256/138‑1 (MQ=255)
                                                                                                                       aTTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAATAAGTGTTTTATTGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGattta                    <  2:428890/140‑1 (MQ=255)
                                                                                                                        ttCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATACGCCACTGAATCAGGTCAAAAATCAGCCTTTGATGGCGatttat                   >  2:171786/1‑140 (MQ=255)
                                                                                                                         tCTAAATCGATTATCTTTCCACAACAAAATTATGTAAGTGAAGTATAATTACTTCACTCTTAAAGTGTTTTATGGGGCGTTAAAGTGAATTTATTCGTCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGatttatt                  <  1:670580/140‑1 (MQ=255)
                                                                                                                                cGATTAGCTACACAGAACCAAGTTATACAAGAGAAGNAGAACAACAACCCTCAAGAAGCATCTAACGGCGGTGTAAAGCTAATTCATTGGACAAAAAATCAAAACAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGc           <  2:1933/140‑1 (MQ=255)
                                                                                                                                 gATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCt          <  1:677040/140‑1 (MQ=255)
                                                                                                                                  aTTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTa         <  1:339599/140‑1 (MQ=255)
                                                                                                                                  aTTAGCTAGACACAACTAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTa         <  2:392039/140‑1 (MQ=255)
                                                                                                                                    tAGCTATACAAAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACa       <  1:427427/140‑1 (MQ=255)
                                                                                                                                      gCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGCATTACTTCGCCACTGCATCCGGACACAACTCAGACTTGGATAGGGATTTAATTGTAGGCTACaaa     >  2:522635/1‑140 (MQ=255)
                                                                                                                                      gCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGCCAATGCATCAGGTCAGCAATCGGCCTTTGAGAGCGAGTTATCTGTAGGCTACaag     >  2:522654/1‑139 (MQ=255)
                                                                                                                                        tataCACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAAcc   >  2:289957/1‑140 (MQ=255)
                                                                                                                                         ataCACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTTAAAACGAATTCATTCGCCACTAAATAAGGTCAAAAAGCAGCCATTGATAGCGGTTTATTTGTAGCCTACAAACCt  >  1:235177/1‑140 (MQ=255)
                                                                                                                                         ataCACAACAAAATTATGTAAGTGAACTAGAACAACATAAAGCAAACAGTGTTTTATGGGGCGTTAAAGCGAATTCATTCGACACTGAATCAGGTCACAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCt  <  2:728629/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
AATTAATTATTTACCTAAAAATAAAATTGAATCTACAAACGTGAGTCAGACATTAGGATACAATATCGGTGGTAATTTCCAATCAGCCCCATCACTCGGTGGTAATGGATCATTTAACTATTCTAAATCGATTAGCTATACACAACAAAATTATGTAAGTGAAGTAGAACAACAAAACTCAAAAAGTGTTTTATGGGGCGTCAAAGCGAATTCATTCGCCACTGAATCAGGTCAAAAATCAGCCTTTGATAGCGATTTATTTGTAGGCTACAAACCT  >  CP000730/2544138‑2544414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: