Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 560,665 G→A 7.0% L323L (TTG→TTA)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913560,6650GA 7.0% 183.9 / 3.4 57L323L (TTG→TTA) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

ATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGG  >  NC_000913/560520‑560801
                                                                                                                                                 |                                                                                                                                        
ataataTGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGc                                                                                                                                         <  7:22238/147‑1 (MQ=255)
          tCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCt                                                                                                                              <  3:161738/148‑1 (MQ=255)
           ccTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCAGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCtt                                                                                                                             >  5:312530/1‑148 (MQ=255)
            cTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTc                                                                                                                            <  6:76932/148‑1 (MQ=255)
               aTAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCaa                                                                                                                          <  7:251798/147‑1 (MQ=255)
                  gCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATgg                                                                                                                      <  5:174063/148‑1 (MQ=255)
                  gCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATg                                                                                                                       <  5:145715/147‑1 (MQ=255)
                    cAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGc                                                                                                                     >  6:251933/1‑147 (MQ=255)
                    cAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGc                                                                                                                     >  8:221015/1‑147 (MQ=255)
                     aGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGAATTGCACCCGACATCTTCAAATGGCGa                                                                                                                   <  3:219977/148‑1 (MQ=255)
                     aGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGAATTGCACCCGACATCTTCAAATGGCGa                                                                                                                   <  3:223886/148‑1 (MQ=255)
                      gCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGAt                                                                                                                  <  7:221015/148‑1 (MQ=255)
                       cAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATc                                                                                                                 >  5:211035/1‑148 (MQ=255)
                           ggTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGa                                                                                                             >  6:42814/1‑148 (MQ=255)
                                   ccAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCa                                                                                                     >  4:19921/1‑148 (MQ=255)
                                          tACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATcg                                                                                                  <  7:126954/144‑1 (MQ=255)
                                          tACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATcg                                                                                                  >  8:126954/1‑144 (MQ=255)
                                                 gATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGAt                                                                                       <  6:307266/148‑1 (MQ=255)
                                                      cccGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCa                                                                                   <  6:312530/147‑1 (MQ=255)
                                                           aCGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAAt                                                                                 <  7:161543/144‑1 (MQ=255)
                                                           aCGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAAt                                                                                 >  8:161543/1‑144 (MQ=255)
                                                                   ccAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTac                                                                     <  5:42814/148‑1 (MQ=255)
                                                                   ccAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTa                                                                      >  1:273825/1‑147 (MQ=255)
                                                                          aCGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTc                                                              >  1:305946/1‑148 (MQ=255)
                                                                           cGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGcgcg                                                                                         <  6:143552/120‑1 (MQ=255)
                                                                           cGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTACCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGcgcg                                                                                         >  5:143552/1‑120 (MQ=255)
                                                                               tCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTAt                                                         >  2:94069/1‑148 (MQ=255)
                                                                                  gCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCa                                                      <  6:211035/148‑1 (MQ=255)
                                                                                    aaaaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAAc                                                    <  2:273825/148‑1 (MQ=255)
                                                                                    aaaaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAAc                                                    <  1:90446/148‑1 (MQ=255)
                                                                                     aaaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCTATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCaa                                                     >  1:100495/1‑146 (MQ=255)
                                                                                     aaaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCaa                                                     >  8:40367/1‑146 (MQ=255)
                                                                                     aaaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACa                                                   <  3:19921/148‑1 (MQ=255)
                                                                                      aaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGAt                                                                                       <  5:314659/111‑1 (MQ=255)
                                                                                      aaTGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGAt                                                                                       >  6:314659/1‑111 (MQ=255)
                                                                                             ttattaTCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTACACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAAt                                           >  1:168322/1‑148 (MQ=255)
                                                                                             ttattaTCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTACACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAAt                                           >  1:173202/1‑148 (MQ=255)
                                                                                                        cAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTAcaca                                                                  <  4:295220/114‑1 (MQ=255)
                                                                                                        cAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTAcaca                                                                  >  3:295220/1‑114 (MQ=255)
                                                                                                         aaaGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGa                               <  1:94069/148‑1 (MQ=255)
                                                                                                         aaaGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGa                               <  7:40367/148‑1 (MQ=255)
                                                                                                           aGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTa                                       <  8:248563/138‑1 (MQ=255)
                                                                                                           aGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTa                                       >  7:248563/1‑138 (MQ=255)
                                                                                                                tGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGt                        >  5:71672/1‑148 (MQ=255)
                                                                                                                 gTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGtt                       >  4:185930/1‑148 (MQ=255)
                                                                                                                    tCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTc                     >  4:15297/1‑147 (MQ=255)
                                                                                                                      ccccGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCaaa                  <  5:164610/148‑1 (MQ=255)
                                                                                                                          ggCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAAttt               >  8:299004/1‑147 (MQ=255)
                                                                                                                            cgcTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTAc                                      <  5:280224/122‑1 (MQ=255)
                                                                                                                            cgcTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTAc                                      >  6:280224/1‑122 (MQ=255)
                                                                                                                            cgcTTTTGAAATTACAGATTTACACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGAc            <  2:168322/148‑1 (MQ=255)
                                                                                                                                 ttGAAATTACAGATTTACACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGAc            <  2:173202/143‑1 (MQ=255)
                                                                                                                                    aaaTTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACaa                                    >  7:82892/1‑116 (MQ=255)
                                                                                                                                    aaaTTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACaa                                    <  8:82892/116‑1 (MQ=255)
                                                                                                                                    aaaTTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGAcgg     <  5:251933/147‑1 (MQ=255)
                                                                                                                                      aTTACAGATTTGCACCCGACATCTTCAAATGGCTATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGAcggcgg  <  2:100495/148‑1 (MQ=255)
                                                                                                                                      aTTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGAcggcgg  >  5:277918/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
ATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTTTTGAAATTACAGATTTGCACCCGACATCTTCAAATGGCGATCTGGACGTCACCATCGACGAGCGCGATGGCAATCAGCAGAATTACACAATTCCGTATTCAACAGTGCCAATTTTACAACGCGAAGGGCGTTTCAAATTTGACCTGACGGCGG  >  NC_000913/560520‑560801

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