Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,678 T→A 10.0% intergenic (+32/‑334) 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,6780TA10.0% 167.7 / 9.4 50intergenic (+32/‑334)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (19/26);  new base A (2/3);  total (21/29)
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.95e-01

GTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTT  >  NC_000913/2462533‑2462823
                                                                                                                                                 |                                                                                                                                                  
gTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGc                                                                                                                                                <  4:723255/149‑1 (MQ=255)
 ttGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCa                                                                                                                                               <  3:147540/149‑1 (MQ=255)
  tGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAg                                                                                                                                              <  2:60643/149‑1 (MQ=255)
  tGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAg                                                                                                                                              <  2:682895/149‑1 (MQ=255)
    gtgtTTCTTATATGCACGGTCAGAGAGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGGTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                            <  3:560443/149‑1 (MQ=255)
      gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGTTAAAGTCACCTGCTATGCAGGTGa                                                                                                                                          <  3:351740/149‑1 (MQ=255)
      gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                          <  4:792944/149‑1 (MQ=255)
      gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                          <  2:766883/149‑1 (MQ=255)
      gtTTCTTAGATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGa                                                                                                                                          <  1:438030/149‑1 (MQ=255)
         tCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAATCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACtt                                                                                                                                       >  4:176459/1‑149 (MQ=255)
                                 aaaaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagt                                                                                                               >  3:509268/1‑149 (MQ=255)
                                   aaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagtag                                                                                                             >  3:458623/1‑149 (MQ=255)
                                                       aGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcctcc                                                                                         <  1:416220/149‑1 (MQ=255)
                                                        gTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcctcct                                                                                        <  4:509268/149‑1 (MQ=255)
                                                               tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                                 >  2:471585/1‑149 (MQ=255)
                                                                      gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGaaa                                                                                                                 <  4:123083/110‑1 (MQ=255)
                                                                      gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGaaa                                                                                                                 >  3:123083/1‑110 (MQ=255)
                                                                             cTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                   <  4:458623/149‑1 (MQ=255)
                                                                                            tAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                    >  2:600484/1‑149 (MQ=255)
                                                                                            tAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                    >  2:271383/1‑149 (MQ=255)
                                                                                             aaCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                                 <  2:412982/119‑1 (MQ=255)
                                                                                             aaCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                                 >  1:412982/1‑119 (MQ=255)
                                                                                                tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                >  4:720073/1‑149 (MQ=255)
                                                                                                         gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                       <  1:271383/149‑1 (MQ=255)
                                                                                                           gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                   <  2:57917/119‑1 (MQ=255)
                                                                                                           gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                   >  1:57917/1‑119 (MQ=255)
                                                                                                           gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTgg                                                                 <  4:34773/121‑1 (MQ=255)
                                                                                                           gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTgg                                                                 >  3:34773/1‑121 (MQ=255)
                                                                                                               tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATc                                              <  2:331931/136‑1 (MQ=255)
                                                                                                               tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATc                                              >  1:331931/1‑136 (MQ=255)
                                                                                                                gATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTcc                                                                                                                                >  2:531819/1‑53 (MQ=255)
                                                                                                                gATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTcc                                                                                                                                <  1:531819/53‑1 (MQ=255)
                                                                                                                 aTAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCca                               <  3:176459/149‑1 (MQ=255)
                                                                                                                     cgcgTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGa                           <  3:720073/149‑1 (MQ=255)
                                                                                                                         ttCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                                 <  1:422804/139‑1 (MQ=255)
                                                                                                                         ttCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                                 >  2:422804/1‑139 (MQ=255)
                                                                                                                                  aTAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAggg              <  1:471585/149‑1 (MQ=255)
                                                                                                                                     aaGTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATaa                                         <  1:350636/119‑1 (MQ=255)
                                                                                                                                     aaGTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATaa                                         >  2:611218/1‑119 (MQ=255)
                                                                                                                                     aaGTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATaa                                         <  1:611218/119‑1 (MQ=255)
                                                                                                                                     aaGTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATaa                                         >  2:350636/1‑119 (MQ=255)
                                                                                                                                      aGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg          <  1:600484/149‑1 (MQ=255)
                                                                                                                                       gTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAAt                                               <  2:380467/111‑1 (MQ=255)
                                                                                                                                       gTCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAAt                                               >  1:380467/1‑111 (MQ=255)
                                                                                                                                          aCCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat      >  2:397939/1‑149 (MQ=255)
                                                                                                                                            cTGCATAGCAGGTGACTTTAACTCCGCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaa                                                 <  4:575625/104‑1 (MQ=255)
                                                                                                                                            cTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaa                                                 >  3:575625/1‑104 (MQ=255)
                                                                                                                                              gCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcaca                             >  2:211082/1‑122 (MQ=255)
                                                                                                                                              gCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcaca                             <  1:211082/122‑1 (MQ=255)
                                                                                                                                              gCATAGCAGGTGACTTTAACTCCCCCACCTCTCCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCGTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAttt  >  4:490561/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTT  >  NC_000913/2462533‑2462823

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