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

CGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTT  >  NC_000913/2462529‑2462824
                                                                                                                                                    |                                                                                                                                                    
cGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCa                                                                                                                                                      <  4:449118/149‑1 (MQ=255)
    gTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGc                                                                                                                                                 <  4:781085/149‑1 (MQ=255)
                       aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCa                                                                                                                                                >  4:351379/1‑131 (MQ=255)
                       aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCa                                                                                                                                                <  3:351379/131‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa                                                                                                                                           <  4:802560/107‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa                                                                                                                                           >  3:802560/1‑107 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                           >  1:664741/1‑107 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                           <  2:664741/107‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt                                                                                                                                         <  2:884977/109‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt                                                                                                                                         <  2:1160221/109‑1 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt                                                                                                                                         >  1:1160221/1‑109 (MQ=255)
                                                    ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACt                                                                                                                                         >  1:884977/1‑109 (MQ=255)
                                                      aTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGAc                                                                                                                                          <  1:1415537/106‑1 (MQ=255)
                                                      aTACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGAc                                                                                                                                          >  2:1415537/1‑106 (MQ=255)
                                                         ccAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGt                                                                                                                                             <  1:1232990/100‑1 (MQ=255)
                                                                   tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa                                                                                                                                           <  2:170030/92‑1 (MQ=255)
                                                                   tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa                                                                                                                                           >  1:170030/1‑92 (MQ=39)
                                                                     tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                             <  3:275882/88‑1 (MQ=255)
                                                                     tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                             >  4:275882/1‑88 (MQ=255)
                                                                      gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAg                                                                                                                                               >  3:209563/1‑85 (MQ=255)
                                                                      gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAg                                                                                                                                               <  4:209563/85‑1 (MQ=255)
                                                                           tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                           <  3:190841/84‑1 (MQ=255)
                                                                           tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                           >  4:190841/1‑84 (MQ=255)
                                                                               gCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCa                                                                                                                                                      >  3:1238030/1‑70 (MQ=255)
                                                                               gCCTGGCTGTTCGGTACTAACTTTAACTACGCGGTCTGATAACGCGTTCGCCTGGATCAAGTCACCTGCa                                                                                                                                                      <  4:1238030/70‑1 (MQ=255)
                                                                                                                                            tCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcac                               >  4:166003/1‑127 (MQ=255)
                                                                                                                                              aCCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat       >  1:502410/1‑149 (MQ=255)
                                                                                                                                                cTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat     >  1:517331/1‑149 (MQ=255)
                                                                                                                                                cTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat     >  2:1257675/1‑149 (MQ=255)
                                                                                                                                                  gCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAttt   >  1:1113078/1‑149 (MQ=255)
                                                                                                                                                   cATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAtttt  >  2:1127104/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTT  >  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: