Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,472,611 G→A 35.3% G144G (GGC→GGT btuB ← vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,472,6110GA35.3% 35.2 / 30.1 34G144G (GGC→GGTbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base A (8/4);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCCACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGAC  >  W3110S.gb/3472550‑3472673
                                                             |                                                              
cccACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCa                                                       >  1:20440/1‑71 (MQ=255)
cccACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCACCTATTGCa                                                       >  1:4526959/1‑71 (MQ=255)
      ctgctgAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGaa                                                 >  1:2895332/1‑71 (MQ=255)
           gAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCAt                                             >  1:1056078/1‑70 (MQ=255)
           gAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCAt                                             >  1:4560620/1‑70 (MQ=255)
            aaaTTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCAt                                                      >  1:1469952/1‑60 (MQ=255)
            aaaTTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATa                                            <  1:5009425/70‑1 (MQ=255)
            aaaTTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATa                                            <  1:1767206/70‑1 (MQ=255)
            aaaTTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATa                                            <  1:1767627/70‑1 (MQ=255)
               tttCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACa                                        <  1:4708347/71‑1 (MQ=255)
                   cgttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCg                                     <  1:505461/70‑1 (MQ=255)
                    gttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCg                                     >  1:350622/1‑69 (MQ=255)
                    gttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCg                                     >  1:1593754/1‑69 (MQ=255)
                     ttccgGGTTCATCGCGCGGCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAg                                  <  1:4236116/71‑1 (MQ=255)
                      tccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAgc                                 >  1:4607978/1‑71 (MQ=255)
                        cgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAgcgc                               <  1:2591507/71‑1 (MQ=255)
                        cgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAgcgc                               <  1:1889115/71‑1 (MQ=255)
                                tCGCGCGTCGTGATTATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACg                       <  1:1142009/71‑1 (MQ=255)
                                tCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATa                                            <  1:1251589/50‑1 (MQ=255)
                                  gcgcgTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGa                                   >  1:3262583/1‑57 (MQ=255)
                                          tgatgaTATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTc              <  1:3155712/70‑1 (MQ=255)
                                          tgatgaTATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCa             <  1:2492819/71‑1 (MQ=255)
                                               atatTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCa                         >  1:4960626/1‑54 (MQ=255)
                                               atatTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACg        >  1:3390439/1‑71 (MQ=255)
                                               atatTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACg        >  1:3197616/1‑71 (MQ=255)
                                               atatTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACg        >  1:163783/1‑71 (MQ=255)
                                               atatTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACg        >  1:798053/1‑71 (MQ=255)
                                                tatTCACCACCCCGCCTATTGCATCGGAACCATAAACAg                                       <  1:5081492/39‑1 (MQ=255)
                                                    caccacCCCGCCTATTGCATCGGAACCATAAACAGCGGat                                  >  1:2809813/1‑39 (MQ=255)
                                                    caccacCCCGCCTATTGCATCGGAACCATAAACAGCGGAg                                  >  1:3097151/1‑40 (MQ=255)
                                                    caccacCCCGCCTATTGCATCGGAACCATAAACAGCGGAGcgcg                              >  1:5047924/1‑44 (MQ=255)
                                                    caccacCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGa   >  1:3537930/1‑71 (MQ=255)
                                                     accaccCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGa   <  1:70522/70‑1 (MQ=255)
                                                     accaccCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGAc  <  1:4785761/71‑1 (MQ=255)
                                                             |                                                              
CCCACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGAC  >  W3110S.gb/3472550‑3472673

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: