Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,317,696 G→A 7.5% H56Y (CAC→TAC)  yhbC ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,317,6960GA 7.5% 343.4 / 12.6 107H56Y (CAC→TAC) yhbCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (52/47);  new base A (4/4);  total (56/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 119 total aligned reads displayed.

GAGGAGACTTCCAGGTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGATGTG  >  W3110S.gb/3317629‑3317766
                                                                   |                                                                      
gaggagACTTCCAGGTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCt                                                                     >  1:3558383/1‑71 (MQ=255)
     gaCTTCCAGGTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACAt                                                                <  1:535255/71‑1 (MQ=255)
           cAGGTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGca                                                           >  1:862458/1‑70 (MQ=255)
            aGGTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGcaca                                                         <  1:1713454/71‑1 (MQ=255)
             ggTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGcaca                                                         >  1:2847723/1‑70 (MQ=255)
              gTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAat                                                       <  1:1376915/71‑1 (MQ=255)
               ttATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAat                                                       >  1:2243371/1‑70 (MQ=255)
               ttATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTAGCTCACATCAGCACAatc                                                      <  1:3109913/71‑1 (MQ=255)
               ttATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTAGCTCACATCAGCACAat                                                       >  1:4141915/1‑70 (MQ=255)
                tataAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTca                                                                   >  1:1435319/1‑57 (MQ=255)
                 ataAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAatcat                                                    >  1:1928424/1‑71 (MQ=255)
                 ataAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAatcat                                                    >  1:3338118/1‑71 (MQ=255)
                  taAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAatcatc                                                   <  1:2647133/71‑1 (MQ=255)
                    aGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:3327443/1‑71 (MQ=255)
                    aGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 <  1:3838311/71‑1 (MQ=255)
                    aGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 <  1:5009833/71‑1 (MQ=255)
                     gCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCt                                                                     >  1:3778362/1‑50 (MQ=255)
                     gCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAAtcatca                                                  <  1:1361030/69‑1 (MQ=255)
                     gCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:1670553/1‑70 (MQ=255)
                        aCGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACAtt                                             <  1:4509584/71‑1 (MQ=255)
                         cGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAatca                                                     <  1:2858406/62‑1 (MQ=255)
                         cGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTg                                            <  1:2588190/71‑1 (MQ=255)
                           gTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGa                                           >  1:4413662/1‑70 (MQ=255)
                             gATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTc                                                                    >  1:600907/1‑43 (MQ=255)
                              aTGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATGAACAtt                                             <  1:657667/65‑1 (MQ=255)
                              aTGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACAtt                                             <  1:3306854/65‑1 (MQ=255)
                              aTGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGc                                        <  1:248496/70‑1 (MQ=255)
                              aTGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGc                                        >  1:1597870/1‑70 (MQ=255)
                                 ggATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTca                                                                   >  1:3716535/1‑40 (MQ=255)
                                 ggATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTca                                                                   >  1:1182354/1‑40 (MQ=255)
                                 ggATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATc                                    >  1:1055890/1‑71 (MQ=255)
                                  gATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATc                                    <  1:3701208/70‑1 (MQ=255)
                                   aTCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCt                                   >  1:2332098/1‑70 (MQ=255)
                                     cTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGcc                                       <  1:4148837/64‑1 (MQ=255)
                                     cTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCtt                                  >  1:22396/1‑69 (MQ=255)
                                     cTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCtt                                  >  1:1197516/1‑69 (MQ=255)
                                     cTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCtt                                  >  1:1664171/1‑69 (MQ=255)
                                      ttCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGca                                                           <  1:2902654/43‑1 (MQ=255)
                                       tCAACATCCAGCACAGCACTTACCTGGTGGCTCACAt                                                                <  1:1672741/37‑1 (MQ=255)
                                       tCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:1601533/1‑52 (MQ=255)
                                       tCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:1333264/1‑52 (MQ=255)
                                       tCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:3748090/1‑52 (MQ=255)
                                       tCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCAc                               >  1:5075243/1‑70 (MQ=255)
                                        gaaCATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCtt                                  <  1:1450376/65‑1 (MQ=255)
                                        cAACATCCTGCACAGCACTTACCTGGTGGCTCACATCAGCACAatcatc                                                   <  1:277000/49‑1 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGa                                           >  1:3467853/1‑57 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTa                             >  1:2513687/1‑71 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTa                             >  1:4982867/1‑71 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTGGCTCAAATCAGCACAATCATCAAc                                                >  1:2184980/1‑52 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTAGCTCACATCAGCACAATCATCAACATTg                                            >  1:2755246/1‑56 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTAGCTCACATCAGCACAATCATCAACATTGATGCCATCtt                                  <  1:171147/66‑1 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTAGCTCACATCAGCACAATCATCAACATTGATGCCATCtt                                  <  1:2330443/66‑1 (MQ=255)
                                        cAACATCCAGCACAGCACTTACCTGGTAGCTCACATCAGCACAATCATCAACATTGATGCCATCtt                                  <  1:4339282/66‑1 (MQ=255)
                                        cAACATCCAGCAAAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTa                             >  1:4910804/1‑71 (MQ=255)
                                         aaCATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:1871360/1‑50 (MQ=255)
                                         aaCATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:2049260/1‑50 (MQ=255)
                                          aCATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATc                           <  1:4035537/71‑1 (MQ=255)
                                           cATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACAtt                                             <  1:3963517/52‑1 (MQ=255)
                                           cATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATc                           >  1:90459/1‑70 (MQ=255)
                                              ccAGCACAGCACTTACCTGGTGGCTCACATCAGCACAatca                                                     >  1:3431835/1‑41 (MQ=255)
                                               cagcacagcaCTTACCTGGTGGCTCACATCAGCACAatcatc                                                   <  1:1576551/42‑1 (MQ=255)
                                               cagcacagcaCTTACCTGGTGGCTCACATCAGCACAatcatc                                                   <  1:1246267/42‑1 (MQ=255)
                                               cagcacagcaCTTACCTGGTGGCTCACATCAGCACAATCATCAACa                                               >  1:639774/1‑46 (MQ=255)
                                               cagcacagcaCTTACCTGGTGGCTCACATCAGCACAATCATCAACa                                               >  1:2493006/1‑46 (MQ=255)
                                               cagcacagcaCTTACCTGGTGGCTCACATCAGCACAATCATCAACa                                               >  1:3512824/1‑46 (MQ=255)
                                                agcacagcacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAAtata                     <  1:3741877/71‑1 (MQ=255)
                                                agcacagcacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAAtata                     <  1:3216429/71‑1 (MQ=255)
                                                 gcacagcacTTACGTGGTGGCTCACATCAGCACAatcatc                                                   <  1:1171578/40‑1 (MQ=255)
                                                 gcacagcacTTACCTGGTGGCTCACATCAGCACAatcatc                                                   <  1:4279523/40‑1 (MQ=255)
                                                 gcacagcacTTACCTGGTGGCTCACATCAGCACAatcatc                                                   <  1:4095607/40‑1 (MQ=255)
                                                 gcacagcacTTACCTGGTGGCTCACATCAGCACAATCATCaa                                                 >  1:785865/1‑42 (MQ=255)
                                                 gcacagcacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCaa                         <  1:1835794/66‑1 (MQ=255)
                                                 gcacagcacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCaa                         <  1:4077399/66‑1 (MQ=255)
                                                  cacagcacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAata                       >  1:3093812/1‑67 (MQ=255)
                                                       cacTTATCTGGTGGCTCTCATCAGCACAATCATCAACATTGATGCCATCTTCACTATCaa                         <  1:4776587/60‑1 (MQ=255)
                                                       cacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATgcgc              <  1:2277724/71‑1 (MQ=255)
                                                       cacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATgcg               >  1:1648054/1‑70 (MQ=255)
                                                         cTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTa                             <  1:2092093/54‑1 (MQ=255)
                                                         cTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTa                             <  1:2086085/54‑1 (MQ=255)
                                                         cTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTa                             <  1:3335509/54‑1 (MQ=255)
                                                          ttACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAAtata                     <  1:686784/61‑1 (MQ=255)
                                                          ttACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAAtata                     <  1:3178831/61‑1 (MQ=255)
                                                           tACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCAc                               >  1:5030965/1‑50 (MQ=255)
                                                           tACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCAc                               >  1:211760/1‑50 (MQ=255)
                                                            aCCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCAt                                     >  1:2475522/1‑43 (MQ=255)
                                                            aCCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAgtgt         >  1:1926294/1‑71 (MQ=255)
                                                            aCCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAgtgt         >  1:4685079/1‑71 (MQ=255)
                                                            aCCTGGTAGCTCACATCAGCACAATCATCAACATTGa                                           >  1:1303917/1‑37 (MQ=255)
                                                             ccTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGtgtg        <  1:3735229/71‑1 (MQ=255)
                                                             ccTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGtgtg        <  1:3684459/71‑1 (MQ=255)
                                                             ccTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGtgtg        <  1:53635/71‑1 (MQ=255)
                                                                ggtggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCa                                >  1:1739579/1‑44 (MQ=255)
                                                                ggtggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCAc                               >  1:1076731/1‑45 (MQ=255)
                                                                ggtggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAata                       >  1:3113511/1‑53 (MQ=255)
                                                                ggtggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATgcgc              >  1:1341420/1‑62 (MQ=255)
                                                                 gtggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTAt                            <  1:643663/47‑1 (MQ=255)
                                                                 gtggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGAt     <  1:4546681/70‑1 (MQ=255)
                                                                   ggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCa                          >  1:472384/1‑47 (MQ=255)
                                                                   ggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGAtgtg  >  1:4311335/1‑71 (MQ=255)
                                                                   ggCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGAtgtg  >  1:3419738/1‑71 (MQ=255)
                                                                   |                                                                      
GAGGAGACTTCCAGGTTATAAGCAACGGTGATGGGATCTTCAACATCCAGCACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGATGTG  >  W3110S.gb/3317629‑3317766

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: