Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,806,448 T→C 58.8% pseudogene (1641/1900 nt)
pseudogene (1157/1416 nt)
arpB →
arpB →
ECK1718:JW5278+JW1710:b4494; hypothetical protein
ECK1718:JW1710:b1721; hypothetical protein, C‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,806,4480TC58.8% 17.2 / 37.8 34pseudogene (1641/1900 nt)
pseudogene (1157/1416 nt)
arpB
arpB
ECK1718:JW5278+JW1710:b4494; hypothetical protein
ECK1718:JW1710:b1721; hypothetical protein, C‑ter fragment
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base C (6/14);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TACAAGTTGAGCAAAGCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCTTTT  >  W3110S.gb/1806381‑1806515
                                                                   |                                                                   
tACAAGTTGAGCAAAGCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATcata                                                                  >  1:982618/1‑71 (MQ=255)
tACAAGTTGAGCAAAGCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATcata                                                                  <  1:438283/71‑1 (MQ=255)
        gAGCAAAGCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGcc                                                           <  1:1341071/70‑1 (MQ=255)
        gAGCAAAGCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGAATCATACATCGcc                                                           <  1:1937583/70‑1 (MQ=255)
              aGCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACACCTGCATCATACATCGCCATGAGc                                                     >  1:2105401/1‑70 (MQ=255)
               gCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATc                                                              >  1:1749157/1‑60 (MQ=255)
               gCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGc                                                     >  1:638418/1‑69 (MQ=255)
                cTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAg                                                  <  1:300006/71‑1 (MQ=255)
                cTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAg                                                  >  1:2621545/1‑71 (MQ=255)
                  aaCATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGcc                                                           >  1:1398623/1‑60 (MQ=255)
                    catcatTGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGt                                               <  1:2428922/70‑1 (MQ=255)
                      tcatGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAAt                                            >  1:936078/1‑71 (MQ=255)
                      tcatGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAAt                                            >  1:2328860/1‑71 (MQ=255)
                       catGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGc                                                            <  1:390312/54‑1 (MQ=255)
                        atGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATg                                           <  1:653470/70‑1 (MQ=255)
                         tGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCa                                                          >  1:741851/1‑54 (MQ=255)
                         tGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATg                                           <  1:2485370/69‑1 (MQ=255)
                         tGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGaa                                         <  1:1697482/71‑1 (MQ=255)
                         tGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGaa                                         <  1:192525/71‑1 (MQ=255)
                         tGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGaa                                         <  1:1966074/71‑1 (MQ=255)
                         tGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGaa                                         <  1:242700/71‑1 (MQ=255)
                         tGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGaa                                         <  1:1657302/71‑1 (MQ=255)
                          ggATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGaa                                         <  1:890878/70‑1 (MQ=255)
                            aTTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGAAGAc                                      <  1:192275/71‑1 (MQ=255)
                                  taAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGAAGAcgtcgt                                 <  1:1259732/70‑1 (MQ=255)
                                   aaaaGGAGCTACAGCACAGGGAACCCCTGCATcata                                                                  >  1:298305/1‑36 (MQ=255)
                                     aaGGAGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTa                             <  1:2486649/71‑1 (MQ=255)
                                         aGCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCt                          >  1:1369768/1‑70 (MQ=255)
                                          gCTACAGCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTAt                            <  1:128887/67‑1 (MQ=255)
                                                gCACAGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATc                  <  1:1086437/71‑1 (MQ=255)
                                                    aGGGAACCCCTGCATCATACATCGCCATGAGCAAGGGt                                               <  1:759712/38‑1 (MQ=255)
                                                     gggAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTtata                      >  1:894405/1‑62 (MQ=255)
                                                     gggAACCCCTGCATCATACATCGCCATGAGCAAGGGt                                               <  1:1100276/37‑1 (MQ=255)
                                                       gAACCCCTGCATCATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGa                <  1:2749112/66‑1 (MQ=255)
                                                            ccTGCATCATACATCGCCATGAGCAAGGGTAATGAAGa                                       >  1:2253288/1‑38 (MQ=255)
                                                               gCATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCttt   >  1:358005/1‑71 (MQ=255)
                                                               gCATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCttt   >  1:965327/1‑71 (MQ=255)
                                                                catcatACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCtttt  >  1:1701203/1‑71 (MQ=255)
                                                                   |                                                                   
TACAAGTTGAGCAAAGCTAACATCATGGATTTATTAAAAGGAGCTACAGCACAGGGAACCCCTGCATTATACATCGCCATGAGCAAGGGTAATGAAGACGTCGTGTTATCTTATATATCGACGCTGGGTGCTTTT  >  W3110S.gb/1806381‑1806515

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