Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,883,501 C→A 52.7% L369L (CTG→CTT cpdB ← 2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,883,5010CA52.7% ‑4.1 / 136.8 112L369L (CTG→CTTcpdB2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (8/45);  new base A (19/40);  total (27/85)
Fisher's exact test for biased strand distribution p-value = 4.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

AATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGA  >  minE/2883435‑2883556
                                                                  |                                                       
aaTGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAg                                                       >  1:572843/1‑69 (MQ=255)
      aTGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAg                                                 <  1:1724597/69‑1 (MQ=255)
         cTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATa                                              >  1:385675/1‑69 (MQ=255)
         cTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATa                                              >  1:3384064/1‑69 (MQ=255)
         cTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATa                                              >  1:2396804/1‑69 (MQ=255)
         cTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGAACGGTCGGATCGTCCTGCACCAGCGCCAGATa                                              >  1:1500595/1‑69 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:1582413/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:3407920/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:2366693/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:3307793/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:2401365/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:1916633/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:864710/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:2512256/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:2988592/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:297009/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:637251/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:86651/68‑1 (MQ=255)
                    cTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACata                                    <  1:1291518/68‑1 (MQ=255)
                     tttttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATAtt                                  <  1:341617/69‑1 (MQ=255)
                     tttttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATAtt                                  <  1:2480401/69‑1 (MQ=255)
                     tttttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATAtt                                  <  1:1630886/69‑1 (MQ=255)
                     tttttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATAtt                                  <  1:294253/69‑1 (MQ=255)
                     tttttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATAtt                                  <  1:194880/69‑1 (MQ=255)
                     tttttGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATAtt                                  <  1:701113/69‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:2907661/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:2197239/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:2269955/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:525115/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1950174/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1271665/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:2578568/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1801021/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1758777/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1433617/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1577309/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:349143/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1510134/67‑1 (MQ=255)
                         tGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCAAAAGCGCCAGATAGCTATACATATTGt                                <  1:813728/67‑1 (MQ=255)
                             cgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAg                          <  1:2682735/69‑1 (MQ=255)
                             cgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAg                          <  1:2043697/69‑1 (MQ=255)
                             cgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAg                          <  1:2315767/69‑1 (MQ=255)
                             cgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAg                          <  1:2084691/69‑1 (MQ=255)
                             cgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAg                          <  1:1398618/69‑1 (MQ=255)
                             cgTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAg                          <  1:587817/69‑1 (MQ=255)
                             cgTTGTTCAACACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAg                          <  1:1541524/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCTGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:1078820/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:3079015/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:2747207/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:2607217/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:2761173/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:3200192/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:827126/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:1959779/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:1103696/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:780197/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:364972/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:2955229/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:2564125/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:1463315/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:1607928/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:559512/69‑1 (MQ=255)
                                tgttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGAtt                       <  1:2141908/69‑1 (MQ=255)
                                 gttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:3087695/68‑1 (MQ=255)
                                 gttCACCACCTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCTGCAGAtt                       <  1:1918473/68‑1 (MQ=255)
                                    caccacCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTAcc                   >  1:3326812/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCtata                                        >  1:1404703/1‑46 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTGCCGa                 >  1:1902842/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCg                  >  1:1050034/1‑68 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:2610771/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:3097981/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:700731/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:1015440/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:2385037/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:2096504/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:200915/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:851025/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:1572698/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:1447775/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:1268044/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:1161926/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGa                 >  1:493176/1‑69 (MQ=255)
                                      ccaccTGCACGGTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTCTCGGCAGATTTACCGa                 >  1:707495/1‑69 (MQ=255)
                                        accTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGAtt               >  1:2474799/1‑69 (MQ=255)
                                             cACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCtt          >  1:2010400/1‑69 (MQ=255)
                                                ggTCGGATCGTCCTGCACAAGCGCCAGATAGCTATACATATTGt                                <  1:1964111/44‑1 (MQ=255)
                                                    ggATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACg   <  1:379210/69‑1 (MQ=255)
                                                    ggATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACg   <  1:3294186/69‑1 (MQ=255)
                                                    ggATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACg   <  1:2873634/69‑1 (MQ=255)
                                                    ggATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACg   <  1:2414561/69‑1 (MQ=255)
                                                    ggATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACg   <  1:1924026/69‑1 (MQ=255)
                                                    ggATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACg   <  1:1823418/69‑1 (MQ=255)
                                                    ggATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGCCTTGCTGACg   <  1:2612837/69‑1 (MQ=255)
                                                      aTCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGa  <  1:3415753/68‑1 (MQ=255)
                                                      aTCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGa  <  1:2556512/68‑1 (MQ=255)
                                                      aTCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGa  <  1:1992588/68‑1 (MQ=255)
                                                      aTCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGa  <  1:857246/68‑1 (MQ=255)
                                                      aTCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGa  <  1:1640878/68‑1 (MQ=255)
                                                      aTCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGa  <  1:1295663/68‑1 (MQ=255)
                                                      aTCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGa  <  1:937819/68‑1 (MQ=255)
                                                                  |                                                       
AATGTAATGCTCGACATACGCTTTTTGCGCGTTGTTCACCACCTGCACGGTCGGATCGTCCTGCACCAGCGCCAGATAGCTATACATATTGTCGGCAGATTTACCGATTGGCTTGCTGACGA  >  minE/2883435‑2883556

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: