Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,986 A→T 7.3% intergenic (+49/+27) matP → / ← ompA macrodomain Ter protein/outer membrane porin A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9860AT 7.3% 145.6 / 2.1 53intergenic (+49/+27)matP/ompAmacrodomain Ter protein/outer membrane porin A
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base T (2/2);  total (22/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 = 7.71e-01

CAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCA  >  NC_000913/1018845‑1019121
                                                                                                                                             |                                                                                                                                       
cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGcagcag                                                                                                                                        >  2:586854/1‑143 (MQ=255)
cAGCGGATTGAAGATGCAGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGCATAACTGACGAAAGTCAGTTCAATTTACTAAAAGAAAAAAAAACCCCGcagcag                                                                                                                                        <  2:497161/143‑1 (MQ=255)
   cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTGCAATTTACTAAAGGC‑AAAAAAACCCCGCAGCAGCggg                                                                                                                                    <  1:597666/143‑1 (MQ=255)
    tGAGTGAAGATGAGGAAAAGCGAGAGAACGACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAAATTACTAAAGACAAAAAAAACCCCGCTGCTGCggg                                                                                                                                    <  2:747984/143‑1 (MQ=255)
       ttGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgggggt                                                                                                                                 >  1:359770/1‑141 (MQ=255)
          aaGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                              <  1:117380/143‑1 (MQ=255)
                     aaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTTAGGATTTACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                                                   <  1:738763/143‑1 (MQ=255)
                        aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                >  1:321313/1‑143 (MQ=255)
                             gaAATACGCGAATAAAATGTCTTCTCTGAAGCAGTATCTGCATGCATTTCTGGTTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                           <  1:643796/143‑1 (MQ=255)
                              aaaTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCACGCAGCAGCGGGGCTTTGCGACCAAACGAGAACGTaac                                                                                                          >  2:507012/1‑142 (MQ=255)
                                    gcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCg                                                                                                    >  2:421355/1‑143 (MQ=255)
                                     cgAATAAAACGTCTACTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCgg                                                                                                   <  1:388672/143‑1 (MQ=255)
                                       aaTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCt                                                                                                 <  2:370959/143‑1 (MQ=255)
                                              tGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAc                                                                                          <  2:169963/143‑1 (MQ=255)
                                                gcttcgCTGAAGCAGGATCTGCGGGCTGTGCTGTGTAAGGAATAACTCACGGCAGTCAGTTCAATTTACTAAAGGAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACaa                                                                                        <  2:103662/142‑1 (MQ=255)
                                                  ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACg                                                                                      >  1:208524/1‑143 (MQ=255)
                                                       tcaaGCTGGACCTGCAGGCAGTGCTGGGTAAGGACTAACTGACGACAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                 <  1:560589/141‑1 (MQ=255)
                                                       tGACGCAGGATCGGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGATGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                 <  1:52195/143‑1 (MQ=255)
                                                       tGAAGCAGGATCTGCAGGCATTGCTGGGTAACGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                 <  1:615297/143‑1 (MQ=255)
                                                       tGAAGCAGGATCGGCAGGCAAGGCTGGGAAAGGAAGAACGGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                 <  1:678899/143‑1 (MQ=255)
                                                            cAGGACCTTCAGGCATGGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                            <  2:300650/143‑1 (MQ=255)
                                                                 tCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAACAAACCCCGCAGCAGCGGGGTTTTTCTACCAGCCGAGAACTTAAGGATGAGGCTAAGTTACAACGTCTTTGATACCTTTa                                                                       >  1:56309/1‑143 (MQ=255)
                                                                  cTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAACTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                      <  1:914632/143‑1 (MQ=255)
                                                                  cTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTAAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                      <  1:586854/143‑1 (MQ=255)
                                                                      aGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCGACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTc                                                                  <  1:574165/143‑1 (MQ=255)
                                                                            tGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtctc                                                            <  1:914293/143‑1 (MQ=255)
                                                                               tGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTCCGATCCCTAc                                                         >  2:857823/1‑143 (MQ=255)
                                                                                   tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAACATAGCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                                                     <  2:359770/143‑1 (MQ=255)
                                                                                         aTAACTTATGAAAATCAGTACAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATcc                                               <  1:882363/143‑1 (MQ=255)
                                                                                          tAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCg                                              <  1:353536/143‑1 (MQ=255)
                                                                                            aCTGACGAAAGTCAGTACAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACCTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                            <  1:163267/143‑1 (MQ=255)
                                                                                                aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGCCGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCa                                        >  2:879568/1‑143 (MQ=255)
                                                                                                aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGCGCCa                                        >  1:198090/1‑143 (MQ=255)
                                                                                                aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGCTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACCATCCGGAgggc                                        >  1:199212/1‑140 (MQ=255)
                                                                                                   aaaGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTGTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGc                                     >  1:643761/1‑143 (MQ=255)
                                                                                                       gcAGTTCAAGCTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGGTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                 <  2:293862/142‑1 (MQ=255)
                                                                                                          gTTCACTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt                              <  2:556571/143‑1 (MQ=255)
                                                                                                           ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATc                             >  2:338346/1‑143 (MQ=255)
                                                                                                           ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATc                             >  1:631539/1‑143 (MQ=255)
                                                                                                              aagttACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGt                          <  1:143050/140‑1 (MQ=255)
                                                                                                              aaTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTGGATCAGt                          >  1:655545/1‑143 (MQ=255)
                                                                                                              aaTTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGACCCGGAGCCAGGCAGCCGATCCGt                          >  1:875268/1‑143 (MQ=255)
                                                                                                               aTTTACTAAAGGGAAAAAAAACCCCGCAGCAGCGGGGTGTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTg                         <  1:431502/143‑1 (MQ=255)
                                                                                                                 ttACTAAAGGCAAAAAAAACCACGCAGCTGCGGGGTTTTTCTTCCACACTAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                       <  1:904390/143‑1 (MQ=255)
                                                                                                                 ttACTAAACGCAAAAAAAGCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                       <  1:93877/143‑1 (MQ=255)
                                                                                                                     tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGCGCCAGGCAGTCGATCAGTgccgccc                   >  2:624664/1‑141 (MQ=255)
                                                                                                                       aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGCGCCAGGCAGTCGACCCGTGCAGCACGc                 >  2:764352/1‑143 (MQ=255)
                                                                                                                       aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGCCTCCGAGACCTTTAACTTCGATCTCTACGCGACGATCCGGCGCCAGGCAGTCGATCAGTGCCGCACGc                 >  1:420930/1‑143 (MQ=255)
                                                                                                                              aactacgcccGCAGCAGCGGGGTTTTTCTACCAGACGACAACTTAAGCCTGCGGCTGAGTCACAACGTCTTTGATACCTTTACCTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa          <  2:675192/136‑1 (MQ=255)
                                                                                                                              aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTcc          >  1:844061/1‑142 (MQ=255)
                                                                                                                               aaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAACAGTGCAGCACGCTGTTTCAc         >  1:437693/1‑143 (MQ=255)
                                                                                                                                  atccccGCGGCGCCGGAGGTCTTCTACCAGCCCAGAACTTACGCCTGCGGCTGAGTTACAGCCTCTGTGATCCCTTTAACTTCGATCTCTACGCGACGAACCGGATCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt      <  2:276260/141‑1 (MQ=255)
                                                                                                                                  aaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTCTTTCACGtt      >  1:264003/1‑143 (MQ=255)
                                                                                                                                   aCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAGCGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGCCAGTCGATCAGTACACCACGCGGTTCAACGttt     >  2:780503/1‑142 (MQ=255)
                                                                                                                                      ccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGATGTca  <  2:852679/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                       
CAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCA  >  NC_000913/1018845‑1019121

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