Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,793,668 C→A 11.4% E148* (GAA→TAA) ‡ btuE ← glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,6680CA11.4% 113.4 / 6.7 35E148* (GAA→TAA) ‡btuEglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base C (16/15);  new base A (2/2);  total (18/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 = 9.36e-01

TTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAA  >  NC_000913/1793558‑1793814
                                                                                                              |                                                                                                                                                  
ttattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTcc                                                                                                                                            <  3:160656/119‑1 (MQ=255)
ttattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTcc                                                                                                                                            >  4:160656/1‑119 (MQ=255)
 tattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTc                                                                                                                                             <  1:21363/117‑1 (MQ=255)
 tattTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTc                                                                                                                                             >  2:21363/1‑117 (MQ=255)
                                tCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATatcat                                                                                                                                <  3:120313/99‑1 (MQ=255)
                                tCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATatcat                                                                                                                                >  4:120313/1‑99 (MQ=255)
                                 cacaATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCg                                                                              >  4:418147/1‑148 (MQ=255)
                                     aTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaaa                                                                                                                                                 >  1:1037665/1‑77 (MQ=255)
                                     aTGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaaa                                                                                                                                                 <  2:1037665/77‑1 (MQ=255)
                                       gggATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGc                                                                                                                     <  1:822970/103‑1 (MQ=255)
                                       gggATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGc                                                                                                                     >  2:822970/1‑103 (MQ=255)
                                       gggATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctg                                                                                                                                          <  2:207468/82‑3 (MQ=255)
                                       gggATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTGAAAAATTcctg                                                                                                                                          >  1:207468/1‑80 (MQ=255)
                                         gATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaaa                                                                                                                                                 <  3:399172/73‑1 (MQ=255)
                                         gATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCaaa                                                                                                                                                 >  4:399172/1‑73 (MQ=255)
                                                           ccGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCAt                                                                                      >  4:266273/1‑114 (MQ=255)
                                                           ccGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCAt                                                                                      >  4:313291/1‑114 (MQ=255)
                                                           ccGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCAt                                                                                      <  3:266273/114‑1 (MQ=255)
                                                           ccGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCAt                                                                                      <  3:313291/114‑1 (MQ=255)
                                                                        cTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctc                                                                           <  2:775968/112‑1 (MQ=255)
                                                                        cTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGctc                                                                           >  1:775968/1‑112 (MQ=255)
                                                                                  tttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTg                                                                                                   >  3:570945/1‑78 (MQ=255)
                                                                                  tttCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTg                                                                                                   <  4:570945/78‑1 (MQ=255)
                                                                                        tCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCg                                                         <  3:737287/114‑1 (MQ=255)
                                                                                        tCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCg                                                         >  4:737287/1‑114 (MQ=255)
                                                                                         cccTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGgcg                                                                 <  2:562081/105‑1 (MQ=255)
                                                                                         cccTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGgcg                                                                 >  1:562081/1‑105 (MQ=255)
                                                                                           cTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAAt                                                                                 <  3:423066/87‑1 (MQ=255)
                                                                                           cTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAAt                                                                                 >  4:423066/1‑87 (MQ=255)
                                                                                            tGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATg                   >  3:1082474/1‑148 (MQ=255)
                                                                                                       aaTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGcc                                                                                                          <  3:314024/50‑1 (MQ=255)
                                                                                                       aaTTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGcc                                                                                                          >  4:314024/1‑50 (MQ=255)
                                                                                                         tttttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAttttt                                >  1:582991/1‑122 (MQ=255)
                                                                                                         tttttCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAAttttt                                <  2:582991/122‑1 (MQ=255)
                                                                                                             tCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTaa  <  4:1082474/148‑1 (MQ=255)
                                                                                                              |                                                                                                                                                  
TTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTTCAAAATTCCATAAAATATCATCCGGGTACAGCGGTGCACGGCCTTTGCTGACCATACGGGCATAGAATCCGCTCTCTTCCGGCGCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAA  >  NC_000913/1793558‑1793814

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: