Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,462,889 (C)9→8 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,8810C.94.1% 55.6 / ‑0.0 18intergenic (+235/‑131)fadL/yfdFlong‑chain fatty acid outer membrane channel/bacteriophage T2 receptor/protein YfdF
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (6/10);  minor base T (0/1);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTG  >  NC_000913/2462749‑2462982
                                                                                                                                    |                                                                                                     
aCAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTcccccccc                                                                                                <  1:154092/140‑1 (MQ=255)
  aGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTcccccccctg                                                                                              >  1:24030/1‑138 (MQ=255)
                 tGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCa                                                                              >  2:196815/1‑140 (MQ=255)
                      tttCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATacac                                                                         <  1:100224/140‑1 (MQ=255)
                                  aaGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGa                                                             <  1:304408/140‑1 (MQ=255)
                                   aGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt                                                            <  1:199202/140‑1 (MQ=255)
                                   aGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATATTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGAt                                                            <  1:178559/140‑1 (MQ=255)
                                          tCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCt                                                     >  2:81144/1‑140 (MQ=255)
                                                  cAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTaaataaat                                             >  2:90955/1‑140 (MQ=255)
                                                                 aGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATa                              <  1:90955/140‑1 (MQ=255)
                                                                 aGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATa                              >  1:17810/1‑140 (MQ=255)
                                                                  ggCATATTTTCAGGCGTTCTGCTCGTACTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAg                             <  1:248558/140‑1 (MQ=255)
                                                                         tttaaGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTAttt                      >  2:73347/5‑140 (MQ=255)
                                                                                cGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg               <  1:327831/140‑1 (MQ=255)
                                                                                cGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg               <  1:196815/140‑1 (MQ=255)
                                                                                  ttCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTa             <  2:17810/140‑1 (MQ=255)
                                                                                          gTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAatatatat     <  2:234705/140‑1 (MQ=255)
                                                                                            ccTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCC‑‑GGCCATTTGGGCCCCATAAACTCAAATTAAAAATTGGGTCCAAATAAATTAGCGGTTTTTGATTATCTATTTTCTTTGGGAATATATAtttt  >  2:351491/1‑139 (MQ=255)
                                                                                               tctcAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg               >  1:197290/1‑125 (MQ=255)
                                                                                               tctcAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg               <  2:197290/125‑1 (MQ=255)
                                                                                                                                    |                                                                                                     
ACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTG  >  NC_000913/2462749‑2462982

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