Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 352,772 G→A 10.8% T297T (ACG→ACA codA → cytosine/isoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273352,7720GA10.8% 93.4 / 4.8 37T297T (ACG→ACAcodAcytosine/isoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (16/17);  new base A (2/2);  total (18/19)
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.06e-01

AAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCGGCTGAAAATGG  >  CP009273/352500‑353035
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                       
aaaTACGACCGTCTCAGCGACGTTCACTGGGCGGAGATCGAGGACGAGCAGTCGCGGTTTGTCGAAACCGGTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACCCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCAGCAcgcg                                                                                                                                                                                                                                                                       <  1:157333/275‑1 (MQ=255)
                                   gaTCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGACAACCAGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGACGCGGCCTCACGCGCGTTAACGAGATGCTGGAGTCCGGCATTAACGCCt                                                                                                                                                                                                                                    >  1:316173/1‑275 (MQ=255)
                                      ctgtGACGATCTGTCGTGCTTTGTCGAAACCGTTGCTGGCCTGGCGCTACATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCt                                                                                                                                                                                                                                 <  1:300251/272‑1 (MQ=255)
                                              aGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGCGTCTACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTCACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTAGAGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGCCTGCATTGGTCa                                                                                                                                                                                                                          >  1:244811/1‑274 (MQ=255)
                                                                   ccGTTGCTGCCCTGGAGAACCATGAGGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTACGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTgg                                                                                                                                                                                                    <  1:127870/275‑1 (MQ=255)
                                                                                                   ggCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTANCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGt                                                                                                                                                                    >  1:416056/1‑275 (MQ=255)
                                                                                                                    cagccaCACCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACACGCGtt                                                                                                                                                                                                                                                                   >  2:402941/1‑163 (MQ=255)
                                                                                                                        cacaCCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACACGCGTTAAag                                                                                                                                                                                                                                                               <  1:402954/163‑1 (MQ=255)
                                                                                                                                    ttGCACTCCGCTAACGTGGCTTCTACCTCACTCCTTTTGCGCTTGCTGAAAATGTCCGGTTTTAACTTTGTCGCCAACCCGCTGGTCAATATGCAGCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGAt                                                                                                                                   <  2:444599/274‑1 (MQ=255)
                                                                                                                                     tGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCtt                                                                                                                                                                                                                                >  2:229375/1‑181 (MQ=255)
                                                                                                                                         cTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGt                                                                                                                                                                                                                            <  1:229382/181‑1 (MQ=255)
                                                                                                                                            cTATAACGGGGCGTATACCTCACGCCTGATCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCAGCATTAACGTCTGCTTTGGTCACGANGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTATGACAGTTGATGGGGTACg                                                                                                                           >  1:370735/1‑275 (MQ=255)
                                                                                                                                               taACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATAtg                                                                                                                                                                            >  2:75375/1‑223 (MQ=255)
                                                                                                                                                aaCGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTc                                                                                                                                                                                                                                                  <  1:15364/152‑1 (MQ=255)
                                                                                                                                                   ggggCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATAtgctgc                                                                                                                                                                        <  1:75377/223‑1 (MQ=255)
                                                                                                                                                      gCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATATCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTaa                                                                                                                 >  1:277773/1‑275 (MQ=255)
                                                                                                                                                                cACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCtt                                                                                                                                                                                                                                >  2:98216/1‑154 (MQ=255)
                                                                                                                                                                    ccTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGt                                                                                                                                                                                                                            <  1:98218/154‑1 (MQ=255)
                                                                                                                                                                        ttCCTCTTGCTGAAAATGTCCGTTTTTAACTTTGTCGCCAACCCGCTGGTCCATATTCATCTGCAAGGACGTTTCGATCCGNAGCCAAAANGTCGCGGCATCCTGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGGTCCGCTGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATCGGCTACGTGCAGATTAACGATGGCCTGAATTTaa                                                                                                <  1:329274/274‑1 (MQ=255)
                                                                                                                                                                        ttCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACACGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTa                                                                                                                             >  1:147900/1‑245 (MQ=255)
                                                                                                                                                                            gCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCAcg                                                                                                                                                                                                                                                                         <  1:172078/101‑1 (MQ=255)
                                                                                                                                                                            gCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACACGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACggg                                                                                                                         <  2:147894/245‑1 (MQ=255)
                                                                                                                                                                             cTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGACCGGCGAATATGCTGCAAg                                                                                                                                                                     >  2:299793/1‑200 (MQ=255)
                                                                                                                                                                                gCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATAtg                                                                                                                                                                            >  1:70481/1‑190 (MQ=255)
                                                                                                                                                                                 cTGAAAATGGCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGtgct                                                                                                                                                                 <  1:299802/200‑1 (MQ=255)
                                                                                                                                                                                    aaaaTGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATAtgctgc                                                                                                                                                                        <  2:70479/190‑1 (MQ=255)
                                                                                                                                                                                        ggaaatgTATTAATTTTGGCGCCAACCCGCTGGTCAATATTCGTCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCa                                                                               <  2:416043/269‑1 (MQ=255)
                                                                                                                                                                                           ggtgttttAACTTTNCCCCCCACCAGGTGGTCAATATTCATCTGCGACGTCGTTTCTATACGTATCCAAAACGTCACTGCATTACGCGCTTTAAATAGGTTCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTTCAAGTGCTGCATATGAGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATTGCCTGAATTTAATCACCCACCACAGCGCAAgg                                                                            <  2:370722/269‑1 (MQ=255)
                                                                                                                                                                                                         gTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTc                                                                                                                                                                                                                                     >  1:108951/1‑108 (MQ=255)
                                                                                                                                                                                                             ccAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCt                                                                                                                                                                                                                                 <  2:108946/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTc                                                                                                                                                                                                                                     >  1:42654/1‑48 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCtt                                                                                                                                                                                                              >  1:298514/1‑71 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGTGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCGGCTGAAAATgg  >  1:30044/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTg                                                                                                                                                    >  1:135468/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTa                                                                                                                             >  1:178186/1‑152 (MQ=255)
                                                                                                                                                                                                                                                                         gcagcaCGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACggg                                                                                                                         <  2:178180/152‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCt                                                                                                                                                                                                                                 <  2:42652/48‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGAt                                                                                                                                                                                                          <  2:298505/71‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATg                                                                                                                                                <  2:135462/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                       
AAATACGACCGTCTCATCGACGTTCACTGTGATGAGATCGATGACGAGCAGTCGCGCTTTGTCGAAACCGTTGCTGCCCTGGCGCACCATGAAGGCATGGGCGCGCGAGTCACCGCCAGCCACACCACGGCAATGCACTCCTATAACGGGGCGTATACCTCACGCCTGTTCCGCTTGCTGAAAATGTCCGGTATTAACTTTGTCGCCAACCCGCTGGTCAATATTCATCTGCAAGGACGTTTCGATACGTATCCAAAACGTCGCGGCATCACGCGCGTTAAAGAGATGCTGGAGTCCGGCATTAACGTCTGCTTTGGTCACGATGATGTCTTCGATCCGTGGTATCCGCTGGGAACGGCGAATATGCTGCAAGTGCTGCATATGGGGCTGCATGTTTGCCAGTTGATGGGCTACGGGCAGATTAACGATGGCCTGAATTTAATCACCCACCACAGCGCAAGGACGTTGAATTTGCAGGATTACGGCATTGCCGCCGGAAACAGCGCCAACCTGATTATCCTGCCGGCTGAAAATGG  >  CP009273/352500‑353035

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: