Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,074,945 T→G 18.1% N1313T (AAC→ACC)  putA ← fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,074,9450TG18.1% 32.6 / 2.5 20N1313T (AAC→ACC) putAfused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base T (10/6);  new base G (2/2);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCG  >  NC_000913/1074803‑1075087
                                                                                                                                              |                                                                                                                                              
ttGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGCGGTTGTCACATCCTCCAGCTACCTGTTTAAACTATAGTCATTAAGCTGGCGt                                                                                                                                                >  1:311052/1‑143 (MQ=255)
         gcggTATGTGTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcag                                                                                                                                       <  2:15784/141‑1 (MQ=255)
         gAGGTATTTTTGGATCGGTCGGTTGGAATCGTGCGGCGGTTTTTTATTGAGGGGGTTCGCCCACGTTAAAGATGCCGGTGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcag                                                                                                                                       <  2:189377/143‑1 (MQ=255)
          gggTATTTTTCAAACAAACAGAAAATAACGTGCCGGGGGTTTTCCCTTACTGCATTCGACCACGTTAAAGAATCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagc                                                                                                                                      <  1:389577/142‑1 (MQ=255)
             gaTTGTTCATTCAATGGGTTGAAAACGTGCAGCGGATTATTAGGGAGGGGATTCGACCACGTTAAAGGTGGCGGAGGAGGTTGTAACGTCCTCCGGCTACCTGTTTAACCTATAGTCATTCAGCTGGCGGTAACGCcagctta                                                                                                                                   <  2:90559/142‑4 (MQ=255)
                       tcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCAGGAGGAGGTTGTAACGCGCTCCGGCTACCTGTTTAACCTATAGTCATTCAGCTGGCGGTACCGCCAGGGGCAGCGGTATTc                                                                                                                         >  2:62582/1‑143 (MQ=255)
                        caaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGGTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTca                                                                                                                        >  2:17492/1‑143 (MQ=255)
                               ttGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAg                                                                                                                 <  1:218844/143‑1 (MQ=255)
                                          cAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCAATAGTCATTAAACTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGa                                                                                                      >  1:468155/1‑143 (MQ=255)
                                                   attattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTCTAGTCATTAAGCTGGCGTTACCGCCAGCGGCCGCGGTGTTCACACTCAGCGAACGCTCGATATACAGcc                                                                                             >  2:313758/1‑143 (MQ=255)
                                                            gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCGTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGCCCGATCGGGGCACGGCCGCTCCAGaa                                                                                    >  1:119404/1‑143 (MQ=255)
                                                            gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                                                                                    <  1:240893/143‑1 (MQ=255)
                                                                                    cGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTACGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACACCCGTTCCCGCAGGATATTGCTTTCGCCACGGGCaa                                                            >  2:98183/1‑143 (MQ=255)
                                                                                    cGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGCACTCACACTCAGCGAACGCTCGACCTACCGCCGCTCCCGCAGGATCTTGCATCCGCCACGGGcca                                                            >  2:154070/1‑141 (MQ=255)
                                                                                                        tccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACCCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGCTATTGCTTTCGCCACGCGCAACACCCTGCCCCGCAACAAtt                                        >  2:484633/1‑143 (MQ=255)
                                                                                                         ccGGCTACCTGTTTCACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtg                                       <  1:398837/143‑1 (MQ=255)
                                                                                                          cGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgt                                      <  1:515738/143‑1 (MQ=255)
                                                                                                                    tttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCGGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCcgc                            >  2:70519/1‑143 (MQ=255)
                                                                                                                                  aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTTTTCACACTCAGCGAACGCACGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTc              >  2:281440/1‑143 (MQ=255)
                                                                                                                                  aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTCCGCCACAGGCAAGCCCCTGCAACGAACCTAGTGTTCCATCCCGCGCAGCAACttcctc              >  2:494643/1‑138 (MQ=255)
                                                                                                                                        cTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACaa        >  1:515525/1‑143 (MQ=255)
                                                                                                                                              ttACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCACCTGCTTCACACAAGgcgcc  >  2:419942/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                              
TTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCG  >  NC_000913/1074803‑1075087

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