Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,202,356 G→A 11.1% A694T (GCA→ACA)  yehI → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,202,3560GA11.1% 111.5 / 5.0 45A694T (GCA→ACA) yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/3);  total (21/24)
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.19e-01

CCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCTTAAACAGGGTCCAATGGACATTTACCCGGCACGCATTC  >  NC_000913/2202221‑2202503
                                                                                                                                       |                                                                                                                                                   
ccTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCt                                                                                                                                        >  1:796317/1‑149 (MQ=255)
ccTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTAACAAACCAAGCGATAATAGCGAAGCCGCATTACTTGCCCt                                                                                                                                        >  1:200430/1‑149 (MQ=255)
  tGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGc                                                                                                                                      <  1:387707/149‑1 (MQ=255)
    ttGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGt                                                                                                                                    <  1:614466/149‑1 (MQ=255)
    ttGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTCATCTGGCTCGCCAGTGGATTCTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCACTACTTGCCCTGCgc                                                                                                                                    >  2:554458/1‑148 (MQ=255)
       cccACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCgaaaa                                                                                                                                 >  2:799101/1‑145 (MQ=255)
       cccACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTa                                                                                                                                 >  1:780203/1‑149 (MQ=255)
            gTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATGGCGAAGCCGCGTTACTTGCCCTGCGTTTACTGTa                                                                                                                            <  1:439807/149‑1 (MQ=255)
             tCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTAc                                                                                                                           <  2:684668/149‑1 (MQ=255)
                       tttCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGATTACTGTACGAACAGGGAc                                                                                                                 >  2:532721/1‑149 (MQ=255)
                           gcggcgCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTAGCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTCAGTTTACTGTACGAACAGGGACATgg                                                                                                             >  1:492324/1‑149 (MQ=255)
                                  aGCGTGATCTGGCTGGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTg                                                                                                      <  1:46525/149‑1 (MQ=255)
                                  aGCGTGATCTGGCTCGCCAGTGGATTTTACCATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGACGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTg                                                                                                      >  1:273052/1‑149 (MQ=255)
                                     gTGATCTGGCTCGCCAGTGGATCTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTa                                                                                                   <  1:380600/149‑1 (MQ=255)
                                     gTGATCGGGCCCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTa                                                                                                   <  1:147859/149‑1 (MQ=255)
                                             gcgcGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGt                                                                                           <  2:180861/146‑1 (MQ=255)
                                                   cAGTGGATTTTACAATGGCCGGAACATACGGCTCGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGCAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAACCCGTGGCaaa                                                                                     >  2:3057/1‑149 (MQ=255)
                                                   cAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGCTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCaaa                                                                                     >  1:331262/1‑149 (MQ=255)
                                                        gATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCACATTACTTGCCCTGCGCTTACTATACGCACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCt                                                                                >  1:560028/1‑149 (MQ=255)
                                                        gATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGAAAATAGCGAAGCCACATTACTTGCCCTGCGTTTACTATACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCt                                                                                >  1:560128/1‑149 (MQ=255)
                                                          ttttACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTgg                                                                              <  2:375249/149‑1 (MQ=255)
                                                          ttttACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTgg                                                                              <  1:553903/149‑1 (MQ=255)
                                                               cAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCg                                                                         <  2:149642/149‑1 (MQ=255)
                                                                aaTGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTCCAAACCGTGGCAAACCGCTGGCAGCGt                                                                        >  2:163803/1‑149 (MQ=255)
                                                                     ccGGAACCTACGGCTAGTGCACTTACCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGa                                                                   <  2:79327/149‑1 (MQ=255)
                                                                     ccGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAATCAGTGGCAAACCGCTGGCAGCGTACAGa                                                                   <  2:306248/149‑1 (MQ=255)
                                                                         aaCATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCCGCGTACAGATGTa                                                               >  2:274626/1‑149 (MQ=255)
                                                                              aCGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTc                                                          >  1:125098/1‑149 (MQ=255)
                                                                              aCGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGCACAGATGTATGGTc                                                          >  1:802371/1‑149 (MQ=255)
                                                                                       gCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCACATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGa                                                 <  1:721452/149‑1 (MQ=255)
                                                                                          cTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTTGCAGCGTACAGATGTATGGTCTGCCCTGGAGCa                                              <  2:568155/149‑1 (MQ=255)
                                                                                              tCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGTCGTATTACTTGCCCGGCGTTTACTGTACGAACAGGGACAGGGCGAATTGCTAAAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTg                                          <  2:136370/149‑1 (MQ=255)
                                                                                               cccTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGc                                         <  1:799101/149‑1 (MQ=255)
                                                                                               cccTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGc                                         >  1:740955/1‑149 (MQ=255)
                                                                                                ccTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCt                                        >  1:684525/1‑149 (MQ=255)
                                                                                                ccTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCt                                        >  1:684498/1‑149 (MQ=255)
                                                                                                 ctcttctcttTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCtt                                       <  1:539072/149‑1 (MQ=255)
                                                                                                 ctctTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCtt                                       <  1:360826/149‑1 (MQ=255)
                                                                                                  tctTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGAGGTATGGTCAGCGCTGGTGCAGTTGCGTa                                      >  2:723182/1‑149 (MQ=255)
                                                                                                           tACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCTTAAACAGGGTc                             <  2:756948/149‑1 (MQ=255)
                                                                                                                     aGCGATAATGGCGAAGCCGCATTACTTGCCCTGCGTTTACCGTACGAACAGGGACATGGCGAATTGCTACAGACCGTGGCAAACCGCTGGCCGCGTACAGATGTATGGTCTGCCCGTGAGCAGCTGCTTAAACAGGGTCCAATGGAccg                   >  2:442557/1‑147 (MQ=255)
                                                                                                                                ggAAGCCACATTACTTGCCCTATTCGGACTACACGGACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCTTAAACAGGGTCCAATGGACATTTACCCGGCAc        <  2:560028/148‑1 (MQ=255)
                                                                                                                                cGAAGCCACATTACTTGCCCTGCGTTTACTATACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCTTAAACAGGGTCCAATGGACATTTACCCGGCAc        <  2:560128/149‑1 (MQ=255)
                                                                                                                                  aaGCCGCATTACTTGCCCTGCGGTTACCGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCTTAAACAGGGTCCAATGGACATTTACCCGGCACGc      <  1:52767/149‑1 (MQ=255)
                                                                                                                                      cGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCTTAAACAGGGTCCAATGGACATTTACCCGGCACGCATTc  <  1:795284/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
CCTGTTGCCCACGTCTGGCGACGTTTCGCGGCGCAGCGTGATCTGGCTCGCCAGTGGATTTTACAATGGCCGGAACATACGGCTAGTGCACTTATCCCTCTTGTCTTTACCAAACCCAGCGATAATAGCGAAGCCGCATTACTTGCCCTGCGTTTACTGTACGAACAGGGACATGGCGAATTGCTACAAACCGTGGCAAACCGCTGGCAGCGTACAGATGTATGGTCTGCCCTGGAGCAGTTGCTTAAACAGGGTCCAATGGACATTTACCCGGCACGCATTC  >  NC_000913/2202221‑2202503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑