Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,753,986 A→T 12.9% intergenic (+40/+22) bamE → / ← ratB lipoprotein component of BamABCDE OM biogenesis complex/UPF0125 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,753,9860AT12.9% 89.0 / 6.7 31intergenic (+40/+22)bamE/ratBlipoprotein component of BamABCDE OM biogenesis complex/UPF0125 family protein
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base T (2/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTA  >  NC_000913/2753838‑2754133
                                                                                                                                                    |                                                                                                                                                   
ccAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCaa                                                                                                                                                     <  3:26532/149‑1 (MQ=255)
   gCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGATCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGa                                                                                                                                                  <  1:1026215/149‑1 (MQ=255)
   gCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                    <  1:205439/147‑1 (MQ=255)
   gCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                    >  2:205439/1‑147 (MQ=255)
     aaCCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGc                                                                                                                                                <  2:927409/149‑1 (MQ=255)
                               cAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCAcc                                                                                                                                             >  4:624894/1‑126 (MQ=255)
                               cAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCAcc                                                                                                                                             <  3:624894/126‑1 (MQ=255)
                                 aaCGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGc                                                                                                                                                <  4:534403/121‑1 (MQ=255)
                                 aaCGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGc                                                                                                                                                >  3:534403/1‑121 (MQ=255)
                                          gctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                    >  1:1772137/1‑108 (MQ=255)
                                          gctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                    <  2:1772137/108‑1 (MQ=255)
                                                 tttAACAGTAGCGGTGTGTTGACCAATATTGATAACAAATCTGAGCTGAGTGGTAACTAATAGTAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGACTGTTATTTATTCGCTAATTTTACTGCTCGt                                                                                                    >  3:1585745/1‑149 (MQ=255)
                                                 tttAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttt                                                                                                                                         <  2:505197/112‑1 (MQ=255)
                                                 tttAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttt                                                                                                                                         >  1:505197/1‑112 (MQ=255)
                                                                   ttGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                    >  4:163608/1‑83 (MQ=255)
                                                                   ttGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                    <  3:163608/83‑1 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATgaga                                                                   <  4:1585745/149‑1 (MQ=255)
                                                                                         cTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTg                                                                                              <  1:1381227/115‑1 (MQ=255)
                                                                                         cTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTg                                                                                              >  2:1381227/1‑115 (MQ=255)
                                                                                                     gTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATg                                                                                                                                                   <  4:1208897/50‑1 (MQ=255)
                                                                                                     gTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATg                                                                                                                                                   >  3:1208897/1‑50 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                  >  2:1757469/1‑106 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                  <  1:1757469/106‑1 (MQ=255)
                                                                                                                     ttGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTc                                                                                                           >  3:927739/1‑74 (MQ=255)
                                                                                                                     ttGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTc                                                                                                           <  4:927739/74‑1 (MQ=255)
                                                                                                                                       ccgatcttGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGc              <  3:1663554/142‑1 (MQ=255)
                                                                                                                                         aaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATaaa                                                         >  3:589153/1‑104 (MQ=255)
                                                                                                                                         aaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATaaa                                                         <  4:589153/104‑1 (MQ=255)
                                                                                                                                               gCTCAATGAGCACCTTTTTTCCGTCTGTTATTTATTGGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCTCCATCATGCCCACTATTGGTTAGTTTTGCCAGatcgg      >  4:1663554/1‑145 (MQ=255)
                                                                                                                                                 tCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGcc             >  2:687949/1‑140 (MQ=255)
                                                                                                                                                 tCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGcc             <  1:687949/140‑1 (MQ=255)
                                                                                                                                                 tCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTg    >  2:1311085/1‑149 (MQ=255)
                                                                                                                                                   aaTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGa                                                             <  3:727847/90‑1 (MQ=255)
                                                                                                                                                   aaTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGa                                                             >  4:727847/1‑90 (MQ=255)
                                                                                                                                                   aaTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTa  >  2:981232/1‑149 (MQ=255)
                                                                                                                                                    aTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCcat                                       >  2:147245/1‑111 (MQ=255)
                                                                                                                                                    aTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCcat                                       <  1:147245/111‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTA  >  NC_000913/2753838‑2754133

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