Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,252 A→T 47.1% intergenic (+18/+16) add → / ← ydgJ adenosine deaminase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,2520AT47.1% 1.2 / 50.5 34intergenic (+18/+16)add/ydgJadenosine deaminase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (10/8);  new base T (8/8);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

ATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCA  >  NC_000913/1703104‑1703400
                                                                                                                                                    |                                                                                                                                                    
ataCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                                                                                                      <  1:1618828/149‑1 (MQ=255)
                      tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                                                                                                      <  4:1025070/127‑1 (MQ=255)
                      tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                                                                                                      >  3:1025070/1‑127 (MQ=255)
                               ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                                                                                                 >  1:1790981/1‑123 (MQ=255)
                               ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                                                                                                 <  2:974984/123‑1 (MQ=255)
                               ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                                                                                                 <  2:1790981/123‑1 (MQ=255)
                               ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                                                                                                 >  1:974984/1‑123 (MQ=255)
                                                      aGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                                                                             >  4:51949/1‑104 (MQ=255)
                                                      aGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                                                                             <  3:51949/104‑1 (MQ=255)
                                                               gTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGa                                                                                              <  3:1826137/142‑1 (MQ=255)
                                                               gTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGa                                                                                              >  4:1826137/1‑142 (MQ=255)
                                                                           cTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTc                                                                           >  1:499687/1‑149 (MQ=255)
                                                                            tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCa                                                                                                                                                  >  4:506942/1‑77 (MQ=255)
                                                                            tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCa                                                                                                                                                  <  3:506942/77‑1 (MQ=255)
                                                                               ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcca                                                                                                                                           >  2:258426/1‑81 (MQ=255)
                                                                               ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcca                                                                                                                                           <  1:258426/81‑1 (MQ=255)
                                                                                      gcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc                                                                                      >  3:314471/1‑127 (MQ=255)
                                                                                      gcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc                                                                                      <  4:314471/127‑1 (MQ=255)
                                                                                            gAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCt                                                          <  2:499687/149‑1 (MQ=255)
                                                                                                   gcgcACTGCGAGAAAAAGTCGCCGCGTAGTAACAAAATTGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGc                                                   >  3:1175997/1‑149 (MQ=255)
                                                                                                                                         ggatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTc                                    >  2:104219/1‑126 (MQ=255)
                                                                                                                                         ggatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTc                                    <  1:104219/126‑1 (MQ=255)
                                                                                                                                          gatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGc                  >  2:380658/1‑143 (MQ=255)
                                                                                                                                          gatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGc                  <  1:380658/143‑1 (MQ=255)
                                                                                                                                           atggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATc               >  3:876761/1‑145 (MQ=255)
                                                                                                                                           atggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATc               <  4:876761/145‑1 (MQ=255)
                                                                                                                                           atggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCgg                                                >  2:938744/1‑112 (MQ=255)
                                                                                                                                           atggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCgg                                                <  1:938744/112‑1 (MQ=255)
                                                                                                                                              gTGCATTTGTACCATCCATTTTTCAGGCAAGGCACAAAGTCGCGCGAt                                                                                                             >  4:1626581/1‑48 (MQ=255)
                                                                                                                                              gTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGAt                                                                                                             <  3:1626581/48‑1 (MQ=37)
                                                                                                                                              gTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATa        >  2:391550/1‑149 (MQ=255)
                                                                                                                                               tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTg                                                         >  2:1053814/1‑99 (MQ=255)
                                                                                                                                               tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTg                                                         <  1:1053814/99‑1 (MQ=255)
                                                                                                                                               tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTaa                     >  3:349743/1‑135 (MQ=255)
                                                                                                                                               tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTaa                     <  4:349743/135‑1 (MQ=255)
                                                                                                                                                gCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCAt                >  4:1822746/1‑139 (MQ=255)
                                                                                                                                                gCAAATGCACCATCCAGTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCGGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCAt                <  3:1822746/139‑1 (MQ=255)
                                                                                                                                                   aaTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaat                                                                                 >  2:1127218/1‑71 (MQ=255)
                                                                                                                                                   aaTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaat                                                                                 <  1:1127218/71‑1 (MQ=255)
                                                                                                                                                    aTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCa  >  3:841623/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCGCA  >  NC_000913/1703104‑1703400

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