Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,764,980 G→A 9.6% V248V (GTC→GTT sufB ← component of SufBCD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,764,9800GA 9.6% 147.5 / 5.2 52V248V (GTC→GTTsufBcomponent of SufBCD complex
Reads supporting (aligned to +/- strand):  ref base G (28/19);  new base A (2/3);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ACTTCCACCACTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGCGTTAA  >  W3110S.gb/1764912‑1765044
                                                                    |                                                                
aCTTCCACCACTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGAc                                                                >  1:2842611/1‑71 (MQ=255)
 cTTCCACCACTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACg                                                               >  1:2953861/1‑71 (MQ=255)
    ccaccaCTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAg                                                            >  1:585775/1‑71 (MQ=255)
       ccacTGCCGCGTGTAACTGATAGCTGTGACGCACCGGAGCGGAACAGCCTTCAATGTAGCTg                                                                  <  1:907913/62‑1 (MQ=255)
       ccacTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTg                                                                  <  1:1251590/62‑1 (MQ=255)
       ccacTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTg                                                                  <  1:1172987/62‑1 (MQ=255)
        cacTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGt                                                        >  1:2509358/1‑71 (MQ=255)
         acTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGAc                                                                >  1:2144836/1‑62 (MQ=255)
          cTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCt                                                      <  1:1851546/71‑1 (MQ=255)
             ccGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTtcg                                                   >  1:1202434/1‑71 (MQ=255)
              cgcgTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTtcgt                                                  >  1:2439457/1‑71 (MQ=255)
                cgTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTtcgtcg                                                >  1:1582235/1‑71 (MQ=255)
                    tAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCtt                                                     >  1:340445/1‑62 (MQ=255)
                    tAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGcca                                            >  1:2228102/1‑71 (MQ=255)
                      aCTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTg                                                         >  1:1542626/1‑56 (MQ=255)
                      aCTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGccac                                           >  1:1502323/1‑70 (MQ=255)
                          aTAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAg                                        <  1:2057048/69‑1 (MQ=255)
                          aTAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAACTGTCTTCGTCGGCCACCAg                                        <  1:3011896/69‑1 (MQ=255)
                            aGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGa                                       >  1:3427909/1‑68 (MQ=255)
                              cTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGt                                  <  1:2299270/71‑1 (MQ=255)
                              cTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGt                                  >  1:455534/1‑71 (MQ=255)
                              cTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTAACGTAGCTGTCTTCGTCGGCCACCAGAATGGt                                  >  1:2603466/1‑71 (MQ=255)
                              cTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTAACGTAGCTGTCTTCGTCGGCCACCAGAATGGt                                  >  1:2833937/1‑71 (MQ=255)
                                 tCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTAACGTAGCTGTCTTCGTCGGCCACCAGaa                                      <  1:2047000/64‑1 (MQ=255)
                                 tCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTAACGTAGCTGTCTTCGTCGGCCACCAGaa                                      <  1:435843/64‑1 (MQ=255)
                                 tCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTAACGTAGCTGTCTTCGTCGGCCACCAGaa                                      <  1:1914306/64‑1 (MQ=255)
                                 tCACGCACCGGAGCCGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTtc                                                    <  1:2357204/50‑1 (MQ=255)
                                  cACGCACCGGAGCGGAACACCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAAt                                     <  1:916225/64‑1 (MQ=255)
                                   aCGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGccacca                                         >  1:508281/1‑59 (MQ=255)
                                   aCGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGccacca                                         >  1:2065934/1‑59 (MQ=255)
                                     gCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTc                            >  1:3516092/1‑70 (MQ=255)
                                     gCACCGGACCGGAACAGCCTTCAATGTAGCTGACTTAGCt                                                          >  1:3878983/1‑40 (MQ=38)
                                      cACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCaa                          <  1:3040487/71‑1 (MQ=255)
                                      cACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCaa                          <  1:3045648/71‑1 (MQ=255)
                                       aCCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAg                                                            >  1:3662243/1‑36 (MQ=255)
                                       aCCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCaaa                         <  1:3118890/71‑1 (MQ=255)
                                          ggAGCGGAACAGCCTTCAATGTAGGTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTg                      >  1:1389173/1‑71 (MQ=255)
                                          ggAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTtc                                                    <  1:2058112/41‑1 (MQ=255)
                                          ggAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTtc                                                    <  1:4365692/41‑1 (MQ=255)
                                          ggAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTg                      >  1:2920349/1‑71 (MQ=255)
                                          ggAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTg                      >  1:3410474/1‑71 (MQ=255)
                                             gCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGccc                   >  1:3501558/1‑71 (MQ=255)
                                             gCGGAACAGCCTACAATGTAGCTGACGTACCTGTCTtcgtcg                                                >  1:3930418/1‑42 (MQ=255)
                                               ggAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAAt                                     <  1:3837552/51‑1 (MQ=255)
                                                    aGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGttttt            <  1:2632909/71‑1 (MQ=255)
                                                      ccTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCa                           <  1:643636/54‑1 (MQ=255)
                                                      ccTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTc          >  1:3516022/1‑71 (MQ=255)
                                                          cAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGaaa                                     >  1:2548050/1‑39 (MQ=255)
                                                           aaTGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCaaa                         >  1:2958807/1‑51 (MQ=255)
                                                            aTGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGCGtt    <  1:4242146/71‑1 (MQ=255)
                                                            aTGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGCGtt    <  1:4402788/71‑1 (MQ=255)
                                                            aTGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGCGtt    <  1:2763415/71‑1 (MQ=255)
                                                             tGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGCGTTa   <  1:1385056/71‑1 (MQ=255)
                                                               tAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGCGTTaa  <  1:3740327/70‑1 (MQ=255)
                                                                   tGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGc       >  1:2583044/1‑61 (MQ=255)
                                                                    |                                                                
ACTTCCACCACTGCCGCGTGTAACTGATAGCTGTCACGCACCGGAGCGGAACAGCCTTCAATGTAGCTGACGTAGCTGTCTTCGTCGGCCACCAGAATGGTGCGCTCAAACTGCCCGGTTTTTTCTGCGTTAA  >  W3110S.gb/1764912‑1765044

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: