Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 510,145 A→T 12.1% E54D (GAA→GAT trpD → anthranilate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947510,1450AT12.1% 90.7 / 3.0 33E54D (GAA→GATtrpDanthranilate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (13/16);  new base T (2/2);  total (15/18)
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.92e-01

CCGTATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGCTTCGTCCTTCGTCCTGGCG  >  NC_002947/510007‑510289
                                                                                                                                          |                                                                                                                                                
ccGTATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGCGGGCCCCAAACGGGCCTTCCTTAGTGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTcggtgc                                                                                                                                       >  1:228654/1‑150 (MQ=255)
   tATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgcttt                                                                                                                                    >  1:365659/1‑147 (MQ=255)
    aTCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgctttc                                                                                                                                   <  2:365659/150‑5 (MQ=255)
        tcgGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgcttt                                                                                                                                    <  1:202042/145‑4 (MQ=255)
        tcgGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgcttt                                                                                                                                    >  2:202042/1‑142 (MQ=255)
          gggCAGTTGGACCTGACTAGAGAAGAAATGACAGAGGTCATGCTCCAGATCATGACCTGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTg                                                                                                                             >  2:315866/1‑150 (MQ=255)
          gggCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTg                                                                                                                             <  2:179288/150‑1 (MQ=255)
            gCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGAACCAGCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGAt                                                                                                                           <  1:469364/150‑1 (MQ=255)
              aGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGc                                                                                                                         <  1:2455/150‑1 (MQ=255)
                      cTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                                       <  2:512958/144‑1 (MQ=255)
                      cTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                                       >  1:512958/1‑144 (MQ=255)
                           tACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGa                                                                                                            >  2:488989/1‑150 (MQ=255)
                             ccGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGaaa                                                                                                          >  2:176117/1‑150 (MQ=255)
                               gaagaaATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGaaaaa                                                                                                        >  1:71844/1‑149 (MQ=255)
                                aagaaATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAgg                                                                                                       <  1:488989/150‑1 (MQ=255)
                                         gcgAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGc                                                                                              <  1:455975/150‑1 (MQ=255)
                                                 aTGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAg                                                                                                                    >  1:424297/1‑120 (MQ=255)
                                                 aTGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAg                                                                                                                    <  2:424297/120‑1 (MQ=255)
                                                               gACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                                       >  1:372415/1‑103 (MQ=255)
                                                               gACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                                       <  2:372415/103‑1 (MQ=255)
                                                                aCCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGat                                                                       <  1:9998/150‑1 (MQ=255)
                                                                       aGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCg                                                                <  2:228654/150‑1 (MQ=255)
                                                                        gTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGAATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCgg                                                               <  1:315866/150‑1 (MQ=255)
                                                                                     ctgATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGAt                                                  >  2:471632/3‑150 (MQ=255)
                                                                                       gATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATgg                                                <  1:471632/150‑1 (MQ=255)
                                                                                           ggCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGcc                                            <  2:267731/150‑1 (MQ=255)
                                                                                                       aTGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGACGAGCTGCAAAACCTCGATCGTGTCGTCGATATTGTATGGACCGGTCGTGATGTCGCCTAAATCTTCAAc                                >  2:61793/1‑150 (MQ=255)
                                                                                                                            gAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGcttcgtccttc           <  1:176117/150‑1 (MQ=255)
                                                                                                                                 cATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGcttcgtccttc           <  2:592367/145‑1 (MQ=255)
                                                                                                                                 cATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGcttcgtccttc           >  1:592367/1‑145 (MQ=255)
                                                                                                                                     gacgaAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCaa                                 <  2:216165/119‑1 (MQ=255)
                                                                                                                                     gacgaAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCaa                                 >  1:216165/1‑119 (MQ=255)
                                                                                                                                     gacgaAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGCTTCGTCCTTCGTCCTggcg  >  1:371404/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                
CCGTATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGCTTCGTCCTTCGTCCTGGCG  >  NC_002947/510007‑510289

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