Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 849,783 (A)6→7 coding (909/912 nt)
coding (1/996 nt)
yvcJ →
yvcK →
ATP‑binding protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730849,7771.A89.5% 52.0 / 1.5 19coding (903/912 nt)yvcJATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (10/7);  total (11/8)
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.52e-01

TTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTAC  >  CP000730/849640‑849907
                                                                                                                                           |                                                                                                                                  
ttAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGaaa                                                                                                                                 <  1:147638/141‑1 (MQ=255)
ttAGTAATTGACATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGaaa                                                                                                                                 <  2:353414/141‑1 (MQ=255)
     aaTTGCTATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGTAGTATTTGAATATAATGTTTATGTGAATCATAGGTATGCACATGTTGAAATTGGCGAGaaaaaaat                                                                                                                            >  2:92088/1‑139 (MQ=255)
             taggTTGTACGGGTGAACAACATAGATATGTAGCATTAGCAGAAAGACTAGGTAATTATCTAGATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAATAAAAAGAAAAAAATGAGACaaa                                                                                                                     <  2:1691/139‑1 (MQ=255)
                 ttGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATaaa                                                                                                                 >  1:15495/1‑141 (MQ=255)
                         gTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTAc                                                                                                         <  2:15495/141‑1 (MQ=255)
                         gTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTAc                                                                                                         >  1:178849/1‑141 (MQ=255)
                          tGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACt                                                                                                        <  2:217612/141‑1 (MQ=255)
                                     gATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGgtggtg                                                                                             >  1:75985/1‑141 (MQ=255)
                                                ttAGCAGAACGACGAGGTAATAAACTAAATGAAGTATTTGAATATAATGTTAATGTGCATCATAGGGACGCACATATTGAAAAAAAAAAAAAAAAAACGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTa                                                                                  <  2:307222/141‑1 (MQ=255)
                                                                     tATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTa                                                             >  2:166423/1‑141 (MQ=255)
                                                                         tAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAagag                                                         >  2:63312/1‑141 (MQ=255)
                                                                             tGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTTGCGAGAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGCTGGCACTGGCTTCTCAGCTATGCCTACGGGATTAAGAGACTTc                                                    >  2:225279/1‑141 (MQ=255)
                                                                                                     tgtgCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGt                             <  2:178849/141‑1 (MQ=255)
                                                                                                      gtgCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTa                            >  1:21631/1‑141 (MQ=255)
                                                                                                       tgactcatAGGCACGACCATATTGAAAGGAGCTAGAAAAAAATGAGACTAATATAATTTGTACTTATTGGTGGTGGCACTGGCTTATCAGTTATGGCTAAAAGAATAAGAGAATTCCCAAGTCATATTACGGCGATTGTaa                           <  2:271205/137‑1 (MQ=255)
                                                                                                        gCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAAc                          <  1:368420/141‑1 (MQ=255)
                                                                                                                gggACGCACATATTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTg                  >  2:207318/1‑141 (MQ=255)
                                                                                                                 ggACACACATATTGAAAGTGTCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGa                 <  2:354354/141‑1 (MQ=255)
                                                                                                                          tatTGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAAtggtgg        >  2:211374/1‑141 (MQ=255)
                                                                                                                            ttGAAAGTGGCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGa      >  1:50594/1‑141 (MQ=255)
                                                                                                                                aaGTGGCGAGAAAAAAATGAGACAAACAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTAc  >  1:356367/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                  
TTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTGGCGAGAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTAC  >  CP000730/849640‑849907

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