Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,440,755 T→C 7.1% K165E (AAA→GAA)  cpdB ← 2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,440,7550TC 7.1% 152.5 / 2.9 56K165E (AAA→GAA) cpdB2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base C (2/2);  total (30/26)
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

CCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGAC  >  W3110S.gb/4440692‑4440824
                                                               |                                                                     
ccATTATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGAcc                                                                <  1:1784056/71‑1 (MQ=255)
ccATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGAcc                                                                <  1:3529318/71‑1 (MQ=255)
ccATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGAcc                                                                <  1:1918022/71‑1 (MQ=255)
   tGATTTGTGGTGGCACGACGCCAATATAGACAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACt                                                             >  1:2328792/1‑71 (MQ=255)
     aTTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTc                                                           >  1:2456518/1‑71 (MQ=255)
         gtggtgGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGgtg                                                       >  1:1378498/1‑71 (MQ=255)
            gtggCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCt                                                    >  1:1544540/1‑71 (MQ=255)
               gCACGCCGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCttt                                                  >  1:1462559/1‑70 (MQ=255)
               gCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCttt                                                  <  1:4093764/70‑1 (MQ=255)
               gCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGAACACTTCGGTGTCttt                                                  >  1:3761416/1‑70 (MQ=255)
                 acgacgCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTttaat                                               >  1:1831330/1‑71 (MQ=255)
                  cgacgCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCTGTGTCTttattt                                              <  1:1499493/71‑1 (MQ=255)
                   gacgCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTttaatta                                             <  1:1525317/71‑1 (MQ=255)
                         aaTATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGa                                                                  >  1:128415/1‑44 (MQ=255)
                          atatAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTa                                                                      >  1:3091883/1‑39 (MQ=255)
                          atatAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTCAATTaaa                                           <  1:136049/66‑1 (MQ=255)
                             taGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTttaatt                                              <  1:3157728/60‑1 (MQ=255)
                             taGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATa                                         <  1:828849/65‑1 (MQ=255)
                             taGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGgtg                                    >  1:1570651/1‑70 (MQ=255)
                                ccAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTaaa                                           <  1:3272409/60‑1 (MQ=255)
                                ccAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTaaa                                           <  1:2038791/60‑1 (MQ=255)
                                ccAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTaaa                                <  1:4585613/71‑1 (MQ=255)
                                 cAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTaaa                                >  1:822096/1‑70 (MQ=255)
                                 cAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAAc                               <  1:1720295/71‑1 (MQ=255)
                                 cAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAAc                               <  1:420530/71‑1 (MQ=255)
                                  aaTCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCgg                                                         >  1:2025639/1‑44 (MQ=255)
                                    tCTTCAGCGTCTGTTTTTTTCCGTCTTCATCGACCACTTCGGTGTCttt                                                  <  1:2468624/49‑1 (MQ=255)
                                    tCTTCAGCGTCTGTTTTTTTCCGTCTTCATCGACCACTTCGGTGTCttt                                                  <  1:2474035/49‑1 (MQ=255)
                                     cTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTttaa                                                >  1:2225167/1‑50 (MQ=255)
                                     cTTCAGCGTCTGTTTTTTTCCGTCTTCATCGACCa                                                               >  1:2220476/1‑35 (MQ=255)
                                     cTTCAGCGTCTGTTTTTTTCCGTCTTCATCGACCa                                                               >  1:1464639/1‑35 (MQ=255)
                                      ttCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTaaa                                >  1:1084174/1‑65 (MQ=255)
                                      ttCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTgg                          >  1:124945/1‑71 (MQ=255)
                                      ttCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTgg                          >  1:3565208/1‑71 (MQ=255)
                                      ttCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTg                           >  1:1645326/1‑70 (MQ=255)
                                          gCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTaaa                                           <  1:926018/50‑1 (MQ=255)
                                            gTCTGTTTTTTTCCGTCTTTATCGACCACTTCGgt                                                        <  1:4624962/35‑1 (MQ=255)
                                             tCTGTTTTTTTCCGTCTTTATCGACCTCTTCGGTGTCtt                                                   <  1:3938020/39‑1 (MQ=255)
                                              cTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTgg                  >  1:801294/1‑71 (MQ=255)
                                               tGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGtt                     <  1:4503444/67‑1 (MQ=255)
                                               tGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGt                 >  1:3627539/1‑71 (MQ=255)
                                                gTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGtt                >  1:4107143/1‑71 (MQ=255)
                                                gTTTATTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGt                 >  1:3988641/1‑70 (MQ=255)
                                                 tttttttCCGTCTTTATCGACCACTTCGGTGTCTttaa                                                >  1:2010532/1‑38 (MQ=255)
                                                 tttttttCCGTCTTTATCGACCACTTCGGTGTCTttaa                                                <  1:1341302/38‑1 (MQ=255)
                                                 tttttttCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTg                           <  1:322568/59‑1 (MQ=255)
                                                   tttttCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTg             >  1:423193/1‑71 (MQ=255)
                                                    ttttCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGtt                >  1:391076/1‑67 (MQ=255)
                                                    ttttCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTg             <  1:837189/70‑1 (MQ=255)
                                                      ttCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGc           <  1:820181/70‑1 (MQ=255)
                                                        ccGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGtt                >  1:2622370/1‑63 (MQ=255)
                                                        ccGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGtt                >  1:794406/1‑63 (MQ=255)
                                                          gTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTgg            >  1:3469079/1‑65 (MQ=255)
                                                           tCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAAt     <  1:1409814/71‑1 (MQ=255)
                                                           tCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAAt     <  1:1805371/71‑1 (MQ=255)
                                                            cTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCaa      >  1:751792/1‑69 (MQ=255)
                                                             tttATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGtt                     >  1:3708345/1‑53 (MQ=255)
                                                             tttATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTc               >  1:2142819/1‑59 (MQ=255)
                                                             tttATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGa   >  1:4162724/1‑71 (MQ=255)
                                                             tttATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGa   >  1:1776485/1‑71 (MQ=255)
                                                             tttATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGa   >  1:1047903/1‑71 (MQ=255)
                                                              ttATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGAc  >  1:1342314/1‑71 (MQ=255)
                                                               |                                                                     
CCATGATTTGTGGTGGCACGACGCCAATATAGCCAATCTTCAGCGTCTGTTTTTTTCCGTCTTTATCGACCACTTCGGTGTCTTTAATTAAATACGGTGTAAACATTGGCTGTTTGGTTCTGGCGTCAATGAC  >  W3110S.gb/4440692‑4440824

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: