Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,346,872 C→T 41.9% V529I (GTT→ATT)  gmr ← modulator of Rnase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,346,8720CT41.9% 24.0 / 49.4 43V529I (GTT→ATT) gmrmodulator of Rnase II stability
Reads supporting (aligned to +/- strand):  ref base C (12/13);  new base T (9/9);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGA  >  W3110S.gb/1346805‑1346939
                                                                   |                                                                   
cTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATc                                                                  <  1:2747878/71‑1 (MQ=255)
cTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATc                                                                  <  1:473216/71‑1 (MQ=255)
  aaaTTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCta                                                                >  1:2323666/1‑71 (MQ=255)
  aaaTTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCta                                                                >  1:1602908/1‑71 (MQ=255)
   aaTTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAATATCta                                                                <  1:736329/70‑1 (MQ=255)
   aaTTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCta                                                                <  1:188203/70‑1 (MQ=255)
         tgAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACCCTCTGTCAGTTCAACATCTATAggg                                                           >  1:1254889/1‑69 (MQ=255)
             tAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCta                                                                <  1:1936660/60‑1 (MQ=255)
                      gTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTg                                            <  1:1866675/71‑1 (MQ=255)
                      gTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTg                                            <  1:2325698/71‑1 (MQ=255)
                             ttCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATc                                                                  <  1:2122510/42‑1 (MQ=255)
                             ttCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCta                                                                <  1:1831083/44‑1 (MQ=255)
                             ttCATCATTCTCAATCAGACAAATCTCTGTCAGTTCAACATCTATAGGGCAg                                                        <  1:2780451/52‑1 (MQ=255)
                              tcatcaTTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCaa                                                 >  1:1823902/1‑58 (MQ=255)
                               catcatTCTCAATCAGACAACTCTCTGTCAGTTCTATATCTATAGGGCAGTATTCAAAATTGAGTTCCTg                                    <  1:2028606/70‑1 (MQ=255)
                               catcatTCTCAATCAGACAACTCTCTGTCAGTTCAATATCta                                                                >  1:1862820/1‑42 (MQ=255)
                               catcatTCTCAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTGAGTTCCTg                                    <  1:2565786/70‑1 (MQ=255)
                               catcatTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTg                                    <  1:62618/70‑1 (MQ=255)
                               cataatTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTg                                    >  1:593112/1‑70 (MQ=255)
                                atcatTCTCAATCAGACAACTCTCTGTCAGTTCAATATCta                                                                <  1:1266720/41‑1 (MQ=255)
                                 tcatTCTCAATCAGACAACTCTCTGTCAGTTCAAc                                                                     <  1:212382/35‑1 (MQ=255)
                                 tcatTCTCAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTGAGTTCCTg                                    >  1:1194915/1‑68 (MQ=255)
                                 tcatTCTCAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTGAGTTCCTGgag                                 <  1:437015/71‑1 (MQ=255)
                                  catTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAg                                                        <  1:1471953/47‑1 (MQ=255)
                                      ctcAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTAt                                                     >  1:1246008/1‑46 (MQ=255)
                                      ctcAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTAt                                                     >  1:348186/1‑46 (MQ=255)
                                      ctcAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTAt                                                     >  1:1429270/1‑46 (MQ=255)
                                      ctcAATCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTAt                                                     >  1:1620567/1‑46 (MQ=255)
                                           tCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTATTCAAAAtt                                             <  1:894875/49‑1 (MQ=255)
                                           tCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTg                                            <  1:1321067/50‑1 (MQ=255)
                                           tCAGACAACTCTCTGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTg                                            <  1:295602/50‑1 (MQ=255)
                                             aGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAAt                                              >  1:803919/1‑46 (MQ=255)
                                                    tctctGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGt               >  1:2158990/1‑70 (MQ=255)
                                                      tctGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGa             >  1:2010497/1‑70 (MQ=255)
                                                      tctGTCAGTTCAATATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGa             >  1:504412/1‑70 (MQ=255)
                                                      tctGTCAGTTCAACATCTATAGGGCAGTATTCAAAAt                                              >  1:898473/1‑37 (MQ=255)
                                                      tctGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCt                                     >  1:1351894/1‑46 (MQ=255)
                                                            aGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCa                      >  1:1575112/1‑55 (MQ=255)
                                                            aGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAg                     >  1:213063/1‑56 (MQ=255)
                                                              ttCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGaa            <  1:1355963/63‑1 (MQ=255)
                                                                cAATATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGa  >  1:1598317/1‑71 (MQ=255)
                                                                cAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGga                                  >  1:225577/1‑39 (MQ=255)
                                                                cAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTg   <  1:240303/70‑1 (MQ=255)
                                                                 aaCATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTg   <  1:506820/69‑1 (MQ=255)
                                                                   |                                                                   
CTAAATTGTTGAATAACAGACAGTGCCAGTTCATCATTCTCAATCAGACAACTCTCTGTCAGTTCAACATCTATAGGGCAGTATTCAAAATTGAGTTCCTGGAGAACCTGTTTCAGGGCGGTGAAAATGGTTTGA  >  W3110S.gb/1346805‑1346939

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: