Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 506,820 T→C 9.3% D437D (GAT→GAC)  trpE → anthranilate synthase component 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947506,8200TC 9.3% 123.0 / 4.1 43D437D (GAT→GAC) trpEanthranilate synthase component 1
Reads supporting (aligned to +/- strand):  ref base T (22/17);  new base C (2/2);  total (24/19)
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.81e-01

GCGATTCTGCCGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGGAAGAAACCATCAACAAGCGCCGGGCGATGTTCCG  >  NC_002947/506686‑506954
                                                                                                                                      |                                                                                                                                      
gCGATTCTGCCGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATa                                                                                                                                       >  1:607460/1‑136 (MQ=255)
      cTGCCGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCaa                                                                                                                                 >  1:513461/1‑136 (MQ=255)
        gCCGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGACACCGCAAtt                                                                                                                               >  2:557095/1‑136 (MQ=255)
        gCCGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGACACCGCAAtt                                                                                                                               >  2:567601/1‑136 (MQ=255)
          cGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGc                                                                                                                             <  2:485300/136‑1 (MQ=255)
          cGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGc                                                                                                                             <  2:668737/136‑1 (MQ=255)
                    cTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTAcc                                                                                                                   <  2:437044/136‑1 (MQ=255)
                    cTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTAcc                                                                                                                   <  2:119178/136‑1 (MQ=255)
                          gggTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTg                                                                                                             >  1:543120/1‑136 (MQ=255)
                                gcCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATTGATACCGCAATTGCCATTCGTACCGCAGTGATCAAg                                                                                                       >  1:14535/1‑136 (MQ=255)
                                gcCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAg                                                                                                       >  2:460005/1‑136 (MQ=255)
                                 cccAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAgg                                                                                                      >  1:510038/1‑136 (MQ=255)
                                  ccAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGa                                                                                                     >  2:339568/1‑136 (MQ=255)
                                  ccAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGa                                                                                                     <  1:467007/136‑1 (MQ=255)
                                    aaGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACg                                                                                                   <  1:327266/136‑1 (MQ=255)
                                          cGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGc                                                                                             <  1:579706/136‑1 (MQ=255)
                                          cGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGACACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGc                                                                                             <  1:557095/136‑1 (MQ=255)
                                          cGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGACACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGc                                                                                             <  1:567601/136‑1 (MQ=255)
                                               gATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCAt                                                                                        >  1:432055/1‑136 (MQ=255)
                                                 tGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCAtgt                                                                                      >  1:517382/1‑136 (MQ=255)
                                                 tGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCAtgt                                                                                      >  1:438020/1‑136 (MQ=255)
                                                   gagaTCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGc                                                                                    <  1:138483/136‑1 (MQ=255)
                                                    agaTCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGcc                                                                                   >  2:643337/1‑135 (MQ=255)
                                                    agaTCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCa                                                                                   >  2:341444/1‑136 (MQ=255)
                                                           cgacgaGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGt                                                                            <  2:513461/136‑1 (MQ=255)
                                                                              aaGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTggcg                                                                        >  1:44479/1‑121 (MQ=255)
                                                                              aaGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTggcg                                                                        <  2:44479/121‑1 (MQ=255)
                                                                              aaGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGAca                                                         >  2:266152/1‑135 (MQ=255)
                                                                              aaGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACt                                                         >  1:366502/1‑136 (MQ=255)
                                                                                 cGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATTGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCgg                                                      <  2:14535/136‑1 (MQ=255)
                                                                                 cGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCgg                                                      <  1:554971/136‑1 (MQ=255)
                                                                                             ggcggTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTgg                                          >  1:653515/1‑136 (MQ=255)
                                                                                                   gCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGc                                                                                                                             >  1:616672/1‑47 (MQ=255)
                                                                                                   gCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGc                                                                                                                             <  2:616672/47‑1 (MQ=255)
                                                                                                   gCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGCCTCGGTGCCGGCGCTGGAATggg                                    >  1:435883/1‑136 (MQ=255)
                                                                                                    ccGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGga                                   >  1:240231/1‑136 (MQ=255)
                                                                                                     cGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGgaa                                  >  1:587786/1‑136 (MQ=255)
                                                                                                        gggCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGgaagaa                               <  2:671466/135‑1 (MQ=255)
                                                                                                                     tGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGGAAGAAACCATCAACAAgc                  <  2:110888/136‑1 (MQ=255)
                                                                                                                     tGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGGAAGAAACCATCAACAAgc                  <  1:462311/136‑1 (MQ=255)
                                                                                                                                cATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGGAAGAAACCATCAACAAGCGCCGGGCGATg       >  1:571249/1‑136 (MQ=255)
                                                                                                                                 aTGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGGAAGAAACCATCAACAAGCGCCGGGCGATGt      <  1:565595/136‑1 (MQ=255)
                                                                                                                                     aTACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGGAAGAAACAATCAACAAGCGCCGGGCGATGTTCCg  >  2:558192/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GCGATTCTGCCGGCTGGCACCTTGTCGGGTGCGCCCAAGATCCGTGCGATGGAGATCATCGACGAGCTGGAACCGGTCAAGCGTGGGGTTTACGGCGGTGCCGTCGGCTATTTTGCCTGGAACGGCAACATGGATACCGCAATTGCCATTCGTACCGCAGTGATCAAGGACGGCGAGCTGCATGTGCAGGCCGGTGGCGGCATCGTTGCCGACTCGGTGCCGGCGCTGGAATGGGAAGAAACCATCAACAAGCGCCGGGCGATGTTCCG  >  NC_002947/506686‑506954

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