Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,434,349 T→A 100% 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,3490TA81.1% 84.3 / 15.2 37intergenic (+112/‑197)fklB/cycAFKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (7/0);  new base A (10/20);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.89e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAA  >  W3110S.gb/4434288‑4434407
                                                              |                                                          
aaaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCt                                                   >  1:899866/1‑71 (MQ=255)
 aaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                  >  1:807829/1‑71 (MQ=255)
 aaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                  >  1:219373/1‑71 (MQ=255)
 aaCCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTg                                                  >  1:217121/1‑71 (MQ=255)
   cccTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATtg                                                         >  1:1830331/1‑62 (MQ=255)
            ttGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTAt                                       >  1:1341775/1‑71 (MQ=255)
             tGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTcg                                             >  1:1629753/1‑64 (MQ=255)
                    cGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                               <  1:1142345/71‑1 (MQ=255)
                    cGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                               <  1:867573/71‑1 (MQ=255)
                        gcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                               <  1:1304568/67‑1 (MQ=255)
                        gcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                               <  1:1387459/67‑1 (MQ=255)
                        gcTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGa                               <  1:1676389/67‑1 (MQ=255)
                            ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGc                                                    <  1:10088/42‑1 (MQ=255)
                            ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGc                                                    <  1:124547/42‑1 (MQ=255)
                            ttttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGc                                                    <  1:1375813/42‑1 (MQ=255)
                              ttGTGTTGCAAAGACCCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                      >  1:1198019/1‑71 (MQ=255)
                              ttGTGTTGCAAAGACCCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                      >  1:280163/1‑71 (MQ=255)
                              ttGTGTTGCAAAGACCCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                      >  1:450784/1‑71 (MQ=255)
                              ttGTGTTGCAAAGACCCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                      >  1:967817/1‑71 (MQ=255)
                              ttGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                      >  1:1224264/1‑70 (MQ=255)
                              ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGacac                       >  1:617809/1‑69 (MQ=255)
                              ttGTGTTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcact                      >  1:7023/1‑70 (MQ=255)
                                gtgtTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcactcac                   >  1:1197399/1‑71 (MQ=255)
                                 tgtTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcactcac                   <  1:317967/70‑1 (MQ=255)
                                 tgtTGCAAAGACCCCGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGAcactcac                   <  1:315758/70‑1 (MQ=255)
                                                cGTAAGCGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAAtt    <  1:454850/71‑1 (MQ=255)
                                                    aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGata                             >  1:1087769/1‑42 (MQ=255)
                                                    aGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGata                             >  1:1813661/1‑42 (MQ=255)
                                                      cGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTaa  <  1:1708536/67‑1 (MQ=255)
                                                      cGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTaa  <  1:1519238/67‑1 (MQ=255)
                                                      cGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTaa  <  1:1379394/67‑1 (MQ=255)
                                                      cGTGTATTATTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTaa  <  1:822424/67‑1 (MQ=255)
                                                         gtATTATTGTGAGCTGTTTCGCGTTATCACCGTGatat                            <  1:1768252/38‑1 (MQ=255)
                                                         gtATTATTGTGAGCTGTTTCGCGTTATCACCGTGatat                            <  1:1754005/38‑1 (MQ=255)
                                                         gtATTATTGTGAGCTGTTTCGCGTTATCACCGTGatat                            <  1:1610192/38‑1 (MQ=255)
                                                         gtATTATTGTGAGCTGTTTCGCGTTATCACCGTGatat                            <  1:896978/38‑1 (MQ=255)
                                                         gtATTATTGTGAGCTGTTTCGCGTTATCACCGTGatat                            <  1:1196385/38‑1 (MQ=255)
                                                          tattatTGTGAGCTGTTTCGCGTTATCAGCGTGATATGACACTCACTTTAAACATaa        >  1:1255835/1‑57 (MQ=39)
                                                          tattatTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaa        >  1:666945/1‑57 (MQ=255)
                                                          tattatTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaa        >  1:1392242/1‑57 (MQ=255)
                                                          tattatTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaa        >  1:719704/1‑57 (MQ=255)
                                                          tattatTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATaa        >  1:725824/1‑57 (MQ=255)
                                                          tattatTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATa         >  1:1284462/1‑56 (MQ=255)
                                                              |                                                          
AAACCCTTCAGTTTGCCAAACGGCGCTATTTTGTGTTGCAAAGACCCCGTAAGCGTGTATTTTTGTGAGCTGTTTCGCGTTATCACCGTGATATGACACTCACTTTAAACATAAAATTAA  >  W3110S.gb/4434288‑4434407

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: