Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 519,340 C→G 7.4% I242M (ATC→ATG)  argC → N‑acetyl‑gamma‑glutamyl‑phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947519,3400CG 7.4% 165.5 / 3.5 54I242M (ATC→ATG) argCN‑acetyl‑gamma‑glutamyl‑phosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base G (2/2);  total (27/27)
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.29e-01

TTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTG  >  NC_002947/519194‑519486
                                                                                                                                                  |                                                                                                                                                  
ttCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAGGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGCCGCCGATGATCCGTGGCATCCAt                                                                                                                                                 <  2:127109/150‑1 (MQ=255)
    gCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCa                                                                                                                                             <  2:1191953/150‑1 (MQ=255)
            ggCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcgga                                                                                                                                         <  1:136245/146‑4 (MQ=255)
            ggCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcgga                                                                                                                                         >  2:136245/1‑143 (MQ=255)
                     aTGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTg                                                                                                                            <  2:1132769/150‑1 (MQ=255)
                       gAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGt                                                                                                                          >  1:480889/1‑150 (MQ=255)
                        aaaGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTc                                                                                                                         <  2:480889/150‑1 (MQ=255)
                        aaaGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTc                                                                                                                         <  1:67090/150‑1 (MQ=255)
                             cTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACCTCGGCCTGCCTTTCGTGCCGACTCCCCCGCAGCTGATCGGTGGCATCCGTGCCGCCCTGGACTCGAACGTGGTCGATAc                                                                                                                    <  2:374971/150‑1 (MQ=255)
                                          ggCCATCTCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTg                                                                                                       >  1:1194790/1‑150 (MQ=255)
                                              aTCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAgg                                                                                                   <  1:1297808/150‑1 (MQ=255)
                                                     cTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGtt                                                                                            <  2:1212850/150‑1 (MQ=255)
                                                           ggAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATAc                                                                                                                    <  2:225216/120‑1 (MQ=255)
                                                           ggAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATAc                                                                                                                    >  1:225216/1‑120 (MQ=255)
                                                           ggAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGaa                                                                                      >  1:854637/1‑150 (MQ=255)
                                                                tCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCt                                                                                 <  1:417426/150‑1 (MQ=255)
                                                                    ccAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGc                                                                             >  1:364960/1‑150 (MQ=255)
                                                                           cTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGa                                                                          <  1:232446/146‑1 (MQ=255)
                                                                           cTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGa                                                                          >  2:232446/1‑146 (MQ=255)
                                                                           cTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGaa                                                                      <  2:1126566/150‑1 (MQ=255)
                                                                                  tGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTa                                                                                                                                      <  1:808320/79‑1 (MQ=255)
                                                                                  tGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTa                                                                                                                                      >  2:808320/1‑79 (MQ=255)
                                                                                    gCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTtcgtc                                                             >  2:169466/1‑150 (MQ=255)
                                                                                       gcgGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACg                                                                                   <  1:73914/125‑1 (MQ=255)
                                                                                       gcgGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACg                                                                                   >  2:73914/1‑125 (MQ=255)
                                                                                       gcgGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGAc                                                          <  1:268361/150‑1 (MQ=255)
                                                                                             aaGGACATCGGCCTGACTTTCGTGCCACACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATg                                                    >  1:710511/1‑150 (MQ=255)
                                                                                                 aCATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATg                                                                                                                                                <  1:1038493/54‑1 (MQ=255)
                                                                                                 aCATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATg                                                                                                                                                >  2:1038493/1‑54 (MQ=255)
                                                                                                    tCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCcagcag                                             >  1:159292/1‑150 (MQ=255)
                                                                                                       gCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggca                                          >  1:354233/1‑150 (MQ=255)
                                                                                                        ccTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                         <  2:991010/150‑1 (MQ=255)
                                                                                                        ccTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                         <  1:850686/150‑1 (MQ=255)
                                                                                                         cTGACTTTCGTGCCGCAACTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGc                                        >  1:1047627/1‑150 (MQ=255)
                                                                                                          tGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTCGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCCACGCCGAGCAACCGATCGTCGACGGGATGCCAGCCGGCAGcc                                       >  1:288984/1‑150 (MQ=255)
                                                                                                              tttCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATcc                                   >  2:1253211/1‑150 (MQ=255)
                                                                                                              tttCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATcc                                   >  1:1068090/1‑150 (MQ=255)
                                                                                                                tCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCgg                                 <  2:159292/150‑1 (MQ=255)
                                                                                                                    gCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAAc                             >  1:1158009/1‑150 (MQ=255)
                                                                                                                         aCCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTa                        <  2:354233/150‑1 (MQ=255)
                                                                                                                           cTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGAAAACCCGTAgc                      >  2:454024/1‑150 (MQ=255)
                                                                                                                             gACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGGTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGAAAACCCGTAgcgt                    <  1:454024/150‑1 (MQ=255)
                                                                                                                              aCGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGAGCTGTTCGAGAAACGCTACGCCGAt                                                                         <  1:1324539/96‑1 (MQ=255)
                                                                                                                              aCGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGAGCTGTTCGAGAAACGCTACGCCGAt                                                                         >  2:1324539/1‑96 (MQ=255)
                                                                                                                                   gatgatCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                         >  1:1049046/1‑123 (MQ=255)
                                                                                                                                   gatgatCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                         <  2:1049046/123‑1 (MQ=255)
                                                                                                                                   gatgatCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCgtgg              <  2:854637/150‑1 (MQ=255)
                                                                                                                                     tgatCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCgtggtg            >  2:802568/1‑150 (MQ=255)
                                                                                                                                        tCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCa         >  1:82940/1‑150 (MQ=255)
                                                                                                                                          ggTGGCATGGATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGc                                                                             >  1:1371072/2‑80 (MQ=255)
                                                                                                                                          ggTGGCATGGATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGc                                                                             <  2:1371072/79‑1 (MQ=255)
                                                                                                                                          cGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAAc       >  2:421478/1‑150 (MQ=255)
                                                                                                                                            tGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGt     >  1:1161235/1‑150 (MQ=255)
                                                                                                                                               cATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTg  <  2:1161235/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTG  >  NC_002947/519194‑519486

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