Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,413,585 C→T 8.3% V26V (GTG→GTA ydaO ← predicted C32 tRNA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,413,5850CT 8.3% 181.1 / 6.3 60V26V (GTG→GTAydaOpredicted C32 tRNA thiolase
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base T (2/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTG  >  W3110S.gb/1413523‑1413647
                                                              |                                                              
gagaGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTac                                                         >  1:511075/1‑70 (MQ=255)
gagaGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTac                                                         >  1:1229049/1‑70 (MQ=255)
gagaGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTac                                                         >  1:595685/1‑70 (MQ=255)
gagaGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTac                                                         >  1:535618/1‑70 (MQ=255)
          ccATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGt                                              >  1:966376/1‑71 (MQ=255)
          ccATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGt                                              >  1:1283939/1‑71 (MQ=255)
            atgatgCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGttt                                            >  1:2135406/1‑71 (MQ=255)
            atgatgCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCAACGTTACGACGAAGACGttt                                            >  1:1580248/1‑71 (MQ=255)
                tgCGATCGCCTTCTTCTATCATATTGAAGTCAGCAATGGCTTCGCCTACGTTa                                                          >  1:2017206/1‑53 (MQ=255)
                tgCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCTACGTTa                                                          >  1:528872/1‑53 (MQ=255)
                    aTCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAAt                                     <  1:57508/70‑1 (MQ=255)
                    aTCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAAt                                     <  1:866514/70‑1 (MQ=255)
                          ttcttcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGtt                                >  1:2163668/1‑69 (MQ=255)
                            cttcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAgg                            >  1:1291366/1‑71 (MQ=255)
                            cttcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAgg                            >  1:868185/1‑71 (MQ=255)
                            cttcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAgg                            >  1:395315/1‑71 (MQ=255)
                            cttcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAgg                            >  1:2385578/1‑71 (MQ=255)
                            cttcAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAg                             <  1:513769/70‑1 (MQ=255)
                               cAATCATTTTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAgg                            <  1:557548/68‑1 (MQ=255)
                               cAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAAt                                     <  1:2291947/59‑1 (MQ=255)
                               cAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                         >  1:332267/1‑71 (MQ=255)
                               cAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                         >  1:99134/1‑71 (MQ=255)
                               cAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                         >  1:1008852/1‑71 (MQ=255)
                                aaTCATATTGAAGTCAGCAATGGCTTCGCCTACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGt                        <  1:229051/71‑1 (MQ=255)
                                aaTCATATTGAAGTCAGCAATGGCTTCGCCTACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGt                        <  1:1570840/71‑1 (MQ=255)
                                aaTCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                         <  1:526011/70‑1 (MQ=255)
                                 aTCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTa                       <  1:306963/71‑1 (MQ=255)
                                  tCATATTGAAGTCAGCAATGGCTTCGCCTACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGt                        <  1:38767/69‑1 (MQ=255)
                                  tCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAg                             >  1:718727/1‑64 (MQ=255)
                                  tCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAg                             >  1:1505532/1‑64 (MQ=255)
                                  tCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAg                             >  1:1607870/1‑64 (MQ=255)
                                       ttGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttc                 >  1:1761149/1‑71 (MQ=255)
                                       ttGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGttc                 >  1:1662792/1‑71 (MQ=255)
                                             tCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGt                                        <  1:854688/42‑1 (MQ=255)
                                             tCAGCAATGGCTCCGCCCACGTTACGACGCAGACGTTTTTGt                                        <  1:1874834/42‑1 (MQ=255)
                                              cAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAAt                                     <  1:2364950/44‑1 (MQ=255)
                                                gCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                         <  1:1953828/54‑1 (MQ=255)
                                                gCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                         <  1:1142413/54‑1 (MQ=255)
                                                gCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTg                         <  1:1331940/54‑1 (MQ=255)
                                                gCAATGGCTTTCCCACGTTACGACGCAGACGTTTTTGTAAttt                                   >  1:1471283/1‑43 (MQ=37)
                                                 ctaTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTATTTTGTTCAGGTTg                         <  1:1238530/51‑1 (MQ=255)
                                                 cAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTg                                  <  1:2650113/44‑1 (MQ=255)
                                                  aaTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTa                       >  1:2146746/1‑54 (MQ=255)
                                                  aaTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTa                       >  1:114482/1‑54 (MQ=255)
                                                  aaTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTaa      <  1:1417975/71‑1 (MQ=255)
                                                  aaTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTaa      <  1:2073719/71‑1 (MQ=255)
                                                   aTGGCTTCGCCCACTTTACGACGCAGACGTTTTTGTAATTTGtt                                <  1:2283003/44‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGtt                                <  1:1008019/44‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGtt                                <  1:2544373/44‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGtt                                <  1:1628671/44‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGtt                                <  1:668697/44‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGt                   <  1:2228349/57‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAAt     <  1:1772125/71‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAAt     <  1:1985650/71‑1 (MQ=255)
                                                   aTGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAAt     <  1:2298006/71‑1 (MQ=255)
                                                      gCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATttg  >  1:633307/1‑71 (MQ=255)
                                                      gCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATttg  >  1:115865/1‑71 (MQ=255)
                                                      gCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATttg  >  1:834448/1‑71 (MQ=255)
                                                      gCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATttg  >  1:613367/1‑71 (MQ=255)
                                                      gCTTAGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATttg  >  1:1523484/1‑71 (MQ=255)
                                                              |                                                              
GAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTG  >  W3110S.gb/1413523‑1413647

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: