Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,796,459 G→A 11.6% A136A (GCC→GCT btuC ← vitamin B12 transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,796,4590GA11.6% 162.1 / 10.5 60A136A (GCC→GCTbtuCvitamin B12 transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (2/5);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

TTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTG  >  W3110S.gb/1796389‑1796528
                                                                      |                                                                     
ttAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACgg                                                                       >  1:1070313/1‑71 (MQ=255)
ttAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACgg                                                                       >  1:1065070/1‑71 (MQ=255)
 tAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACggg                                                                      >  1:3520762/1‑71 (MQ=255)
    tGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGa                                                                   <  1:1654795/71‑1 (MQ=255)
      cgcTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGaaa                                                                 >  1:36785/1‑71 (MQ=255)
       gcTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGaa                                                                  <  1:4254852/69‑1 (MQ=255)
          aCAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGt                                                              >  1:444850/1‑70 (MQ=255)
          aCAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGt                                                              >  1:2698123/1‑70 (MQ=255)
           cAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTaa                                                            <  1:1616057/71‑1 (MQ=255)
           cAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTa                                                             <  1:4388818/70‑1 (MQ=255)
            aGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGt                                                              >  1:2076198/1‑68 (MQ=255)
              aTAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAagag                                                         <  1:2247007/71‑1 (MQ=255)
                aaTCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTa                                                       >  1:3020388/1‑71 (MQ=255)
                  tCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTAtt                                                     <  1:1251284/71‑1 (MQ=255)
                    ccTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTa                                                    <  1:1768771/70‑1 (MQ=255)
                    ccTAATCCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTa                                                    <  1:2855740/70‑1 (MQ=255)
                           cAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAagag                                                         >  1:3302723/1‑58 (MQ=255)
                           cAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAagag                                                         >  1:552981/1‑58 (MQ=255)
                           cAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTaaa                                                  >  1:3987734/1‑65 (MQ=255)
                            aaCGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGtgatga                                           >  1:3274916/1‑71 (MQ=255)
                             aCGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTgatgat                                          <  1:799057/71‑1 (MQ=255)
                                    ccagcAATAACCGACTGGTCGAAAGATGACGACGAGCGAAACGTAAGAGTATTAAAg                                                 >  1:1602998/1‑57 (MQ=255)
                                    ccagcAATAACCGACTGGTCGAAAGATGACGACGAGCGAAACGTAAGAGTATTAAAg                                                 >  1:3041411/1‑57 (MQ=255)
                                       gcAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAgcgcgc                                  >  1:1575404/1‑69 (MQ=255)
                                         aaTAACCGACTGGTCGAAAGATGACGACGGGCGAAACTTAAGAGTATTAAAGTGATGATAAGCGCGCCAg                               >  1:1308495/1‑70 (MQ=255)
                                         aaTAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCa                                >  1:1596463/1‑69 (MQ=255)
                                         aaTAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAg                               >  1:2200280/1‑70 (MQ=255)
                                           tAACCGACTGGTCGAAAGATGACGACGGGCGAAACGt                                                              <  1:3985051/37‑1 (MQ=255)
                                           tAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCa                             >  1:1568246/1‑70 (MQ=255)
                                             aCCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTaaa                                                  >  1:3958269/1‑47 (MQ=255)
                                              ccGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCg                         <  1:2729797/71‑1 (MQ=255)
                                              ccGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCg                         <  1:3763479/71‑1 (MQ=255)
                                              ccGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCg                         <  1:1911316/71‑1 (MQ=255)
                                               cGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTAt                                                      >  1:3734936/1‑41 (MQ=255)
                                               cGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCg                         <  1:1307177/70‑1 (MQ=255)
                                                gACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGca                       <  1:563919/71‑1 (MQ=255)
                                                gACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGca                       <  1:3551057/71‑1 (MQ=255)
                                                gACTGGTCGAAAGATGAAGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAg                                       <  1:2120181/55‑1 (MQ=255)
                                                 aCTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTa                                                    <  1:52164/41‑1 (MQ=255)
                                                 aCTGGTCGAAAGATGACGACGAGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCaa                            <  1:3551292/65‑1 (MQ=255)
                                                 aCTGGTCGAAAGATGACGACGAGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCaa                            <  1:757780/65‑1 (MQ=255)
                                                 aCTGGTCGAAAGATGACGACGAGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCaa                            <  1:2379426/65‑1 (MQ=255)
                                                 aCTGGTCGAAAGATGACGACGAGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCaa                            <  1:3909338/65‑1 (MQ=255)
                                                     gTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGtga                                              >  1:566239/1‑43 (MQ=255)
                                                     gTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGcaca                     >  1:4431569/1‑68 (MQ=255)
                                                      tCGAAAGATGACGACGGGCGAAACGTAAGAGTATTa                                                    <  1:3304851/36‑1 (MQ=255)
                                                      tCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAAt                           <  1:3602127/61‑1 (MQ=255)
                                                       cGAAAGATGACGAAGGGCGAAACGTAAGAGTATTAAAGtga                                              >  1:2834547/1‑41 (MQ=255)
                                                          aaGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATTAg                                       <  1:1164456/45‑1 (MQ=255)
                                                          aaGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAgc                                      <  1:2374500/46‑1 (MQ=255)
                                                          aaGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTa             <  1:617636/71‑1 (MQ=255)
                                                           aGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAg            >  1:2236385/1‑71 (MQ=255)
                                                           aGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAg            >  1:2009422/1‑71 (MQ=255)
                                                            gATGACGACGAGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCaa                            <  1:3154495/54‑1 (MQ=255)
                                                             aTGACGACGGGCGAAACGTAAGAGTATTAAAGTgatgat                                          <  1:4434879/39‑1 (MQ=255)
                                                              tGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAg                  <  1:832778/62‑1 (MQ=255)
                                                              tGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgcg          >  1:1558343/1‑70 (MQ=255)
                                                               gacgacGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGc                 <  1:1877813/62‑1 (MQ=255)
                                                                acgacgGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAgc           >  1:1235628/1‑67 (MQ=255)
                                                                acgacgGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGcc        <  1:1156567/70‑1 (MQ=255)
                                                                   acgGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGt    <  1:394750/71‑1 (MQ=255)
                                                                    cgGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGcaca                     >  1:1721196/1‑53 (MQ=255)
                                                                    cgGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGtt   <  1:2015252/71‑1 (MQ=255)
                                                                     gggCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTg  <  1:3406631/71‑1 (MQ=255)
                                                                      |                                                                     
TTAGTGCGCTACAGATAATCCCTAATGCAACGCCAGCCAGCAATAACCGACTGGTCGAAAGATGACGACGGGCGAAACGTAAGAGTATTAAAGTGATGATAAGCGCGCCAGCAATCGCACACAGCCCTAGCGCCCAGTTG  >  W3110S.gb/1796389‑1796528

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: