Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,606,841 C→T 48.3% A26V (GCC→GTC)  lsrB → AI2 transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,606,8410CT48.3% 1.7 / 36.6 29A26V (GCC→GTC) lsrBAI2 transporter
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base T (9/5);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTTTAAGAAAATCGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGATTTTTTACCAGCG  >  W3110S.gb/1606776‑1606897
                                                                 |                                                        
tCGCTTTAAGAAAATCGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCg                                                        <  1:17659/68‑1 (MQ=255)
        aGAAAATCGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCg                                                 <  1:950362/67‑1 (MQ=255)
        aGAAAATCGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCg                                                 <  1:3087226/67‑1 (MQ=255)
        aGAAAATCGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGAGCAGGCCGCAGAGCg                                                 <  1:1433210/67‑1 (MQ=255)
          aaaaTCGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGTCGCAGAGCGTAt                                              <  1:714852/68‑1 (MQ=255)
               cGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCAtt                                        <  1:3677148/69‑1 (MQ=255)
                     acttaGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTc                                   >  1:2086228/1‑68 (MQ=255)
                        taGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGcc                                                         >  1:2248804/1‑43 (MQ=255)
                         aGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGTCGCAGAGCg                                                 >  1:1305833/1‑50 (MQ=255)
                           cgcTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCAAAc                             <  1:1197372/68‑1 (MQ=255)
                            gcTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTg                           >  1:1088400/1‑69 (MQ=255)
                            gcTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTg                           >  1:1872208/1‑69 (MQ=255)
                                  ggCATTGCCGCAATCTCTATGAATGTGCAGGTCGCAGAGCGTATTg                                            >  1:82401/1‑46 (MQ=255)
                                  ggCATTGCCGCAATCTCTATGAATGTGCAGGTCGCAGAGCGTATTg                                            >  1:1542424/1‑46 (MQ=255)
                                  ggCATTGCCGCAATCTCTATGAATGTGCAGGTCGCAGAGCGTATTg                                            >  1:523641/1‑46 (MQ=255)
                                        gCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGa               >  1:1043575/1‑69 (MQ=255)
                                        gCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGa               >  1:570427/1‑69 (MQ=255)
                                        gCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGa               >  1:468544/1‑69 (MQ=255)
                                        gCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGa               >  1:3426677/1‑69 (MQ=255)
                                        gCCGCAATCTCTAAGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTggg                >  1:1563691/1‑68 (MQ=255)
                                              aTCTCTATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGAtttttt         <  1:3163382/69‑1 (MQ=255)
                                              aTCTCTATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGAtttttt         <  1:160220/69‑1 (MQ=255)
                                                ctctATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCa                                >  1:3479923/1‑44 (MQ=255)
                                                ctctATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCa                                >  1:3481406/1‑44 (MQ=255)
                                                ctctATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCa                                >  1:1798813/1‑44 (MQ=255)
                                                ctctATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCa                                >  1:2199520/1‑44 (MQ=255)
                                                ctctATGAATGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCa                                >  1:18743/1‑44 (MQ=255)
                                                     tGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGATTTTTTACCAGCg  >  1:3271506/1‑69 (MQ=255)
                                                       aaTGTGCAGGTCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGAttttt          <  1:3237222/59‑1 (MQ=255)
                                                             cAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGAttttt          <  1:2893091/53‑1 (MQ=255)
                                                                 |                                                        
TCGCTTTAAGAAAATCGCCTTACTTAGCGCTCTTGGCATTGCCGCAATCTCTATGAATGTGCAGGCCGCAGAGCGTATTGCATTTATTCCCAAACTGGTTGGCGTGGGATTTTTTACCAGCG  >  W3110S.gb/1606776‑1606897

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: