Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,447,891 G→A *1260* (TAG→TAA tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,447,8910GA98.5% 214.5 / NA 67*1260* (TAG→TAAtamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/34);  minor base T (1/0);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGC  >  NC_000913/4447782‑4448024
                                                                                                               |                                                                                                                                     
aCTCTATAGCAACACTCACGTTACGTTATAGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACGAGGCACTGGATTGCTCTATCAGTGCGAGTTTTAACAATGCGAATATTTATCTACGGCAGTTTAc                                                                                                         >  2:535299/1‑139 (MQ=255)
    tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                                     >  2:308918/1‑139 (MQ=255)
    tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                                     >  1:27411/1‑139 (MQ=255)
    tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                                     >  1:412833/1‑139 (MQ=255)
    tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCATGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                                     >  1:214976/1‑139 (MQ=255)
      taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                                                   >  2:425397/1‑139 (MQ=255)
      taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                                                   >  2:700976/1‑139 (MQ=255)
        gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                                                 <  2:490612/139‑1 (MQ=255)
          aaCACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                               <  1:700976/139‑1 (MQ=255)
          aaCACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                               <  1:434344/139‑1 (MQ=255)
            cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                                             >  2:351087/1‑139 (MQ=255)
              cTCGCGTTACGTTATCGCCTGATGCCTAAGACATATCTGGCAGCCGTGTCTG‑‑GAGTAGGCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAAGAATGCGAATATTTCTCTACGGACCTTTACGCCACAAACAAgtc                                                                                           >  1:238079/1‑137 (MQ=255)
              cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                                           >  2:410411/1‑139 (MQ=255)
              cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                                           >  1:590370/1‑139 (MQ=255)
                cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         <  2:461084/139‑1 (MQ=255)
                 acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                        >  2:539962/1‑139 (MQ=255)
                    ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                                     <  2:27411/139‑1 (MQ=255)
                       cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                  >  1:635629/1‑139 (MQ=255)
                            tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                             <  1:354855/139‑1 (MQ=255)
                            tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                             <  1:593026/139‑1 (MQ=255)
                                 tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                                                          <  1:70688/139‑1 (MQ=255)
                                      ccTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                   <  1:522650/139‑1 (MQ=255)
                                                 tCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGACACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                        >  2:327489/1‑139 (MQ=255)
                                                   tGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                                                      <  2:241720/139‑1 (MQ=255)
                                                    ggAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                                                     >  1:679999/1‑139 (MQ=255)
                                                      aaGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAtt                                                   >  2:89056/1‑139 (MQ=255)
                                                         ccGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                                >  2:542513/1‑139 (MQ=255)
                                                         ccGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                                >  1:439789/1‑139 (MQ=255)
                                                         ccGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                                >  2:704936/1‑139 (MQ=255)
                                                         ccGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACCGCAGTCTACGCCACAAACATGGCAACAGTCACTGTATGACCAATGCCCAGTTACTCGGCGATTGCa                                                >  1:50075/1‑139 (MQ=255)
                                                          cGTGTCTG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAg                                               <  2:797880/139‑1 (MQ=255)
                                                                tG‑‑GTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa                                         <  1:410411/139‑1 (MQ=255)
                                                                 g‑‑gTGTAGACCATGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                        <  2:519351/139‑1 (MQ=255)
                                                                         gACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                  <  2:256637/139‑1 (MQ=255)
                                                                          aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                                                          >  1:715696/1‑98 (MQ=255)
                                                                          aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAc                                                                          <  2:715696/98‑1 (MQ=255)
                                                                           ccAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACGAGGCAACAGGCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                                <  1:79288/139‑1 (MQ=255)
                                                                           ccAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                                <  1:285787/139‑1 (MQ=255)
                                                                           ccAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                                <  2:50075/139‑1 (MQ=255)
                                                                           ccAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                                <  1:542513/139‑1 (MQ=255)
                                                                            cAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg                               <  2:429972/139‑1 (MQ=255)
                                                                              ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                   >  2:235430/1‑37 (MQ=255)
                                                                              ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                   <  1:235430/37‑1 (MQ=255)
                                                                              ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                                >  2:199/1‑136 (MQ=255)
                                                                              ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                                <  1:199/136‑1 (MQ=255)
                                                                               gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         <  2:89436/78‑1 (MQ=255)
                                                                               gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         >  1:89436/1‑78 (MQ=255)
                                                                               gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         <  1:412077/78‑1 (MQ=255)
                                                                               gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         >  2:412077/1‑78 (MQ=255)
                                                                                 aCTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATTACTACCAGTTGta                          <  1:539962/139‑1 (MQ=255)
                                                                                     ggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         <  1:225448/70‑1 (MQ=255)
                                                                                     ggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         >  2:225448/4‑73 (MQ=255)
                                                                                         tGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTg                   <  2:777661/139‑1 (MQ=255)
                                                                                         tGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTg                   <  1:670471/139‑1 (MQ=255)
                                                                                           ctctATCAGTTCGAGTTTTATCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTggg                 >  1:469785/1‑139 (MQ=255)
                                                                                              tATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCa              <  1:351087/139‑1 (MQ=255)
                                                                                                tCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                                                   <  2:538502/52‑1 (MQ=255)
                                                                                                tCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                                                   >  1:538502/1‑52 (MQ=255)
                                                                                                     tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                    <  1:531943/110‑1 (MQ=255)
                                                                                                     tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                    >  2:531943/1‑110 (MQ=255)
                                                                                                     tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCa       >  2:587374/1‑139 (MQ=255)
                                                                                                        aGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                          <  1:286460/117‑1 (MQ=255)
                                                                                                        aGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                          >  2:286460/1‑117 (MQ=255)
                                                                                                         gTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcg   <  1:327489/139‑1 (MQ=255)
                                                                                                          ttttAACGATGCGAATATTTGTCTACGGCAGTTTACGCCATAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc  >  2:488540/1‑139 (MQ=255)
                                                                                                          ttttAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         <  2:237088/52‑1 (MQ=255)
                                                                                                          ttttAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                         >  1:237088/1‑52 (MQ=255)
                                                                                                               |                                                                                                                                     
ACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTG‑‑GTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGC  >  NC_000913/4447782‑4448024

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