Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,683,336 T→A 9.2% L216Q (CTG→CAG)  folM → dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,683,3360TA 9.2% 176.9 / 8.1 65L216Q (CTG→CAG) folMdihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (24/35);  new base A (3/3);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCT  >  W3110S.gb/1683266‑1683405
                                                                      |                                                                     
gAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCt                                                                       >  1:1071571/1‑71 (MQ=255)
gAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCt                                                                       >  1:210307/1‑71 (MQ=255)
     tCGACAACAGCCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCg                                                                  <  1:2241962/71‑1 (MQ=255)
      cGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGa                                                                 >  1:342961/1‑71 (MQ=255)
       gACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGa                                                                 <  1:2568627/70‑1 (MQ=255)
           acaGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcact                                                            <  1:1770696/71‑1 (MQ=38)
             aGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactta                                                          <  1:1806978/71‑1 (MQ=255)
             aGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttactta                                                          >  1:5597229/1‑71 (MQ=255)
                cgcTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttact                                                        >  1:5227999/1‑70 (MQ=255)
                  cTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttacttac                                                     <  1:5126771/71‑1 (MQ=255)
                     aaTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTAcc                                                    >  1:4047369/1‑69 (MQ=255)
                     aaTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAg                                                  >  1:4793214/1‑71 (MQ=255)
                      aTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTAcc                                                    <  1:6808258/68‑1 (MQ=255)
                      aTAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGt                                                 >  1:734542/1‑71 (MQ=255)
                       tAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAt                                                                <  1:130760/55‑1 (MQ=255)
                           tCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttact                                                        <  1:7017093/59‑1 (MQ=255)
                           tCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttact                                                        <  1:4836137/59‑1 (MQ=255)
                           tCACTGATGAAAACCGCGCCTGGAGAGAGAGAAGTGATCGACCTGGTCGAttacttact                                                        <  1:2691018/59‑1 (MQ=255)
                               tgatgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttact                                                        <  1:5671488/55‑1 (MQ=255)
                                 atgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTg                                               >  1:2846874/1‑62 (MQ=255)
                                 atgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCAc                                      >  1:3012278/1‑71 (MQ=255)
                                  tgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttac                                                         <  1:7778372/51‑1 (MQ=255)
                                  tgaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttac                                                         <  1:7915812/51‑1 (MQ=255)
                                   gaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttactta                                                      <  1:3991789/53‑1 (MQ=255)
                                   gaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCg                                    >  1:5828251/1‑71 (MQ=255)
                                   gaAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCg                                    >  1:4518804/1‑71 (MQ=255)
                                    aaaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCt                                                                       >  1:2450936/1‑35 (MQ=255)
                                    aaaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTgg                                                                     >  1:5910059/1‑37 (MQ=255)
                                    aaaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttact                                                            <  1:4872397/46‑1 (MQ=255)
                                    aaaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCg                                    <  1:375413/70‑1 (MQ=255)
                                     aaaCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCt                                             >  1:6273206/1‑60 (MQ=255)
                                       aCCGCTCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACg                                <  1:4135084/71‑1 (MQ=255)
                                        ccGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTAcc                                                    <  1:5681772/50‑1 (MQ=255)
                                         cgcgCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCa                              <  1:5829692/71‑1 (MQ=255)
                                          gcgcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAg                             >  1:4163008/1‑71 (MQ=255)
                                           cgcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGa                                  >  1:1023030/1‑65 (MQ=255)
                                            gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttactta                                                      <  1:4415420/44‑1 (MQ=255)
                                            gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttactta                                                      <  1:6629190/44‑1 (MQ=255)
                                            gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGtt                           <  1:7699561/71‑1 (MQ=255)
                                            gcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGtt                           <  1:5771522/71‑1 (MQ=255)
                                               tGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAg                             <  1:2828774/66‑1 (MQ=255)
                                                 gCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCa                      <  1:7916406/71‑1 (MQ=255)
                                                 gCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCa                      <  1:3492613/71‑1 (MQ=255)
                                                  cGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTACCa                      >  1:1145420/1‑70 (MQ=255)
                                                  cGAGAAAGAAGTGATCGACATGGTTGATTACTTACTTACCAGTTGCTTTGTCAcc                                     <  1:700524/55‑1 (MQ=255)
                                                   gagaAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGc                                              <  1:7132107/45‑1 (MQ=255)
                                                     gaaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTc                                        >  1:5083430/1‑49 (MQ=255)
                                                     gaaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACtt                   >  1:7175647/1‑70 (MQ=255)
                                                     gaaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACtt                   >  1:717490/1‑70 (MQ=255)
                                                        agaagTGATCGACCTGGTCGATTACTTACTTACCAGt                                                 >  1:2702889/1‑37 (MQ=255)
                                                           agTAATCGACCAGGTCGATCACTTACTTACCAGTTGCt                                             <  1:7218221/38‑1 (MQ=14)
                                                             tGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgt          <  1:2725350/71‑1 (MQ=255)
                                                              gATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAg                             >  1:5792018/1‑51 (MQ=255)
                                                              gATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGTCGgtc         >  1:6229160/1‑71 (MQ=255)
                                                              aaTCGACCAGGTCGATCACTTACTGACCAGTTGCTTTGTCACCGGACGCa                              >  1:7475028/2‑50 (MQ=25)
                                                               aTCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATgg              <  1:7720455/65‑1 (MQ=255)
                                                               aTCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATgg              <  1:1145172/65‑1 (MQ=255)
                                                               aTCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATgg              <  1:4105053/65‑1 (MQ=255)
                                                                tCGACCAGGTCGATCACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACt                    >  1:5647434/1‑58 (MQ=39)
                                                                 cGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgt       >  1:4686814/1‑70 (MQ=255)
                                                                 cGACCAGGTCGATCACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgt       >  1:3209235/1‑70 (MQ=255)
                                                                  gACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgtc      <  1:7021083/70‑1 (MQ=255)
                                                                  gACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgtc      <  1:734081/70‑1 (MQ=255)
                                                                  gACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCa     >  1:6389292/1‑71 (MQ=255)
                                                                   aCCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCAt    <  1:6570227/71‑1 (MQ=255)
                                                                    ccTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATc   <  1:5096888/71‑1 (MQ=255)
                                                                     cTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATg               >  1:3443898/1‑58 (MQ=255)
                                                                     cTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCt  <  1:327235/71‑1 (MQ=255)
                                                                      tGGTCGATTACTTACTTTCCAGTTGCTTTGTCACCGGa                                  >  1:4436466/1‑38 (MQ=255)
                                                                      tGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGc                               >  1:1923306/1‑41 (MQ=255)
                                                                      tGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACt                    <  1:197967/52‑1 (MQ=255)
                                                                      tGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCg            <  1:7389270/60‑1 (MQ=255)
                                                                      |                                                                     
GAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCT  >  W3110S.gb/1683266‑1683405

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: