Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 436,250 C→T 12.2% I146I (ATC→ATT pgpA → phosphatidylglycerophosphatase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb436,2500CT12.2% 126.8 / 7.7 49I146I (ATC→ATTpgpAphosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base T (4/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATATGTGGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCAC  >  W3110S.gb/436183‑436307
                                                                   |                                                         
atatGTGGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATcgac                                                        <  1:1552887/71‑1 (MQ=255)
   tgtgGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATcgacga                                                      <  1:2441151/70‑1 (MQ=255)
   tgtgGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATcgacga                                                      <  1:1260081/70‑1 (MQ=255)
   tgtgGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGat                                                     <  1:3323104/71‑1 (MQ=255)
    gtgGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATcgacga                                                      <  1:3034693/69‑1 (MQ=255)
     tgGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGatat                                                   >  1:2359885/1‑71 (MQ=255)
     tgGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGatat                                                   >  1:173031/1‑71 (MQ=255)
       gAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATAtt                                                  <  1:3156763/70‑1 (MQ=255)
       gAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATAtt                                                  <  1:2151609/70‑1 (MQ=255)
           ccGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGc                                             >  1:1616094/1‑71 (MQ=255)
                 ccgctccgctGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGt                                       >  1:1576412/1‑71 (MQ=255)
                 ccGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGt                                       >  1:293759/1‑71 (MQ=255)
                 ccGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGt                                       >  1:690136/1‑71 (MQ=255)
                   gATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATcgacga                                                      <  1:1874970/54‑1 (MQ=255)
                   gATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATcgacga                                                      <  1:2645513/54‑1 (MQ=255)
                   gATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATcgacga                                                      <  1:3783038/54‑1 (MQ=255)
                   gATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGa                                     >  1:3517142/1‑71 (MQ=255)
                    aTCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCgg                                          <  1:1252554/65‑1 (MQ=255)
                    aTCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGAt                                    <  1:2074440/71‑1 (MQ=255)
                      cctcTGGTTTGATCGCAATGTACATGGCGGCATGGGGATCATGATTGACGATATTGTCGCCGGGGTGAt                                    <  1:3488890/66‑1 (MQ=255)
                      ccGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCggg                                         >  1:1943827/1‑64 (MQ=255)
                      ccGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCgg                                          >  1:2548833/1‑63 (MQ=255)
                      ccGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGAACATGATCGACGATATTGTCGCCGGGGTGAt                                    >  1:2877369/1‑69 (MQ=255)
                      ccGCTGGTTTGATCGCAATGTACATGGCGGCATGGGGATCATGATTGACGATATTGTCGCCGGGGTGAt                                    <  1:3120539/69‑1 (MQ=255)
                         cTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCggg                                         <  1:1244088/61‑1 (MQ=255)
                         cTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCg                               <  1:3676262/71‑1 (MQ=255)
                          tGGTTTGATCGCAATGTACATGGCGGCATGGGGATCATGATTGACGATATTGTCg                                              >  1:1456380/1‑55 (MQ=255)
                          tGGTTTGATCGCAATGTACATGGCGGCATGGGGATCATGATTGACGATATTGTCGc                                             >  1:3652166/1‑56 (MQ=255)
                          tGGTTTGATCGCAATGTACATGGCGGCATGGGGATCATGATTGACGATATTGTCGc                                             >  1:992397/1‑56 (MQ=255)
                          tGGTTTGATCGCAATGTACATGGCGGCATGGGGATCATGATTAACGATATTGTCGc                                             >  1:960356/1‑56 (MQ=255)
                               tGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCa                         >  1:3384376/1‑71 (MQ=255)
                                 aTCGCAATGTGCATGGCGGCATGGGGATCCTGATCGACGATAtt                                                  <  1:3459997/44‑1 (MQ=255)
                                 aTCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATc                       <  1:1968557/71‑1 (MQ=255)
                                       aTGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCg                               <  1:2773662/57‑1 (MQ=255)
                                       aTGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCg                               <  1:2682461/57‑1 (MQ=255)
                                          tgCATGGCGGCATGGGGATCATGATCGACGATAtt                                                  >  1:3261602/1‑35 (MQ=255)
                                          tgCATGGCGGCATGGGGATCATGATCGACGATAtt                                                  >  1:1558567/1‑35 (MQ=255)
                                            cATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTa             <  1:3433998/70‑1 (MQ=255)
                                            cATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTa             <  1:651459/70‑1 (MQ=255)
                                              tGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGt                                       <  1:3254182/42‑1 (MQ=255)
                                              tGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGt                                       <  1:345392/42‑1 (MQ=255)
                                              tGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATc                       <  1:1763642/58‑1 (MQ=255)
                                               ggcggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCATGTATTTTATCgg         >  1:1899373/1‑71 (MQ=255)
                                                gcggcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGt        >  1:2089676/1‑71 (MQ=255)
                                                   gcATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGtcat     <  1:1740284/71‑1 (MQ=255)
                                                     aTGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGtcatca   <  1:967468/71‑1 (MQ=255)
                                                        gggATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCg          >  1:821572/1‑61 (MQ=255)
                                                        gggATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCAc  >  1:3853186/1‑69 (MQ=255)
                                                         ggATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGtc       >  1:2774530/1‑63 (MQ=255)
                                                                gATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTAt                 >  1:1011874/1‑46 (MQ=255)
                                                                gATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTAt                 >  1:3371960/1‑46 (MQ=255)
                                                                gATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCAGTa                  >  1:3569840/1‑45 (MQ=255)
                                                                 aTCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTAt            <  1:258947/50‑1 (MQ=255)
                                                                 aTCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTAt            <  1:3375732/50‑1 (MQ=255)
                                                                 aTCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTAt            <  1:2301899/50‑1 (MQ=255)
                                                                 aTCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCgg         <  1:2368652/53‑1 (MQ=255)
                                                                   |                                                         
ATATGTGGAAGCCGTGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCAC  >  W3110S.gb/436183‑436307

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: