Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,462,889 (C)9→10 intergenic (+243/‑123) fadL → / → yfdF long‑chain fatty acid outer membrane channel/bacteriophage T2 receptor/protein YfdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,8801.C78.6% 76.0 / 15.4 28intergenic (+234/‑132)fadL/yfdFlong‑chain fatty acid outer membrane channel/bacteriophage T2 receptor/protein YfdF
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base C (9/13);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATG  >  NC_000913/2462743‑2463001
                                                                                                                                          |                                                                                                                          
ttGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑cc                                                                                                                         >  2:370368/1‑140 (MQ=255)
      aCAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑cccccccc                                                                                                                   <  2:2236/140‑1 (MQ=255)
      aCAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑cccccccc                                                                                                                   <  2:64190/140‑1 (MQ=255)
        aGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑cccccccccc                                                                                                                 >  1:280641/1‑139 (MQ=255)
                                    cAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAAt                                                                                      <  1:320417/140‑1 (MQ=255)
                                       tAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCCAGCCTATTGTCCGCCCATACCCCCAAATTa                                                                                    >  2:326557/1‑140 (MQ=255)
                                          gCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt                                                                               <  1:350786/140‑1 (MQ=255)
                                                tCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGtt                                                                          <  1:104134/140‑1 (MQ=255)
                                                tCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGtt                                                                          <  2:131951/140‑1 (MQ=255)
                                                     aGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCt                                                                        >  2:48265/1‑137 (MQ=255)
                                                     aGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCt                                                                        <  1:48265/137‑1 (MQ=255)
                                                        cAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCtaaataa                                                                  >  2:29316/1‑140 (MQ=255)
                                                           aCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTTTTGGGGCCCCATACACTCAATTTAAAGATAGGTTCTAATTAATTg                                                               >  2:287495/1‑140 (MQ=255)
                                                           aCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATg                                                               >  2:93578/1‑140 (MQ=255)
                                                           aCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGa                                                              >  2:157920/1‑140 (MQ=255)
                                                            ccTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGa                                                              >  1:351716/1‑140 (MQ=255)
                                                            ccTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAg                                                             >  2:310089/1‑140 (MQ=255)
                                                                 aGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGtt                                                         <  1:370368/140‑1 (MQ=255)
                                                                 aGGGAAACGCATATTTTCAGGCCTTCTGCTCGCCCTTCTCAACTAGTTTACTTTTCTGCATTTCCAGGATATTCCCCCCCCCCCGGCTATTGCGGGCCCATCCTCACAAATTAAAGATGGGTTCTGAATAAATGAACGtt                                                         >  2:284389/1‑140 (MQ=255)
                                                                      aaGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGa                                                   <  1:157920/140‑1 (MQ=255)
                                                                       aGGCATATTTTCCGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATTCTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGa                                                   <  1:287495/140‑1 (MQ=255)
                                                                       aGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGa                                                   <  2:280641/140‑1 (MQ=255)
                                                                       aGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGa                                                   <  1:93578/140‑1 (MQ=255)
                                                                                      cGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCAtt                                     <  2:389943/140‑1 (MQ=255)
                                                                                        ttCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg                                  <  2:102417/140‑1 (MQ=255)
                                                                                        ttCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg                                  <  1:380954/140‑1 (MQ=255)
                                                                                        ttCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg                                  >  1:299744/1‑140 (MQ=255)
                                                                                           tGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTaa                               <  1:29316/140‑1 (MQ=255)
                                                                                             cTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAatat                            <  1:310089/140‑1 (MQ=255)
                                                                                                gTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAatatat                          <  2:351716/140‑1 (MQ=255)
                                                                                                          aaGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACaa               <  1:218591/140‑1 (MQ=255)
                                                                                                               tttACTTTTCTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAaatca           >  2:183614/1‑140 (MQ=255)
                                                                                                                        cTGCATTTCCAGGATACTCCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATg  >  1:360393/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                          
TTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑‑CCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATG  >  NC_000913/2462743‑2463001

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