Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,174,888 C→A 13.8% P218P (CCG→CCT mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,174,8880CA13.8% 103.3 / 3.7 36P218P (CCG→CCTmfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base C (18/13);  new base A (3/2);  total (21/15)
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.65e-01

CGCTTCACTTCGAAGGTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCA  >  W3110S.gb/1174826‑1174957
                                                              |                                                                     
cgcCTCACTTCGAAGGTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTc                                                               <  1:2281843/71‑1 (MQ=255)
 gcTTCACTTCGAAGGTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCg                                                              >  1:466090/1‑71 (MQ=255)
  cTTCACTTCGAAGGTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGt                                                             >  1:2345125/1‑71 (MQ=255)
     cACTTCGAAGGTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGaaa                                                                  >  1:1828670/1‑63 (MQ=255)
     cACTTCGAAGGTATCACGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTAGGAAATTCGTgcg                                                          >  1:3760135/1‑71 (MQ=255)
          cGAAGGTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCggg                                                      >  1:494257/1‑70 (MQ=255)
            aaGGTATCACGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTAGGAAATTCGt                                                             <  1:2961361/61‑1 (MQ=255)
            aaGGTATCACGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTAGGAAATTCGt                                                             <  1:860779/61‑1 (MQ=255)
               gTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGgcagca                                                 <  1:153341/70‑1 (MQ=255)
                  tCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGAtt                                             <  1:1844787/71‑1 (MQ=255)
                    gcgcCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGa                                           >  1:2774547/1‑71 (MQ=255)
                    gcgcCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGa                                           >  1:1813199/1‑71 (MQ=255)
                          cTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGc                                       >  1:253612/1‑69 (MQ=255)
                          cTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGc                                       >  1:658405/1‑69 (MQ=255)
                            ggCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTc                                                               >  1:702314/1‑43 (MQ=255)
                            ggCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTc                                                               >  1:156092/1‑43 (MQ=255)
                             gatGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTc                                    <  1:130557/67‑1 (MQ=255)
                             gCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTc                                    <  1:770083/69‑1 (MQ=255)
                             gCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCt                                   <  1:1128110/70‑1 (MQ=255)
                                     aCAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTtcctc                           >  1:3245665/1‑70 (MQ=255)
                                        gTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTg                                            <  1:1088933/50‑1 (MQ=255)
                                         ttCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGgcagca                                                 <  1:662845/44‑1 (MQ=255)
                                         ttCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGcagcag                                                <  1:3142819/45‑1 (MQ=255)
                                           cAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGt                     >  1:3252121/1‑70 (MQ=255)
                                           cAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTg                    >  1:969550/1‑71 (MQ=255)
                                           cAATTGCCGCTTTATCGGTAGGAAATTCGTGCGCGGGAAGCAGATTGATCTCTTCTACTTCCTCCAGCGTg                    >  1:511917/1‑71 (MQ=255)
                                           cAATTGCCGCTTTATCGGTAGGAAATTCGTGCGCGGGAAGCAGATTGATCTCTTCTACTTCCTCCAGCGTg                    >  1:130270/1‑71 (MQ=255)
                                            aaTTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTg                    <  1:2213151/70‑1 (MQ=255)
                                              ttGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGAt                                          >  1:2309705/1‑46 (MQ=255)
                                              ttGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGAt                                          >  1:2220202/1‑46 (MQ=255)
                                               tGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTtcctcc                          >  1:1594337/1‑61 (MQ=255)
                                                gCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTa                                  >  1:259809/1‑52 (MQ=255)
                                                 ccGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCgctg               <  1:179395/70‑1 (MQ=255)
                                                       tATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTc         <  1:283131/70‑1 (MQ=255)
                                                       tATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTc         <  1:2656492/70‑1 (MQ=255)
                                                        aTCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGa       >  1:3413133/1‑71 (MQ=255)
                                                             tcggAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCa  <  1:3412274/71‑1 (MQ=255)
                                                              |                                                                     
CGCTTCACTTCGAAGGTATCGCGCCACTGGCTGCGGAACAGTTCAATTGCCGCTTTATCGGTCGGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCA  >  W3110S.gb/1174826‑1174957

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: