Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,237,778 C→T 86.7% R347R (CGG→CGA ycgB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,237,7780CT86.7% 109.1 / 9.6 45R347R (CGG→CGAycgBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (20/19);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCAGTCGGAACCGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCGCTGTACCACGGGCTGT  >  W3110S.gb/1237713‑1237845
                                                                 |                                                                   
agccagTCGGAACCGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGttt                                                                <  1:729183/71‑1 (MQ=255)
agccagTCGGAACCGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGttt                                                                <  1:2743248/71‑1 (MQ=255)
 gccagTCGGAACCGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTa                                                               <  1:1297359/71‑1 (MQ=255)
   cagTCGGAACCGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTAat                                                             <  1:259842/71‑1 (MQ=255)
    agTCGGAACCGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAata                                                            >  1:132984/1‑71 (MQ=255)
            ccGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGaa                                                    <  1:361967/71‑1 (MQ=255)
              ggCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACa                                                  >  1:496460/1‑71 (MQ=255)
              ggCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACa                                                  >  1:2778196/1‑71 (MQ=255)
                  atatCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACa                                                  >  1:1650937/1‑67 (MQ=255)
                     tCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGa                                            <  1:2500638/70‑1 (MQ=255)
                        gggAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAAccc                                        <  1:1598446/71‑1 (MQ=255)
                          gAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCg                                                                   <  1:41066/42‑1 (MQ=255)
                               aGTATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGaa                                                    <  1:2338575/52‑1 (MQ=255)
                                 tATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGttt                                                                <  1:2463432/38‑1 (MQ=255)
                                 tATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAAccc                                        <  1:2266428/62‑1 (MQ=255)
                                 tATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCa                                 <  1:76166/69‑1 (MQ=255)
                                 tATTTGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCa                                 <  1:1237410/69‑1 (MQ=255)
                                 tATTTGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTAATATCCTGGaa                                                    <  1:960926/50‑1 (MQ=255)
                                  aTTTGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATAc                              >  1:760523/1‑71 (MQ=255)
                                   tttGTCTTCTTCCGTTGGCGACTGACAAATTCGtt                                                                 <  1:2716881/35‑1 (MQ=255)
                                   tttGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACg                             >  1:1952668/1‑71 (MQ=255)
                                   tttGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAAAATGGCGAACCCGAGGGCATACg                             >  1:2039081/1‑71 (MQ=255)
                                   tttGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACg                             >  1:1464/1‑71 (MQ=255)
                                     tGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAat                                                             >  1:2138587/1‑37 (MQ=255)
                                     tGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAat                                                             >  1:2135281/1‑37 (MQ=255)
                                     tGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAat                                                             >  1:2797078/1‑37 (MQ=255)
                                     tGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATc                                                          >  1:1697999/1‑40 (MQ=255)
                                     tGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACggg                           >  1:1593695/1‑71 (MQ=255)
                                     tGTCTTCTTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACgg                            <  1:802403/70‑1 (MQ=255)
                                     tGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTAata                                                            >  1:19901/1‑38 (MQ=255)
                                       tcttctTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGtt                         <  1:2455141/71‑1 (MQ=255)
                                       tcttctTCCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGtt                         <  1:605227/71‑1 (MQ=255)
                                         ttcttcCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGaa                                           >  1:442485/1‑51 (MQ=255)
                                         ttcttcCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGa                       >  1:1043325/1‑71 (MQ=255)
                                         ttcttcCGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGa                       >  1:2523129/1‑71 (MQ=255)
                                         ttcttcAGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGa                       >  1:1877121/1‑71 (MQ=255)
                                              ccGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCg                  >  1:581492/1‑71 (MQ=255)
                                              ccGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCg                  >  1:2075576/1‑71 (MQ=255)
                                              ccGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCg                  >  1:2124461/1‑71 (MQ=255)
                                              ccGTTGGCGACTGACAAATTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCg                  >  1:2265979/1‑71 (MQ=255)
                                                 ttGGCGACTGACAAATTCGTTTAATATCCTGGAACATGTCGa                                            <  1:1986909/42‑1 (MQ=38)
                                                              aaTTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGtt                         <  1:2662847/48‑1 (MQ=255)
                                                              aaTTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCGCTGTACGACGGGCTGt  >  1:331051/1‑71 (MQ=255)
                                                              aaTTCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCGCTGTACCACGGGCTGt  <  1:1415949/71‑1 (MQ=255)
                                                              aaTTCGTTTAATATCATGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCGCTGTACCACGGGCTGt  <  1:2641641/71‑1 (MQ=255)
                                                                 |                                                                   
AGCCAGTCGGAACCGGCGATATCCGGGAACCAGTATTTGTCTTCTTCCGTTGGCGACTGACAAATCCGTTTAATATCCTGGAACATGGCGAACCCGAGGGCATACGGGTTGATGCCGCTGTACCACGGGCTGT  >  W3110S.gb/1237713‑1237845

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