Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,489,179 T→C 21.8% F77F (TTT→TTC ydcF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,489,1790TC21.8% 141.4 / 37.3 72F77F (TTT→TTCydcFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base C (8/7);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

CATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAC  >  W3110S.gb/1489109‑1489247
                                                                      |                                                                    
cATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACtttt                                                                      <  1:3192212/71‑1 (MQ=255)
 aTGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACttttt                                                                     >  1:4158188/1‑71 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:3403892/70‑1 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:3538726/70‑1 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:1815633/1‑70 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:1732174/1‑70 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:2957406/1‑70 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   >  1:3325749/1‑70 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:1112960/70‑1 (MQ=255)
    tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTAGACAACTTTCTTg                                                                   >  1:3421982/1‑70 (MQ=255)
      aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGtat                                                                <  1:3261484/71‑1 (MQ=255)
      aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGtat                                                                <  1:286262/71‑1 (MQ=255)
      aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGtat                                                                <  1:4338747/71‑1 (MQ=255)
      aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGta                                                                 <  1:2485692/70‑1 (MQ=255)
            cccGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAgcgc                                                           <  1:2539955/70‑1 (MQ=255)
                cgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACttttt                                                                     >  1:2725896/1‑56 (MQ=255)
                 gATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGGCATCGc                                                     <  1:807075/71‑1 (MQ=255)
                    cAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGca                                                    >  1:3089327/1‑69 (MQ=255)
                    aaGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGtat                                                                >  1:4296812/2‑57 (MQ=255)
                      gCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:3342240/52‑1 (MQ=255)
                       cAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAg                                                 <  1:924054/69‑1 (MQ=255)
                        aaaTTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCAt                                              >  1:3184283/1‑71 (MQ=255)
                         aaTTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATc                                             <  1:197511/71‑1 (MQ=255)
                          aTTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACtttt                                                                      <  1:2229858/45‑1 (MQ=255)
                           ttCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCg                                           <  1:624192/71‑1 (MQ=255)
                           ttCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCg                                           >  1:261250/1‑71 (MQ=255)
                            tCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGta                                                                 >  1:3275483/1‑48 (MQ=255)
                            tCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGta                                                                 >  1:292639/1‑48 (MQ=255)
                              cTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCAc                                        <  1:3957958/71‑1 (MQ=255)
                              cTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCAc                                        <  1:2868064/71‑1 (MQ=255)
                               tttACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTg                                                                   <  1:2857017/43‑1 (MQ=255)
                                  aCTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTac                                     <  1:105746/70‑1 (MQ=255)
                                  aCTGATTAGTGGTGGTATAGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTac                                     >  1:3592176/1‑70 (MQ=255)
                                   cTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTac                                     <  1:4270661/69‑1 (MQ=255)
                                      aTTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAacac                                <  1:1702818/71‑1 (MQ=255)
                                      aTTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAacac                                <  1:1349508/71‑1 (MQ=255)
                                      aTTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGtata                                                               >  1:1150502/1‑40 (MQ=255)
                                      aTTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTacaaca                                 <  1:4452820/70‑1 (MQ=255)
                                      aTTAGTGGTGGTATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAacac                                >  1:4113082/1‑71 (MQ=255)
                                         aGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTa                              <  1:2488775/70‑1 (MQ=255)
                                         aGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTa                              <  1:788635/70‑1 (MQ=255)
                                          gtggtgGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCAt                                                        >  1:1610947/1‑43 (MQ=255)
                                           tggtggTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAg                                                 <  1:4136814/49‑1 (MQ=255)
                                             gtggtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                          >  1:1201926/1‑70 (MQ=255)
                                             gtggtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGc                         <  1:137547/71‑1 (MQ=255)
                                             gtggtATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                          <  1:3328088/70‑1 (MQ=255)
                                             gtggtATCGGTCACTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                          <  1:3327998/70‑1 (MQ=255)
                                                gtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCAt                                              <  1:3270567/47‑1 (MQ=255)
                                                gtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGcacc                      <  1:542408/71‑1 (MQ=255)
                                                gtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGcac                       <  1:946873/70‑1 (MQ=255)
                                                  aTCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCAt                                              >  1:778857/1‑45 (MQ=255)
                                                  aTCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGcaccac                    >  1:1017508/1‑71 (MQ=255)
                                                   tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCg                                                      <  1:1847980/36‑1 (MQ=255)
                                                   tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAg                                                 <  1:2452753/41‑1 (MQ=255)
                                                   tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACt                                       <  1:3618263/51‑1 (MQ=255)
                                                   tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACt                   <  1:1261838/71‑1 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCATCGGACAGCAt                                              >  1:3181162/1‑42 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATc                                             >  1:2753080/1‑43 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCAc                                        >  1:1032952/1‑48 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTac                                     <  1:1833164/51‑1 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATc                            >  1:2107266/1‑60 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 >  1:2911474/1‑71 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 >  1:4024399/1‑71 (MQ=255)
                                                     ggTCACTCGACAACTTTTTTGTATAGCGCCAACGCACAGCATCCGCACTACAACACTATc                            >  1:1558871/1‑60 (MQ=255)
                                                      gTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCa                                               >  1:1336547/1‑40 (MQ=255)
                                                        cACTCGACAACTTTTTTGTATAGCGCCATCGCACAg                                                 <  1:586868/36‑1 (MQ=255)
                                                        cACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 >  1:3214061/1‑68 (MQ=255)
                                                        cACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 >  1:551018/1‑68 (MQ=255)
                                                         aCTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTgg                 <  1:510116/67‑1 (MQ=255)
                                                         aCTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgca               >  1:3408586/1‑69 (MQ=255)
                                                         aCTCGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcag              <  1:1635102/70‑1 (MQ=255)
                                                            cGACAACTTTCTTGTATAGCGCCATCGCACAGCATCCGCa                                         >  1:2690239/1‑40 (MQ=255)
                                                             gACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTacaaca                                 >  1:3846094/1‑47 (MQ=255)
                                                             gACAACTTTATTGTATAGCGCCATCGCAAAGCATCCGCACTACAACACTATCCGCACCACTGgcagag            >  1:1988552/1‑68 (MQ=255)
                                                              aCAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagcaga        <  1:1833837/71‑1 (MQ=255)
                                                                  cTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGc     >  1:2312659/1‑70 (MQ=255)
                                                                   tttcttGTATAGCGCCATCGCACAGCATCCGCACTac                                     <  1:37842/37‑1 (MQ=255)
                                                                    tttttGTATAGCGCCATCGCACAGCATCCGCACTacaac                                  >  1:214557/1‑39 (MQ=255)
                                                                    tttttGTATAGCGCCATCGCACAGCATCCGCACTacaac                                  >  1:2972638/1‑39 (MQ=255)
                                                                     ttttGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAc  >  1:183364/1‑70 (MQ=255)
                                                                     ttttGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAc  >  1:1224645/1‑70 (MQ=255)
                                                                      |                                                                    
CATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGAC  >  W3110S.gb/1489109‑1489247

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: