Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,434,349 T→A 40.9% intergenic (+112/‑197) fklB → / → cycA FKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,434,3490TA40.9% 41.4 / 80.6 66intergenic (+112/‑197)fklB/cycAFKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base A (16/11);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCT  >  W3110S.gb/4434281‑4434417
                                                                      |                                                                    
cAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTAtttt                                                                     <  1:637142/70‑1 (MQ=255)
    aaTAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTtgt                                                                  >  1:146923/1‑69 (MQ=255)
     aTAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTtgtg                                                                 >  1:1863760/1‑69 (MQ=255)
     aTAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTtgtg                                                                 >  1:1784681/1‑69 (MQ=255)
      tAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTtgtg                                                                 <  1:4192346/68‑1 (MQ=255)
      tAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGc                                                              <  1:2464377/71‑1 (MQ=255)
          cccTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGtt                                                          >  1:4535418/1‑71 (MQ=255)
               cAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTcgcg                                                     >  1:1793936/1‑71 (MQ=255)
               cAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTcg                                                       <  1:2588855/69‑1 (MQ=255)
                aGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTcgcg                                                     <  1:4331295/70‑1 (MQ=255)
                 gTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGtt                                                   <  1:4178063/71‑1 (MQ=255)
                  tttGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGtt                                                   <  1:3209798/71‑1 (MQ=255)
                    tGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTtg                                                                   >  1:1573160/1‑52 (MQ=255)
                     gCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTAt                                                 <  1:2284368/70‑1 (MQ=255)
                     gCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATc                                                <  1:3513384/70‑1 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGt                                           <  1:3213866/71‑1 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCAcc                                             >  1:4230630/1‑69 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCg                                            >  1:2803287/1‑70 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGt                                           <  1:1516472/71‑1 (MQ=255)
                          aCGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTg                                          <  1:4412172/71‑1 (MQ=255)
                          aCGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCAc                                              <  1:606181/67‑1 (MQ=255)
                            ggCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTAttttt                                                                    >  1:4424758/1‑43 (MQ=255)
                              cgcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGatat                                      <  1:3330204/71‑1 (MQ=255)
                               gcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGatat                                      <  1:2853945/70‑1 (MQ=255)
                                cTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                    <  1:619261/71‑1 (MQ=255)
                                  aTTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATtg                                                                   >  1:1894304/1‑38 (MQ=255)
                                  aTTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTcg                                                       >  1:1352934/1‑50 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTg                                                            <  1:4795263/44‑1 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGacac                                 >  1:3635062/1‑71 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGtt                                                   <  1:4109109/53‑1 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGt                                                    >  1:3452131/1‑52 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGt                                                    >  1:3097042/1‑52 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                                 >  1:548497/1‑55 (MQ=255)
                                   ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGacac                                 >  1:251960/1‑71 (MQ=255)
                                    tttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                >  1:4330823/1‑71 (MQ=255)
                                    tttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                            >  1:3219152/1‑43 (MQ=255)
                                     ttGTGTTGCAAAGACTCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                >  1:4357217/1‑71 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                    >  1:2910569/1‑66 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                    >  1:2904758/1‑66 (MQ=255)
                                     ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                >  1:3967624/1‑70 (MQ=255)
                                           tGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACttt                          >  1:823897/1‑70 (MQ=255)
                                            gCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTa                         <  1:3648060/70‑1 (MQ=255)
                                            gCAAAGACCCAGTAAGAGTGTATTTTTGTGAGCTTTTTCGCGTTATCACCGTGatat                                      <  1:595538/57‑1 (MQ=255)
                                             cAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTaaa                       >  1:3982237/1‑71 (MQ=255)
                                                    ccccGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                                 <  1:707617/39‑1 (MQ=255)
                                                    ccccGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                                 <  1:2155268/39‑1 (MQ=255)
                                                      cccGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaaaa                >  1:982532/1‑71 (MQ=255)
                                                       ccGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGtt                                                   <  1:35904/35‑1 (MQ=255)
                                                       ccGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGtt                                                   <  1:4405150/35‑1 (MQ=255)
                                                        cGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAt               >  1:1806911/1‑70 (MQ=255)
                                                        cGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTa                                                  <  1:3948552/35‑1 (MQ=255)
                                                        cGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcactc                               >  1:455203/1‑54 (MQ=255)
                                                         gTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCa                                               >  1:283677/1‑37 (MQ=255)
                                                         gTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCAc                                              >  1:3259120/1‑38 (MQ=255)
                                                         gTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAtt              >  1:4280798/1‑70 (MQ=255)
                                                          tAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACa                     <  1:1826626/62‑1 (MQ=255)
                                                          tAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTaa            >  1:1680523/1‑71 (MQ=255)
                                                          tAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGt                                           <  1:696512/40‑1 (MQ=255)
                                                            aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAAc           <  1:4142510/70‑1 (MQ=255)
                                                             gCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACAt         <  1:644850/71‑1 (MQ=255)
                                                               gtgtATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                                         >  1:1836701/1‑37 (MQ=255)
                                                               gtgtATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                                         >  1:693404/1‑37 (MQ=255)
                                                               gtgtATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                                         >  1:4475818/1‑37 (MQ=255)
                                                               gtgtATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                    >  1:2343300/1‑42 (MQ=255)
                                                                tgtATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAg      <  1:4508686/71‑1 (MQ=255)
                                                                 gtATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACAt         >  1:1041129/1‑67 (MQ=255)
                                                                  tATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaaa                 <  1:1048289/58‑1 (MQ=255)
                                                                   aTTTTTGTGAGCTGTTTCGCGTTATCAACGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATc   >  1:1189746/1‑71 (MQ=255)
                                                                    tttttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaa                  <  1:3379463/55‑1 (MQ=255)
                                                                     ttttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGa     >  1:1521308/1‑67 (MQ=255)
                                                                     ttttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGa     >  1:4615048/1‑67 (MQ=255)
                                                                     ttttGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCt  >  1:3263172/1‑70 (MQ=255)
                                                                      |                                                                    
CAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCT  >  W3110S.gb/4434281‑4434417

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: