Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,320,907 C→T D160N (GAT→AAT)  trpD ← anthranilate synthase subunit TrpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,320,9070CT89.5% 48.5 / NA 19D160N (GAT→AAT) trpDanthranilate synthase subunit TrpD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/9);  minor base G (0/2);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTAAGCGTCTGCGCCTGATACAGTTTTTCCAGAATCGGTTGCAGCGTGTTGGCTGGCTCTAGTTTCTGCTGCGCCCAGGCCAGCGTTTGTTCCAGCAGGCGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATCGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCGCCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGC  >  NC_000913/1320679‑1321140
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
ccAGTTGTTCCGGCTTCAGCCCGCCACGCCCCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTAAGCGTCTGCGCCTGATACAGTTGTTCCAGAATCGGTTGCAGCGTGCTGGCTGGCTCTAGTTTCTGCTGCGCCCAGGCCAGCGTTTGTTCCATCAGGCGAGCGCCCTTGGTGGTGAGAATGGAGTCCGGATGGGACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCGCCATg                                                                                                                                                                                                                      <  1:262276/251‑1 (MQ=255)
                 aGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTAAGAGTCTGCGCCTGATACAGTTTTTACAGAATCGGATGCAGCGTGTTGGCTGGCTCTAGTTTCTGCTGCGCCCAGGCCAGCGTTTGTTCCAGCAGGCGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGCACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGtt                                                                                                                                                                                                     >  1:109427/1‑251 (MQ=255)
                            caccacCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTAAGCGTCTGCGCCTGATACAGTTTTTCCAGAATCGGTTGCAGCGTGTTGGCTGGCTCTAGTTTCTGCTGCGCCCAGGCCAGCGTTTGTTCCAGCAGGCGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAAc                                                                                                                                                                                          >  1:116801/1‑251 (MQ=255)
                                                                                                                   gnnnnnctCNGGTTTCTGCTGCGCCCAGGCCAGCGTTGGTTCCAGCAGGCGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTaa                                                                                                                                                                                                               <  2:303152/143‑1 (MQ=255)
                                                                                                                                    gctgcGCCCAGGCCAGCGTTGGTTCCAGCAGGCGAGCGCCCTGGGTGGTGAGAATTGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATGGTGACGTACTGCCATCACCATGCCATTGAAATGGGCGTGGATGTTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCATCGGCAGCTGGTTTGTTAATCCAGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCt                                                                                  <  1:449197/251‑1 (MQ=255)
                                                                                                                                        agacgtgggctctgccgttcgcagcagGGGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAATCCGTCCGTAATGTTACTGCCAACCAGCGACTGATAACGCGCCACCGGCAGCGGTTTTTTTAATCCGGCAAACATCGCCTGACCGTCAGGTTCAATGCTGGAGGCTTTAc                                                                              <  1:203156/230‑1 (MQ=255)
                                                                                                                                         ncccAGGCCAGCGTTTGTTCCAGCAGGTGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGc                                                                                                                       >  1:99718/2‑209 (MQ=255)
                                                                                                                                         gcCCAGGCCAGCGTTTGTTCCAGCAGGCGAGCGCCCTGGGTGGGGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGNNTTGTGANNt                                                                                                                                                                                                                                     >  2:5071/1‑99 (MQ=255)
                                                                                                                                            cAGGCCAGCGTTTGTTGCAGCAGGCGAGCGCCCTGGGTGGTGAGAATTGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGc                                                                                                                                                                                                                                <  1:322871/101‑1 (MQ=255)
                                                                                                                                                             nagcagGCGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACCAACGCGATCCGCATTGCGACGTCCTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATACCGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTgg                                                                                       >  1:55989/2‑221 (MQ=255)
                                                                                                                                                             cagcagGCGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGTCGTTGATGGTTAAACCGTCCGGAATGTTACTg                                                                                                                                                                         <  1:425080/139‑1 (MQ=255)
                                                                                                                                                                               ggtggtGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCGCCTGACCGACATAGCCCCCGt                                      <  1:5071/251‑1 (MQ=255)
                                                                                                                                                                                           gnATTNCNNANNGAANTGGAATCCACAAACGCGATCCNNATGGTGACGNACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATNTTACTGCCANNCANCGAGTNANNNCNNNCCNCCGGCAGCGGGTTTGTGAANCCNNCAAACANCGCCTGACCGTCATGTTCAATGCTgg                                                                                       <  2:116801/191‑1 (MQ=255)
                                                                                                                                                                                            naNTNCGGATGGAACTGGAATCCACAAACGCGATNCGCATTGTGACGTACTGCCATCANNATNNNAGNNAAATg                                                                                                                                                                                                           <  2:366389/73‑1 (MQ=255)
                                                                                                                                                                                            gATTCCGAATGGAACTGGAATCCACAAACGCGATCCGCATTCTGACGTAATGCCATCACCATGCCATTAAAATg                                                                                                                                                                                                           >  1:366389/1‑74 (MQ=255)
                                                                                                                                                                                                   gATGGAACTGGAATCCACAAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAg                                                                                     >  1:431611/1‑185 (MQ=255)
                                                                                                                                                                                                            ggAATNCACAANCNNGATCCGCATTGTGACGTACNGCCATCACCATGCCANTAAANTGGGCNTTNa                                                                                                                                                                                                   >  2:196885/1‑66 (MQ=255)
                                                                                                                                                                                                            ggAATCCACAAACGCGATCCGCATTGTGACGGACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCGCCTGACCGACATAGCCCCCGTAAGCTTCg                              <  1:196885/230‑1 (MQ=255)
                                                                                                                                                                                                                    caAACGCGATCCGCATTGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCCGCATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGATGCTTTACCAGTGAGAATTTCGCCCTCCTGACCGACATAGCCCCCGTAAGGTTCGACAATCGCCTGATGTCCGAGGCAAATGc  <  1:381993/250‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
CCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTAAGCGTCTGCGCCTGATACAGTTTTTCCAGAATCGGTTGCAGCGTGTTGGCTGGCTCTAGTTTCTGCTGCGCCCAGGCCAGCGTTTGTTCCAGCAGGCGAGCGCCCTGGGTGGTGAGAATGGATTCCGGATGGAACTGGAATCCACAAACGCGATCCGCATCGTGACGTACTGCCATCACCATGCCATTAAAATGGGCGTTGATGGTTAAACCGGCCGGAATGTTACTGCCAACCAGCGAGTGATAACGCGCCACCGGCAGCGGGTTTGTTAATCCGGCAAACATCGCCTGACCGTCATGTTCAATGCTGGAGGCTTTACCGTGGAGAATTTCGCCCGCCTGACCGACATAGCCCCCGTAAGCTTCGACAATCGCCTGATGTCCGAGGCAAATGC  >  NC_000913/1320679‑1321140

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