Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,056,972 G→T 12.5% P159Q (CCG→CAG)  PP_3571 ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,056,9720GT12.5% 94.5 / 2.5 32P159Q (CCG→CAG) PP_3571acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (14/14);  new base T (2/2);  total (16/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 = 8.00e-01

TCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCCAACACCGAGGCGA  >  NC_002947/4056838‑4057107
                                                                                                                                      |                                                                                                                                       
tCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACgg                                                                                                                                        <  1:127238/136‑1 (MQ=255)
tCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACgg                                                                                                                                        <  2:85450/136‑1 (MQ=255)
tCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACgg                                                                                                                                        <  2:30006/136‑1 (MQ=255)
  aCCCCATACGGCGCGTACGCCGAATGGCCCTCTGCGCCCTGCACCTCGCAACGCAATGCCAGCTTGCCCTAGTGCCCCAGTACAGGTTTGAGTTCGGTGGGTTAACCGATCAGGCACAGCCGAGGTTTACACGGgc                                                                                                                                      >  1:497152/1‑136 (MQ=255)
          cGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTgg                                                                                                                              <  1:472582/136‑1 (MQ=255)
          cGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTGTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTgg                                                                                                                              <  1:372782/136‑1 (MQ=255)
            gcgcGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAg                                                                                                                            <  1:62237/136‑1 (MQ=255)
               cgTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTTGGTTCACCGATCAGGCACAGCCGTGGTTTACACTGGCGCTGTTGGAGGGc                                                                                                                         <  2:270321/136‑1 (MQ=255)
                gTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcc                                                                                                                        >  1:32881/1‑135 (MQ=255)
                gTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcc                                                                                                                        >  2:304295/1‑135 (MQ=255)
                 tACGCCGAATTGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACTGGCGCTGTTGGAg                                                                                                                            <  1:443456/131‑1 (MQ=255)
                 tACGCCGAATTGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACTGGCGCTGTTGGAg                                                                                                                            >  2:443456/1‑131 (MQ=255)
                 tACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcag                                                                                                                       >  1:412799/1‑136 (MQ=255)
                                          gCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAaaca                                                                                              <  2:322082/136‑1 (MQ=255)
                                              cTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCCCGCCCAAACAACCg                                                                                          >  1:445697/1‑136 (MQ=255)
                                                   aaCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAgcg                                                                                                          >  2:489411/1‑115 (MQ=255)
                                                   aaCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAgcg                                                                                                          <  1:489411/115‑1 (MQ=255)
                                                   aaCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGActtc                                                                                     >  1:141431/1‑136 (MQ=255)
                                                   aaCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGActtc                                                                                     >  2:472750/1‑136 (MQ=255)
                                                             cAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAg                                                                           >  2:372110/1‑136 (MQ=255)
                                                             cAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAg                                                                           >  1:470321/1‑136 (MQ=255)
                                                                  tGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGTGAAGCCAGCATCGAGCGCACGCCCAAACAACCTACTTCTTCGTCATAGGaaaa                                                                      <  2:177042/136‑1 (MQ=255)
                                                                     ccTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGc                                                                   <  1:372110/136‑1 (MQ=255)
                                                                            ccccAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACTGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGgtgt                                                            >  1:100061/1‑136 (MQ=255)
                                                                                cgTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCg                                                        <  1:173998/135‑1 (MQ=255)
                                                                                  tACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGc                                                      <  1:46144/136‑1 (MQ=255)
                                                                                          tGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCaa                                              >  2:49574/1‑136 (MQ=255)
                                                                                            aGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACg                                            >  1:476609/1‑136 (MQ=255)
                                                                                               tCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCt                                         >  2:190416/1‑136 (MQ=255)
                                                                                                         cACCGATCAGGAACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCATGTGTACCGGCATGCGCAACGGCTGCGCGAGaaa                               <  1:306305/136‑1 (MQ=255)
                                                                                                                    cacaGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCAcag                    >  2:423472/1‑136 (MQ=255)
                                                                                                                              gTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGTAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCCAacac          <  2:454684/136‑1 (MQ=255)
                                                                                                                              gTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCCAacac          <  2:18574/136‑1 (MQ=255)
                                                                                                                              gTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGACAAAACACCCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCATCCAacac          <  2:376318/136‑1 (MQ=255)
                                                                                                                               tttACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCCAACAcc         <  2:470321/136‑1 (MQ=255)
                                                                                                                                      gggCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCCAACACCGAGGCGa  >  1:264426/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGCCAGGTGTACCGGCATGCGCAACGGCTGCGCGAGAAACGCGGGCACAGCAGCCAACACCGAGGCGA  >  NC_002947/4056838‑4057107

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