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

CAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCT  >  W3110S.gb/4434281‑4434417
                                                                      |                                                                    
cAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTAtttt                                                                     <  1:637142/70‑1 (MQ=255)
    aaTAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTtgt                                                                  >  1:146923/1‑69 (MQ=255)
     aTAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTtgtg                                                                 >  1:1863760/1‑69 (MQ=255)
     aTAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTtgtg                                                                 >  1:1784681/1‑69 (MQ=255)
      tAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTtgtg                                                                 <  1:4192346/68‑1 (MQ=255)
      tAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGc                                                              <  1:2464377/71‑1 (MQ=255)
          cccTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGtt                                                          >  1:4535418/1‑71 (MQ=255)
               cAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTcgcg                                                     >  1:1793936/1‑71 (MQ=255)
               cAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTcg                                                       <  1:2588855/69‑1 (MQ=255)
                aGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTcgcg                                                     <  1:4331295/70‑1 (MQ=255)
                 gTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGtt                                                   <  1:4178063/71‑1 (MQ=255)
                  tttGCCAAACGGCGCTATTTTGTGTTGCAAAGACC‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGtt                                                   <  1:3209798/71‑1 (MQ=255)
                    tGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTtg                                                                   >  1:1573160/1‑52 (MQ=255)
                     gCCAAACGGCGCTATTTTGTGTTGCAAAGACC‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTAt                                                 <  1:2284368/70‑1 (MQ=255)
                     gCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATc                                                <  1:3513384/70‑1 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGt                                           <  1:3213866/71‑1 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCAcc                                             >  1:4230630/1‑69 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCg                                            >  1:2803287/1‑70 (MQ=255)
                         aaCGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGt                                           <  1:1516472/71‑1 (MQ=255)
                          aCGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTg                                          <  1:4412172/71‑1 (MQ=255)
                          aCGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCAc                                              <  1:606181/67‑1 (MQ=255)
                            ggCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTAttttt                                                                    >  1:4424758/1‑43 (MQ=255)
                              cgcTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGatat                                      <  1:3330204/71‑1 (MQ=255)
                               gcTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGatat                                      <  1:2853945/70‑1 (MQ=255)
                                cTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                    <  1:619261/71‑1 (MQ=255)
                                  aTTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATtg                                                                   >  1:1894304/1‑38 (MQ=255)
                                  aTTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTcg                                                       >  1:1352934/1‑50 (MQ=255)
                                   ttttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTg                                                            <  1:4795263/44‑1 (MQ=255)
                                   ttttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGacac                                 >  1:3635062/1‑71 (MQ=255)
                                   ttttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGtt                                                   <  1:4109109/53‑1 (MQ=255)
                                   ttttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGt                                                    >  1:3452131/1‑52 (MQ=255)
                                   ttttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGt                                                    >  1:3097042/1‑52 (MQ=255)
                                   ttttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                                 >  1:548497/1‑55 (MQ=255)
                                   ttttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGacac                                 >  1:251960/1‑71 (MQ=255)
                                    tttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                >  1:4330823/1‑71 (MQ=255)
                                    tttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTATTGTGAGCTg                                                            >  1:3219152/1‑43 (MQ=255)
                                     ttGTGTTGCAAAGA‑CTCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                >  1:4357217/1‑71 (MQ=255)
                                     ttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                    >  1:2910569/1‑66 (MQ=255)
                                     ttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGa                                    >  1:2904758/1‑66 (MQ=255)
                                     ttGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                                >  1:3967624/1‑70 (MQ=255)
                                           tGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACttt                          >  1:823897/1‑70 (MQ=255)
                                            gCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTa                         <  1:3648060/70‑1 (MQ=255)
                                            gCAAAGA‑C‑CCAGTAAGAGTGTATTTTTGTGAGCTTTTTCGCGTTATCACCGTGatat                                      <  1:595538/57‑1 (MQ=255)
                                             cAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTaaa                       >  1:3982237/1‑71 (MQ=255)
                                                    c‑cccGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                                 <  1:707617/39‑1 (MQ=255)
                                                    c‑cccGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                                 <  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)
                                                                      |                                                                    
CAGCAATAAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGA‑C‑CCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAACATTAGATCT  >  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: ‑