Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,018,976 (A)8→9 intergenic (+39/+37) 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,9681.A92.9% 135.3 / 2.7 42intergenic (+31/+45)matP/ompAmacrodomain Ter protein/outer membrane porin A
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (18/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

GTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAT  >  NC_000913/1018842‑1019093
                                                                                                                               |                                                                                                                              
gTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAacccc                                                                                                                  >  1:1008437/1‑139 (MQ=255)
  tCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCgc                                                                                                                 <  2:1008437/140‑1 (MQ=255)
      cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCgcagca                                                                                                             >  2:1062323/1‑140 (MQ=255)
       tGATTAAAGATGCGAAAAAAAAAGAAAAATACGCAAAAAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCGGGGAAAGAAAAAACTGAAGAAAGTCAGTTAATTTTACAAAAGGAAAAAAAAACCCCCGcagcag                                                                                                            <  2:971377/140‑1 (MQ=255)
        gATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCagcagc                                                                                                           >  1:539784/1‑140 (MQ=255)
                     gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAc                                                                                              <  2:539784/140‑1 (MQ=255)
                     gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCACTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAc                                                                                              <  2:269816/140‑1 (MQ=255)
                       aaaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCa                                                                                            <  1:7478/140‑1 (MQ=255)
                        aaCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAccgg                                                                                           >  1:694501/1‑138 (MQ=255)
                            aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                       <  1:852160/140‑1 (MQ=255)
                             agagAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                       <  1:117903/139‑1 (MQ=255)
                             agagAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                       >  2:117903/1‑139 (MQ=255)
                                         gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACtt                                                                                 >  1:406589/1‑133 (MQ=255)
                                         gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACtt                                                                                 <  2:406589/133‑1 (MQ=255)
                                         gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTg                                                                          >  2:7452/1‑140 (MQ=255)
                                            tAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCgg                                                                       <  2:971676/140‑1 (MQ=255)
                                                      tctcTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCGGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                             <  2:694501/140‑1 (MQ=255)
                                                      tctcTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                             <  1:7452/140‑1 (MQ=255)
                                                      tctcTGAAGCAGGATCCGCAGGCATTGCTGGGTAAGGAAAAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCCGGCTTCTTCTACCAGACGAGAACTTAAGCCTTCGGCTGAGTTACa                                                             <  2:1142295/140‑1 (MQ=255)
                                                       ctctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACaa                                                            >  2:955531/1‑140 (MQ=255)
                                                          tGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGt                                                         >  1:347333/1‑140 (MQ=255)
                                                             aGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCtt                                                      >  2:1070278/1‑140 (MQ=255)
                                                               cAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTg                                                    <  2:473538/140‑1 (MQ=255)
                                                                 ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGAt                                                  <  1:1070278/140‑1 (MQ=255)
                                                                 ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAAc                                                           >  1:147957/1‑130 (MQ=255)
                                                                 ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAAc                                                           <  2:147957/130‑1 (MQ=255)
                                                                             aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                       >  2:871184/1‑139 (MQ=255)
                                                                             aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                       >  2:387834/1‑139 (MQ=255)
                                                                             aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                       <  1:871184/139‑1 (MQ=255)
                                                                             aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                       <  1:387834/139‑1 (MQ=255)
                                                                             aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTc                                      >  1:781055/1‑140 (MQ=255)
                                                                               tGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGGCGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGa                                    >  1:538713/1‑140 (MQ=255)
                                                                                     gTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTa                              >  1:207801/1‑140 (MQ=255)
                                                                                      tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAc                             >  1:659446/1‑140 (MQ=255)
                                                                                          gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                         <  2:659446/140‑1 (MQ=255)
                                                                                               aCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGGCGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATc                    <  2:538713/140‑1 (MQ=255)
                                                                                               aCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                   <  2:121704/125‑1 (MQ=255)
                                                                                               aCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                   >  1:121704/1‑125 (MQ=255)
                                                                                                cTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                         <  1:329412/134‑1 (MQ=255)
                                                                                                cTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                         >  2:329412/1‑134 (MQ=255)
                                                                                                             gTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCa            >  2:805469/1‑134 (MQ=255)
                                                                                                             gTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCa            <  1:805469/134‑1 (MQ=255)
                                                                                                                 aaTTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt  <  2:599136/140‑1 (MQ=255)
                                                                                                                 aaTTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt  <  2:207801/140‑1 (MQ=255)
                                                                                                                               |                                                                                                                              
GTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAT  >  NC_000913/1018842‑1019093

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