Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,154,691 G→A Q53* (CAG→TAG)  argE ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,154,6910GA92.6% 74.2 / ‑0.7 27Q53* (CAG→TAG) argEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (9/16);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

GGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCAC  >  NC_000913/4154444‑4154936
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
ggTTTTTTCAGTTTCGTGACGTCGACAACGCGTAGCGCATCAAGAATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTacac                                                                                                                                                                                                                                                    <  1:30922/251‑1 (MQ=255)
        cAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCcaca                                                                                                                                                                                                                                             <  1:142632/250‑1 (MQ=255)
                   cGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCaa                                                                                                                                                                                                                                  >  1:191870/1‑250 (MQ=255)
                      cGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCACCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAActct                                                                                                                                                                                                                              >  1:126893/1‑251 (MQ=255)
                         cATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTa                                                                                                                                                                                                                           <  1:149606/251‑1 (MQ=255)
                                      cacAAGGATAACCCAAAAAAAGCTTTGCATGTCCGCGGGGCCTAAGCCGAAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGATTCGCCCCTCCAGGGCCCGTCATCAAATGGCCCCGTATCGTTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCca                                                                                                                                                                                                              <  2:199702/249‑1 (MQ=255)
                                      aacAAGGATAAACACAAAAAGGCCTTTCAGGTCGGCGGTGCCTAAGCCTTAAAGCTTGCCGTCACGCTCCGCCAGTGTAAACGGAGCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTACCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATAGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCTCATTGAAGCCCAAATCTTTAAACCAGTACGCca                                                                                                                                                                                                              <  1:201901/249‑1 (MQ=255)
                                       tCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTACCTAAGACGTAAAGCTTGACGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGGCATCAAATGGCACCGTAACGGCATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGGAGGGAGGTCCTGGAACCGGCTAAACTTCCCCATTGAGACCCAACACTTTAAACCAGtacacca                                                                                                                                                                                                              >  2:308/1‑244 (MQ=255)
                                             aaaaaCGCAAAAAAGCCTTCCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGCCACCAAATGGCCCCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCcaca                                                                                                                                                                                                                                             <  2:92000/211‑1 (MQ=255)
                                             aTAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCcaca                                                                                                                                                                                                                                             >  1:92000/1‑213 (MQ=255)
                                               aaaCGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCCACAAGCCGCCAGACCCCTGTCCGCTACTTGCCAGCATATTGAATTTGTTTCGAGTTCCTGGAACAGGCTACCCTTCCACATTGAAGCCCAAATCTTTAAACCAGTCAGCCAGCagagcgt                                                                                                                                                                                                     >  2:102904/1‑248 (MQ=255)
                                                   gCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTa                                                                                                                                                                                                  >  2:162015/1‑250 (MQ=255)
                                                               ttCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCt                                                                                                                                                                                     >  1:176820/1‑251 (MQ=255)
                                                                    gccggcgtTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCTTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGGCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGAt                                                                                                                                                                                <  1:89595/249‑1 (MQ=255)
                                                                        gtgggtcttaacccttacgtcTTGCCGTCATGCTCCGTCAGTGTAAACGGACCGCGCGCCCAGCGACCGCAATCAACAGGCCCCGGACGGTTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAg                                                                                                                                                                            <  2:91611/231‑1 (MQ=255)
                                                                          ggTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTg                                                                                                                                                                          <  1:262860/251‑1 (MQ=255)
                                                                                                                                    catcGACCGTCATCAAATGGCCCCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATaaa                                                                                                                <  1:135967/248‑1 (MQ=255)
                                                                                                                                                 aaaatggacccgattCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGTAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTg                                                                                                   <  2:292555/237‑1 (MQ=255)
                                                                                                                                                               ccGGTATGCCCGGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGa                                                                                     <  1:213248/250‑1 (MQ=255)
                                                                                                                                                                          ccGCCAGCAACAAGCCACCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTg                                                                           <  2:126893/251‑1 (MQ=255)
                                                                                                                                                                                gCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACCCACCTCTGGTCATg                                                                     <  2:191870/251‑1 (MQ=255)
                                                                                                                                                                                           ccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCTCATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGAGTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTGCATTGTTGACACACCTCTTGTCATGATAGTATCAat                                                          <  2:97680/251‑1 (MQ=255)
                                                                                                                                                                                            cAGCCCCCCGTCAGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATTTGCATTGATTTGATCGAGTACCTCTTCCGTGGCGCTTATTGAAGGTGTGGAAATCAGAGCGCGGTACTTCTCGATAAATGGCGGTAATTTGGTTTTCATTGTTTCAACACCCCTGGGTATGATAGTATCAat                                                          >  2:179907/1‑250 (MQ=255)
                                                                                                                                                                                                ccccTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTc                                                      >  1:357468/1‑250 (MQ=255)
                                                                                                                                                                                                         gATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATAAATATTCCTGCAGTAtt                                            >  2:65383/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccGAGTTCCTGGAACATGCTACCCTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCTCCAGCAGAGTTATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGTCGCTTATTTAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGTCGGTAATTTTTTTTTCATTGTTGACACACCTCTGTTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACTTTGac                 <  2:167578/250‑2 (MQ=255)
                                                                                                                                                                                                                                              ggAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGtaataa         <  1:162015/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTGTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACccac  <  2:176820/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
GGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCAC  >  NC_000913/4154444‑4154936

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