Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,424,794 G→T 7.5% P530Q (CCG→CAG)  cpdB ← 2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,424,7940GT 7.5% 139.4 / 3.9 54P530Q (CCG→CAG) cpdB2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (29/21);  new base T (2/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

CACGCTGCCAGCACCGAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCG  >  CP009273/4424662‑4424930
                                                                                                                                    |                                                                                                                                        
cACGCTGCCAGCACCGAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTTGTGGCAAACAGGAACATTGAGTTCGGATCAATCGGCTTGca                                                                                                                                   >  1:1314997/1‑139 (MQ=255)
   gCTGCCAGCACCGAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCAtt                                                                                                                                >  2:1344003/1‑140 (MQ=255)
            aCCGAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAg                                                                                                                       <  2:1248777/140‑1 (MQ=255)
            aCCGAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAg                                                                                                                       <  1:761049/140‑1 (MQ=255)
               gAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCTGCTTGCCATTAAAGGTCAGGtt                                                                                                                    >  2:263540/1‑140 (MQ=255)
               gAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCTGCTTGCCATTAAAGGTCAGGtt                                                                                                                    >  2:94171/1‑140 (MQ=255)
               gAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGtt                                                                                                                    >  2:786416/1‑140 (MQ=255)
                 gCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCt                                                                                                                  >  2:399776/1‑140 (MQ=255)
                 gCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCt                                                                                                                  >  1:1140824/1‑140 (MQ=255)
                 gCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCt                                                                                                                  >  2:907646/1‑140 (MQ=255)
                 gCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCt                                                                                                                  >  2:399788/1‑140 (MQ=255)
                 gCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCt                                                                                                                  >  2:399783/1‑140 (MQ=255)
                 gCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCt                                                                                                                  >  2:205737/1‑140 (MQ=255)
                 gCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCt                                                                                                                  >  1:513927/1‑140 (MQ=255)
                          aTCCGGTGAAGCAAAAGCGAGAAGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCttt                                                                                                         <  2:1178734/140‑1 (MQ=255)
                             cGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCg                                                                                                      >  2:897527/1‑140 (MQ=255)
                             cGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCg                                                                                                      >  2:487590/1‑140 (MQ=255)
                             cGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCg                                                                                                      >  2:462605/1‑140 (MQ=255)
                               gTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCa                                                                                                    <  2:911368/140‑1 (MQ=255)
                               gTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCa                                                                                                    <  2:911367/140‑1 (MQ=255)
                               gTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCa                                                                                                    >  2:640743/1‑140 (MQ=255)
                               gTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCa                                                                                                    >  2:1066303/1‑140 (MQ=255)
                                tGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCAt                                                                                                   <  2:492766/140‑1 (MQ=255)
                                tGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCAt                                                                                                   <  2:263659/140‑1 (MQ=255)
                                         aGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTaa                                                                                          <  2:690259/140‑1 (MQ=255)
                                              tatGGCTGTCGCCCGTACCGGCAAATTTTCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATGCTTTCCGCATTGGAATTAatcatc                                                                                     <  2:1291040/140‑1 (MQ=255)
                                               atGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAAGCCTTTCCGCATTGGCATTAATCATCt                                                                                    >  2:7755/1‑140 (MQ=255)
                                                      tCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTc                                                                             >  2:1098646/1‑140 (MQ=255)
                                                      tCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTc                                                                             >  2:836082/1‑140 (MQ=255)
                                                      tCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTc                                                                             >  2:576506/1‑140 (MQ=255)
                                                      tCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTATTAATCCTTTCCGCATTGGAATTAATCATCTGGCACTc                                                                             >  2:1126277/1‑140 (MQ=255)
                                                             tACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCa                                                                      >  2:150598/1‑140 (MQ=255)
                                                             tACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGGTCGGATCAATCGCCTTGCCATAAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTCGCACTCGCCGTCa                                                                      >  2:762528/1‑140 (MQ=255)
                                                              aCCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCAt                                                                     <  1:795841/140‑1 (MQ=255)
                                                              aCCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCAt                                                                     <  1:1175640/140‑1 (MQ=255)
                                                                         tGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACgggcgggc                                                          <  2:782772/140‑1 (MQ=255)
                                                                              ccgTAAGCGCGATAGTTCTTGGTGGCAACAAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATGCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGt                                                     <  2:317298/140‑1 (MQ=255)
                                                                                          tAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTg                                         <  1:1066303/140‑1 (MQ=255)
                                                                                          tAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTg                                         <  1:640743/140‑1 (MQ=255)
                                                                                                       cAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCa                            >  1:1241194/1‑140 (MQ=255)
                                                                                                           cAGGAACATGGCTTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCaa                        <  2:1314997/140‑1 (MQ=255)
                                                                                                           cAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCaa                        <  2:668836/140‑1 (MQ=255)
                                                                                                                  aTGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCa                                                                                >  1:1256030/1‑77 (MQ=255)
                                                                                                                  aTGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCa                                                                                <  2:1256030/77‑1 (MQ=255)
                                                                                                                  aTGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCACTCTGATAATTCACACCATCAATAACATc                 >  2:93618/1‑140 (MQ=255)
                                                                                                                  aTGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATc                 <  2:1183034/140‑1 (MQ=255)
                                                                                                                  aTGGCGTTCGGATCAATCGGCTTGCCATTAAAGCTCAGGTTCTTAATCCTTTCCGCATTCGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATc                 >  2:514148/1‑140 (MQ=255)
                                                                                                                  aTGGCGTTCGGAACAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGCCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATc                 <  2:760269/140‑1 (MQ=255)
                                                                                                                    ggCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCaa               >  1:323288/1‑140 (MQ=255)
                                                                                                                    ggCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCaa               >  1:807213/1‑140 (MQ=255)
                                                                                                                                aaTCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGc   >  2:855378/1‑140 (MQ=255)
                                                                                                                                 aTCTGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCCTCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCg  <  1:263540/140‑1 (MQ=255)
                                                                                                                                 aTCTGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCg  <  1:94171/140‑1 (MQ=255)
                                                                                                                                 aTCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCg  <  1:786416/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                        
CACGCTGCCAGCACCGAGCGGTTCTCATCCGGTGAAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCG  >  CP009273/4424662‑4424930

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑