Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,017,703 T→C 12.9% N1175D (AAC→GAC)  rhsB ← rhsB element core protein RshB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,017,7030TC12.9% 158.8 / 14.3 62N1175D (AAC→GAC) rhsBrhsB element core protein RshB
Reads supporting (aligned to +/- strand):  ref base T (32/22);  new base C (3/5);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGC  >  W3110S.gb/4017634‑4017760
                                                                     |                                                         
catcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTg                                                          >  1:4319726/1‑71 (MQ=34)
 atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGtt                                                           <  1:3480034/69‑1 (MQ=35)
 atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGc                                                         >  1:3627337/1‑71 (MQ=34)
 atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGc                                                         <  1:1752779/71‑1 (MQ=34)
 atcatACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGc                                                         <  1:4405417/71‑1 (MQ=34)
    atACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGccc                                                       >  1:4517920/1‑70 (MQ=35)
     tACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCa                                                     >  1:2578004/1‑71 (MQ=34)
          ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGccc                                                       >  1:3718923/1‑64 (MQ=35)
          ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcat                                                >  1:3985109/1‑71 (MQ=34)
          ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtca                                                 <  1:4087639/70‑1 (MQ=35)
          ctgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCACtc                                                  <  1:2328411/69‑1 (MQ=25)
           tgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcat                                                >  1:3631792/1‑70 (MQ=35)
            gcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca                                              >  1:678614/1‑71 (MQ=35)
             ccccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca                                              >  1:870433/1‑70 (MQ=37)
             ccccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtca                                                 <  1:3976256/67‑1 (MQ=35)
              cccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCata                                            <  1:4454545/71‑1 (MQ=39)
               ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                            >  1:1092352/1‑70 (MQ=255)
                  gCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGcccc                                                      <  1:2886092/57‑1 (MQ=255)
                      ccGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTgcg                                     >  1:3466368/1‑70 (MQ=37)
                           tAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac                               >  1:1792733/1‑71 (MQ=35)
                           tAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac                               >  1:4471175/1‑71 (MQ=35)
                           tAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac                               >  1:7008/1‑71 (MQ=35)
                            aaGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACg                              >  1:1515996/1‑71 (MQ=38)
                              gctgctGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt                            >  1:1029605/1‑71 (MQ=40)
                              gctgctGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt                            <  1:3989112/71‑1 (MQ=40)
                              gctgctGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCt                            <  1:1847099/71‑1 (MQ=255)
                               ctgcttCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTg                           <  1:2567677/71‑1 (MQ=255)
                                tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcacca                                >  1:786884/1‑65 (MQ=35)
                                 gctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGtt                         >  1:195973/1‑71 (MQ=255)
                                  ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCTTCATATTCTGCGCACCACGCTGt                          <  1:42191/69‑1 (MQ=255)
                                   tgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGc                       >  1:3718220/1‑71 (MQ=255)
                                    gcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCt                                        >  1:1176722/1‑53 (MQ=35)
                                    gcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcaccac                               <  1:455822/62‑1 (MQ=37)
                                        cTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCt                            <  1:4036418/61‑1 (MQ=255)
                                        cTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCtt                  >  1:3893048/1‑71 (MQ=40)
                                         tGATGCGGGTTCTCTTCATTCAGCAGGTTG‑‑CCATTCATCATATTc                                         >  1:2079367/1‑45 (MQ=11)
                                          gATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATtcatca                                              >  1:4422842/1‑41 (MQ=32)
                                          gATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTcc                >  1:3352164/1‑71 (MQ=39)
                                           aTGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCg               <  1:1778066/71‑1 (MQ=39)
                                            tGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTcc                <  1:704746/69‑1 (MQ=40)
                                            tGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGt              <  1:1970642/71‑1 (MQ=39)
                                             gCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGccc                     <  1:1983706/63‑1 (MQ=40)
                                              cGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGc            >  1:1585624/1‑71 (MQ=39)
                                               gggTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCt                                        >  1:1924373/1‑42 (MQ=38)
                                                ggTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                            <  1:3801000/37‑1 (MQ=255)
                                                ggTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                            <  1:1766175/37‑1 (MQ=255)
                                                 gTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATAtt                                          <  1:609393/38‑1 (MQ=40)
                                                 gTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGt                          >  1:210751/1‑54 (MQ=255)
                                                 gTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGt                          >  1:205019/1‑54 (MQ=255)
                                                  ttCTCTTCATTCAGCAGGTTGCCCCATTCATCata                                            <  1:3240154/35‑1 (MQ=40)
                                                  ttCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGcacca                                <  1:3888224/47‑1 (MQ=38)
                                                   tctcTTCATTCAGCAGGTTGCCCCATTCATCATATTCTg                                       >  1:3102495/1‑39 (MQ=40)
                                                   tctcTTCATTCAGCAGGTTGCCCCATTCATCATATTCTg                                       >  1:2718997/1‑39 (MQ=40)
                                                    ctctTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgatg       <  1:3704978/70‑1 (MQ=40)
                                                    ctctTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgat        >  1:252781/1‑69 (MQ=40)
                                                       ttcattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCtgatga      <  1:563/68‑1 (MQ=37)
                                                       ttcattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcg   >  1:1505312/1‑71 (MQ=35)
                                                        tcattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGc                             >  1:1334903/1‑44 (MQ=255)
                                                         cattcaGCAGGTTGCCCCATTCCTCATATTCTGCGca                                   >  1:1014866/1‑37 (MQ=255)
                                                         cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcg   <  1:3246588/69‑1 (MQ=35)
                                                         cattcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgc    <  1:2265494/68‑1 (MQ=35)
                                                         cattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcgc  >  1:157131/1‑70 (MQ=25)
                                                         cattcaGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAgcgc  <  1:1142593/70‑1 (MQ=25)
                                                                     |                                                         
CATCATACTGCTGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGC  >  W3110S.gb/4017634‑4017760

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: