Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,389,156 G→A V77I (GTA→ATA)  gspG → pseudopilin, cryptic, general secretion pathway

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,389,1560GA98.4% 711.8 / ‑0.4 187V77I (GTA→ATA) gspGpseudopilin, cryptic, general secretion pathway
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (95/89);  total (98/89)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 239 total aligned reads displayed.

GGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTCTGCCTGCCG  >  REL606/3389083‑3389228
                                                                         |                                                                        
 gAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTtaatag                                                                       <  1:5985586/76‑1 (MQ=255)
  aaaaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGa                                                                      <  2:1452152/76‑1 (MQ=255)
  aaaaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGa                                                                      >  2:1462940/1‑76 (MQ=255)
  aaaaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGa                                                                      <  1:2461915/76‑1 (MQ=255)
   aaaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATaaaa                                                                     >  2:5694413/1‑76 (MQ=255)
   aaaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGaa                                                                     <  1:1303602/76‑1 (MQ=255)
   aaaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGaa                                                                     >  1:5244153/1‑76 (MQ=255)
    aaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAg                                                                    <  1:1138463/76‑1 (MQ=255)
    aaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAg                                                                    >  2:4029787/1‑76 (MQ=255)
    aaCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAg                                                                    >  1:6800872/1‑76 (MQ=255)
     aCGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAgc                                                                   <  1:1586633/76‑1 (MQ=255)
      cGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAgcg                                                                  <  2:4538014/76‑1 (MQ=255)
      cGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAgcg                                                                  <  1:2219695/76‑1 (MQ=255)
      cGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAgcg                                                                  <  2:6370609/76‑1 (MQ=255)
        cccTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGcc                                                                >  2:1568954/1‑76 (MQ=255)
         ccTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCg                                                               >  1:4362207/1‑76 (MQ=255)
          cTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGa                                                              <  2:322515/76‑1 (MQ=255)
              atatGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACt                                                          >  1:2377055/1‑76 (MQ=255)
              atatGTACAAACTCGACAATCATCACTACCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACt                                                          <  2:869327/76‑1 (MQ=255)
              atatGTAAAACCTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACt                                                          <  2:2403197/76‑1 (MQ=255)
                atGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGc                                                        >  1:5456447/1‑76 (MQ=255)
                atGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGc                                                        <  1:1036118/76‑1 (MQ=255)
                atGTACAAACTCGACAAGCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGc                                                        <  2:102764/76‑1 (MQ=255)
                 tGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGcc                                                       >  1:7309878/1‑76 (MQ=255)
                 tGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGcc                                                       <  1:1470997/76‑1 (MQ=255)
                  gTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGGCGACACTGCCa                                                      >  2:6524725/1‑76 (MQ=255)
                   tACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCAc                                                     >  1:3367688/1‑76 (MQ=255)
                   tACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCAc                                                     <  1:7105176/76‑1 (MQ=255)
                    aCAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCAcc                                                    >  1:481513/1‑76 (MQ=255)
                    aCAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCAcc                                                    >  1:3196616/1‑76 (MQ=255)
                    aCAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCAcc                                                    <  1:7619448/76‑1 (MQ=255)
                      aaaCTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGc                                                  >  1:4925808/1‑76 (MQ=255)
                        aCTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTg                                                <  1:3633438/76‑1 (MQ=255)
                           cGACAACCATCACTACCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGcc                                             <  2:3022021/76‑1 (MQ=255)
                            gaaaaTCAGCACTATCCAACCACAAAGCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCg                                            <  2:3735151/73‑1 (MQ=255)
                             aaaaTCATCCCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGc                                           <  2:3506926/74‑1 (MQ=255)
                              cAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCa                                          >  2:2553664/1‑76 (MQ=255)
                               aaTCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCaa                                         >  1:6901558/1‑76 (MQ=255)
                               aaTCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCaa                                         >  1:2797513/1‑76 (MQ=255)
                               aaTCACCACTACCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCaa                                         <  1:6473209/76‑1 (MQ=255)
                                atcatcACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCaaa                                        >  2:3173946/1‑76 (MQ=255)
                                atcatcACTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCAGCaaa                                        >  2:4623774/1‑76 (MQ=255)
                                atcaacACAACCCAACCCCAAAACGGGGCCTTAAATCCTTAATAGAAGCGCCAACACTGCCACCGCTGGCCGCaaa                                        <  2:6268860/76‑1 (MQ=255)
                                  catccCTACCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACt                                      <  1:5270529/76‑1 (MQ=255)
                                  catcaCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACt                                      <  1:6653266/76‑1 (MQ=255)
                                  catcaCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACt                                      >  1:1423551/1‑76 (MQ=255)
                                  catcaCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACt                                      <  1:2984328/76‑1 (MQ=255)
                                  catcaCTACCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACt                                      <  1:5512182/76‑1 (MQ=255)
                                   atcaCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACta                                     >  2:2791870/1‑76 (MQ=255)
                                   atcaCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACta                                     >  2:654822/1‑76 (MQ=255)
                                   atcaCGACCCAACCACAAATCGGGGGCTTGACTCCTAAATAGAAGCGCCGACCCTGCCACCGCTGGCCGCAAACta                                     <  2:1168101/76‑1 (MQ=255)
                                    tcaCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACtat                                    >  2:7366203/1‑76 (MQ=255)
                                      ccTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATaa                                  <  1:1098654/75‑1 (MQ=255)
                                      ccTACCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATaa                                  <  1:3995016/75‑1 (MQ=255)
                                      aCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATaa                                  >  1:4120488/1‑76 (MQ=255)
                                      aCTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATaa                                  >  1:2811107/1‑76 (MQ=255)
                                       cTATCCAACCACAAATCAGGGGCTTGAATCCTTAATAGGAGGGCCGACACTGCCACCGCTGGCCGCAAAATATaaa                                 >  2:4090551/1‑75 (MQ=255)
                                       cTACCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAAc                                 <  1:4597543/76‑1 (MQ=255)
                                         acccAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaa                               <  2:2949304/74‑1 (MQ=255)
                                         aTCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaa                               >  1:6836948/1‑76 (MQ=255)
                                          tCCAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaag                              >  2:4949255/1‑76 (MQ=255)
                                          cccAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaag                              <  1:6325989/75‑1 (MQ=255)
                                           ccAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaagg                             <  2:5673917/76‑1 (MQ=255)
                                           ccAACCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaagg                             >  1:2542351/1‑76 (MQ=255)
                                             acccACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaaggaa                           <  2:2033001/74‑1 (MQ=255)
                                             aaCCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaaggaa                           >  2:717819/1‑76 (MQ=255)
                                             aaCCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaaggaa                           <  1:6365195/76‑1 (MQ=255)
                                             aaCCACAAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaaggaa                           <  1:6622658/76‑1 (MQ=255)
                                             aaCCACAAAGCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaaggaa                           <  2:4728595/76‑1 (MQ=255)
                                              aCCACAAATCAGGGGCTTGAATCCTTAAAAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACaaggaag                          <  2:831758/76‑1 (MQ=255)
                                                 acaAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGtt                       <  2:7061265/76‑1 (MQ=255)
                                                  caAATCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTta                      <  1:5550549/76‑1 (MQ=255)
                                                   aaaTCGGGGACTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTtat                     <  2:6156630/76‑1 (MQ=255)
                                                   aaaTCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTtat                     >  2:2740465/1‑76 (MQ=255)
                                                   aaaTCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTtat                     <  1:7572268/76‑1 (MQ=255)
                                                    aaTAAGGGCCTTAAATCTTAAATAAAAGCCCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTtata                    <  2:4309860/76‑1 (MQ=255)
                                                     aTCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTatat                   >  2:1891267/1‑76 (MQ=255)
                                                      taaGGGGGTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATc                  <  2:5930427/74‑1 (MQ=255)
                                                      tCAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATc                  <  1:6548431/76‑1 (MQ=255)
                                                      ccAGGGGCTGGAGTCCTGAATAGAAGCGCCGACACTGCCCCCGCTGGCCGCAAACTATAACAAGGAAGGTTATATc                  <  1:6418134/75‑1 (MQ=255)
                                                       cgggggCTGGAACCCTTAAAAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCa                 <  1:5011614/74‑1 (MQ=255)
                                                       cAGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCa                 >  1:1766709/1‑76 (MQ=255)
                                                        aGGGGCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCaa                >  2:6007823/1‑76 (MQ=255)
                                                        aGGGGCTTGAATCCTTAATAGAAGAGCCGACACTGCCACCGCTGGTCGCAAACTTTAACAAGGAAGGCTATATCaa                >  2:4373833/1‑76 (MQ=255)
                                                          gggCTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGc              <  1:3553345/76‑1 (MQ=255)
                                                             cTTGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTc           <  2:5540661/76‑1 (MQ=255)
                                                              ttaaaTCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTct          <  1:6121971/73‑1 (MQ=255)
                                                              ttGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTct          <  1:2746375/76‑1 (MQ=255)
                                                               tGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctg         >  2:6880093/1‑76 (MQ=255)
                                                               tGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctg         >  2:2483132/1‑76 (MQ=255)
                                                               tGAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctg         <  1:7356670/76‑1 (MQ=255)
                                                                gAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgc        >  2:4848501/1‑76 (MQ=255)
                                                                gAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgc        >  1:4451748/1‑76 (MQ=255)
                                                                gAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgc        >  1:376822/1‑76 (MQ=255)
                                                                gAATCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgc        >  1:5746495/1‑76 (MQ=255)
                                                                  aTCCTTAATAGACGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgcct      <  2:6585702/76‑1 (MQ=255)
                                                                  aTCCTTAATAAAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgcct      <  2:5282610/76‑1 (MQ=255)
                                                                   tctttAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgcctg     <  2:1276599/73‑1 (MQ=255)
                                                                   tCCTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTctgcctg     >  2:7520787/1‑76 (MQ=255)
                                                                     cTTAATAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTCtgcctgcc   >  2:5472134/1‑76 (MQ=255)
                                                                         |                                                                        
GGAAAACGCCCTTGATATGTACAAACTCGACAATCATCACTATCCAACCACAAATCAGGGGCTTGAATCCTTAGTAGAAGCGCCGACACTGCCACCGCTGGCCGCAAACTATAACAAGGAAGGTTATATCAAGCGTCTGCCTGCCG  >  REL606/3389083‑3389228

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