Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,226,509 C→T *254* (TAG→TAA yohF ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,226,5090CT97.7% 129.5 / ‑5.4 43*254* (TAG→TAAyohFputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/19);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCA  >  NC_000913/2226391‑2226620
                                                                                                                      |                                                                                                               
tGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTgg                                                                                             >  2:1161180/1‑139 (MQ=255)
 ggATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                            <  1:1161180/139‑1 (MQ=255)
       ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                      >  1:123543/1‑139 (MQ=255)
       ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                      >  2:137002/1‑139 (MQ=255)
            tCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                             <  2:752328/111‑1 (MQ=255)
            tCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                             >  1:752328/1‑111 (MQ=255)
            tCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACAt                                                                                 >  2:303430/1‑139 (MQ=255)
                    cATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccgc                                                                         <  1:137002/139‑1 (MQ=255)
                              ccGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTAt                                                               >  2:1072101/1‑139 (MQ=255)
                                         tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCCCATCCACTATCAACGACTGcc                                                    >  1:639178/1‑138 (MQ=255)
                                         tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCCCATCCACTATCAACGACTGcc                                                    <  2:639178/138‑1 (MQ=255)
                                              gATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGa                                                         <  2:342520/129‑1 (MQ=255)
                                              gATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGa                                                         >  1:342520/1‑129 (MQ=255)
                                                    acgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCAt                                                                      >  2:457060/1‑110 (MQ=255)
                                                    acgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCAt                                                                      <  1:457060/110‑1 (MQ=255)
                                                     cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                   >  1:238810/1‑112 (MQ=255)
                                                     cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                   <  2:238810/112‑1 (MQ=255)
                                                         tgtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                   <  2:315987/108‑1 (MQ=255)
                                                         tgtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                   >  1:315987/1‑108 (MQ=255)
                                                              gTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCt                               <  2:40210/139‑1 (MQ=255)
                                                               tCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATc                                                                     >  2:992875/1‑100 (MQ=255)
                                                               tCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATc                                                                     <  1:992875/100‑1 (MQ=255)
                                                               tCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTc                              >  1:696134/1‑139 (MQ=255)
                                                                   ttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                   <  1:513067/98‑1 (MQ=255)
                                                                   ttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                   >  2:513067/1‑98 (MQ=255)
                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                          <  1:744490/82‑1 (MQ=255)
                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                          >  2:744490/1‑82 (MQ=255)
                                                                            aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGCTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAGCGACTGCCCGGTGGTGTAATTTGCGCCCTCAGAACAAAGCaac                 >  2:677850/1‑139 (MQ=255)
                                                                                         ccAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACcaggcagg    <  1:877861/139‑1 (MQ=255)
                                                                                           aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGa                                                         >  1:887833/1‑84 (MQ=255)
                                                                                           aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGa                                                         <  2:887833/84‑1 (MQ=255)
                                                                                                       aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAg                                                                             <  1:886264/52‑1 (MQ=255)
                                                                                                       aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAg                                                                             >  2:886264/1‑52 (MQ=255)
                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGc                                                                                       >  1:783904/1‑40 (MQ=255)
                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGc                                                                                       <  2:783904/40‑1 (MQ=255)
                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggtggt                                             <  2:977930/82‑1 (MQ=255)
                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggtggt                                             >  1:977930/1‑82 (MQ=255)
                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAg         <  1:801134/118‑1 (MQ=255)
                                                                                                         gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAg         >  2:801134/1‑118 (MQ=255)
                                                                                                              tCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTgg    >  2:481569/1‑118 (MQ=255)
                                                                                                              tCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTgg    <  1:481569/118‑1 (MQ=255)
                                                                                                               cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa  >  2:674155/1‑119 (MQ=255)
                                                                                                               cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa  <  1:674155/119‑1 (MQ=255)
                                                                                                                   gcgtTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc       <  1:232954/110‑1 (MQ=255)
                                                                                                                   gcgtTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc       >  2:232954/1‑110 (MQ=255)
                                                                                                                      |                                                                                                               
TGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCA  >  NC_000913/2226391‑2226620

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: