Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,322,759 G→A A63V (GCT→GTT)  trpE ← anthranilate synthase subunit TrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,322,7590GA93.1% 88.8 / ‑6.7 29A63V (GCT→GTT) trpEanthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/16);  minor base C (1/1);  total (12/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 = 1.00e+00

TCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAGCTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCT  >  NC_000913/1322522‑1322995
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
tctcGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATgcgca                                                                                                                                                                                                                                  <  1:461878/250‑1 (MQ=255)
          ccttcGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCt                                                                                                                                                                                                                       <  2:373936/251‑1 (MQ=255)
                gTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGCGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGCCAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACcagc                                                                                                                                                                                                                 >  1:62665/1‑251 (MQ=255)
                  aCATTCAACAGATTCTGCAATAAACGGAAAGCGTCACAAACCGCAAGGGAGCATAAGCGGGCGTCTTCCTCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCTTTATCCAGTATTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGTTTGTGACAGGGTCACCTAAACCTGTAAGTCGCAGCGCACTTTCTACcagca                                                                                                                                                                                                                >  2:107834/1‑250 (MQ=255)
                   cATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACcagcagc                                                                                                                                                                                                              <  2:442414/251‑1 (MQ=255)
                       cAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGc                                                                                                                                                                                                          <  1:248920/251‑1 (MQ=255)
                                  gCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAatca                                                                                                                                                                                               <  1:261188/251‑1 (MQ=255)
                                  gCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAat                                                                                                                                                                                                 <  1:373687/249‑1 (MQ=255)
                                     aaaaaCGGAAAGCTTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGTAGTGGACTTACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCTAGGAGGGCTTCGCCTGTGCCGGAAAGGGCCTTGATTGCGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGCCTACCAGCAGCAGGCTTTTTAAATCATCt                                                                                                                                                                                            <  2:195684/249‑1 (MQ=255)
                                            gAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGt                                                                                                                                                                                     >  2:218690/1‑251 (MQ=255)
                                               aGCGACAAAAACCGAAAGGGAGCATAAGAGGGCGTCTTCATCCAGCAGGGGACTGACAGGGGGGAAGCGCAGCAAACGGCAGTTTGGTGATTGTTAACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGACGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGAGCAGCGCACTGTCTACCTGCAGCAGGCTTTTTAAATCATCTTTGCTGTCg                                                                                                                                                                                   >  2:393574/1‑250 (MQ=255)
                                                gCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGat                                                                                                                                                                                 >  1:447297/1‑251 (MQ=255)
                                                                                                     gACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGaa                                                                                                                            <  1:441549/251‑1 (MQ=255)
                                                                                                                 gcgcagcgcaCGGCAGTTTTGTGATGGTTCCCCTGTCACACCCGGAGGCAGGGTGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGGGACAGTGTCACCTAAACCTGTAATGCGCAGCGCACTGGCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCCCACAACTGGTGAAAAAGCGCGGTga                                                                                                                <  2:88158/251‑2 (MQ=255)
                                                                                                                           cGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTcgcg                                                                                                      <  1:107834/251‑1 (MQ=255)
                                                                                                                                                  ggccACACCCGCAGGCAGAGCGTTATCCAGTAGGGCCAGGAGGGCTTCGCCCTGGCCGGAAAGTGCCTGGATGGTGACAGTGTCTCCGAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTGGAAATCATCTTTGCTGTCGATACCTGCGGATTCCAGCAGCAGCGTTGCCGGACGAGCCCCACCCAACTGGTGAAAAAGCGCGGTGGGACTGTCGCGATAAGCTCCTTCGCAGGTTAg                                                                                 <  2:447297/247‑1 (MQ=255)
                                                                                                                                                          ccGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGTGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCATCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCgaga                                                                       <  1:28418/251‑1 (MQ=255)
                                                                                                                                                             caggcaggGCGTTATAAAGTCGTGCCAGGAGGGCTTCGCCGTTGCCGGCAAGTGCCTGGGTTGTGACAGGGTCACTTAAAACTGTAATGCGCAGCGCACTGCCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGCTATCTGCGGATTCCAGCAGCAGCGCTGCCGGACGATCCCCACCCAACTGGTGAAAAAGCGCGGTGGGATTGCCGAGGTAAGCGCCTTCGCAGGTTAGCCGTTCGAGAGTc                                                                    >  2:225603/1‑251 (MQ=255)
                                                                                                                                                                   gggCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCGCTGTCTGCCAGCAGCAGGCTTTTTAAAGCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGtttt                                                              <  1:22381/251‑1 (MQ=255)
                                                                                                                                                                     gCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTtg                                                            <  1:218690/251‑1 (MQ=255)
                                                                                                                                                                       gTTATCCAGTAGGGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTtgtg                                                          <  1:104617/251‑1 (MQ=255)
                                                                                                                                                                         tATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTCGCCGTTCGGGAGTCGGTTTTTGGGtt                                                        >  2:11551/1‑251 (MQ=255)
                                                                                                                                                                              agtagtGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCAt                                                   >  1:423161/1‑251 (MQ=255)
                                                                                                                                                                                  gtGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGtt                                               >  1:13393/1‑251 (MQ=255)
                                                                                                                                                                                                  gCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGtt                               <  1:389889/251‑1 (MQ=255)
                                                                                                                                                                                                   ccGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTTGAATGCTCAGCGCACTGTCTTTCAGCATCATGCTTTTTAAATCATCTTTGCTGTCGATATCGGCGGATTCCAGCAGCTGCGTTTCCGGATGATCTTCACACATCTGGTGATAAAGCGTGGTTGGATTGTCGTGATGAGCGCCTTATCAGGTTAGTAGTTCGTGAGGCGGTTTTTGTGGTTGAATTGTTATTCTCTAAGTTTGTTc                              <  2:304389/251‑1 (MQ=255)
                                                                                                                                                                                                            gAAAGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCaaaaaaaa                      >  1:402260/1‑250 (MQ=255)
                                                                                                                                                                                                               aGTGCCTGGATTGTGACAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGccc                  >  2:110979/1‑251 (MQ=255)
                                                                                                                                                                                                                               cAGTGTCACCTAAAACTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCt  >  2:126832/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
TCTCGTTCTTCCTTCGGTACATTCAACAGATTCTGCAATAAACGGAAAGCGTCAAAAACCGAAAGGGAGCATAAGCGGGCGTCTTCATCCAGCAGTGGACTGACAGGGGGGAAGCGCAGCACACGGCAGTTTGGTGATTGTTCACTTTCCACACCCGCAGGCAGGGCGTTATCCAGTAGTGCCAGGAGGGCTTCGCCGTTGCCGGAAAGTGCCTGGATTGTGACAGTGTCACCTAAAGCTGTAATGCGCAGCGCACTGTCTACCAGCAGCAGGCTTTTTAAATCATCTTTGCTGTCGATATCTGCGGATTCCAGCAGCAGCGTTGCCGGACGATCCCCACACAACTGGTGAAAAAGCGCGGTGGGATTGTCGCGATAAGCGCCTTCGCAGGTTAGCAGTTCGAGAGTCGGTTTTTGTGTTTGCATTGTTATTCTCTAATTTTGTTCAAAAAAAAGCCCGCTCATTAGGCGGGCT  >  NC_000913/1322522‑1322995

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