Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,677 A→T 23.0% intergenic (+31/‑335) fadL → / → yfdF long‑chain fatty acid outer membrane transporter/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,6770AT23.0% 49.7 / 13.8 28intergenic (+31/‑335)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (12/8);  new base T (2/4);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

CGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTT  >  NC_000913/2462529‑2462824
                                                                                                                                                    |                                                                                                                                                    
cGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCa                                                                                                                                                      <  4:449118/149‑1 (MQ=255)
    gTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑Gc                                                                                                                                                 <  4:781085/149‑1 (MQ=255)
                       aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                                >  4:351379/1‑131 (MQ=255)
                       aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                                <  3:351379/131‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGa                                                                                                                                           <  4:802560/107‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGa                                                                                                                                           >  3:802560/1‑107 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                           >  1:664741/1‑107 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                           <  2:664741/107‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACt                                                                                                                                         <  2:884977/109‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACt                                                                                                                                         <  2:1160221/109‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACt                                                                                                                                         >  1:1160221/1‑109 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACt                                                                                                                                         >  1:884977/1‑109 (MQ=255)
                                                      aTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGAc                                                                                                                                          <  1:1415537/106‑1 (MQ=255)
                                                      aTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGAc                                                                                                                                          >  2:1415537/1‑106 (MQ=255)
                                                         ccAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGt                                                                                                                                             <  1:1232990/100‑1 (MQ=255)
                                                                   tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGa                                                                                                                                           <  2:170030/92‑1 (MQ=255)
                                                                   tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGa                                                                                                                                           >  1:170030/1‑92 (MQ=39)
                                                                     tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                             <  3:275882/88‑1 (MQ=255)
                                                                     tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                             >  4:275882/1‑88 (MQ=255)
                                                                      gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAg                                                                                                                                               >  3:209563/1‑85 (MQ=255)
                                                                      gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAg                                                                                                                                               <  4:209563/85‑1 (MQ=255)
                                                                           tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                           <  3:190841/84‑1 (MQ=255)
                                                                           tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                           >  4:190841/1‑84 (MQ=255)
                                                                               gCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCa                                                                                                                                                      >  3:1238030/1‑70 (MQ=255)
                                                                               gCCTGGCTGTTCGGTACTAACTTTAACTACGCGGTCTGATAACGCGTTCGCCTGGATCAAGTCACCTGCa                                                                                                                                                      <  4:1238030/70‑1 (MQ=255)
                                                                                                                                            tCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcac                               >  4:166003/1‑127 (MQ=255)
                                                                                                                                              aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat       >  1:502410/1‑149 (MQ=255)
                                                                                                                                                cTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat     >  1:517331/1‑149 (MQ=255)
                                                                                                                                                cTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat     >  2:1257675/1‑149 (MQ=255)
                                                                                                                                                  gCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAttt   >  1:1113078/1‑149 (MQ=255)
                                                                                                                                                   cATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAtttt  >  2:1127104/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTT  >  NC_000913/2462529‑2462824

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