Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,559,538 C→T 10.9% W337* (TGG→TAG) ‡ hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,559,5380CT10.9% 115.2 / 6.0 46W337* (TGG→TAG) ‡hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (2/3);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

GCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATT  >  NC_002947/5559398‑5559679
                                                                                                                                            |                                                                                                                                             
gCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCTGGATgccgg                                                                                                                                      <  1:782720/150‑2 (MQ=255)
  tGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCa                                                                                                                                    <  2:462908/150‑1 (MQ=255)
    gaTCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGt                                                                                                                                  >  1:442308/2‑150 (MQ=255)
     aTCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGtt                                                                                                                                 >  1:475345/1‑150 (MQ=255)
         aTTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCa                                                                                                                             >  2:413019/1‑150 (MQ=255)
                     aaTCCGCCTCGAGAGCAGCAGTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCg                                                                                                                 >  1:867018/1‑150 (MQ=255)
                     aaTCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAg                                                                                                                                             >  1:923468/1‑122 (MQ=255)
                     aaTCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAg                                                                                                                                             <  2:923468/122‑1 (MQ=255)
                     aaTCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCg                                                                                                                 >  2:603968/1‑150 (MQ=255)
                            cTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCNccc                                                                                                          <  2:700121/150‑1 (MQ=255)
                             tCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCTGGATgccgga                                                                                                                                     <  1:487287/122‑3 (MQ=255)
                             tCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCTGGATgccgga                                                                                                                                     >  2:487287/1‑120 (MQ=255)
                              cGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCa                                                                                                        <  1:1022529/150‑1 (MQ=255)
                                agagCAGCAGTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAgt                                                                                                      <  2:867018/150‑1 (MQ=255)
                                     agaaCTCAACGGCTGGCCAGGCATAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGAc                                                                                                 <  2:341975/147‑1 (MQ=255)
                                        aCTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAg                                                                                              >  1:908856/1‑150 (MQ=255)
                                           cAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGtt                                                                                                                                 >  2:815253/1‑112 (MQ=255)
                                           cAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGtt                                                                                                                                 <  1:815253/112‑1 (MQ=255)
                                                    gCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATc                                                                                  >  2:745733/1‑150 (MQ=255)
                                                           agagCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCCCTGTCCGGCATTCAGATGCCGCCAGGTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGgag                                                                           >  1:131852/1‑150 (MQ=255)
                                                             agCCTGAACCACACGCTGGCCCTCGTTGCGCAGATGGGGGACGGGCTGCCAGAGGGCCGCATCGCCACTGTCCGTCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGaa                                                                         <  1:1096730/150‑1 (MQ=255)
                                                               ccTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAAc                                                                       >  2:933903/1‑150 (MQ=255)
                                                               ccTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAAc                                                                       <  2:135966/150‑1 (MQ=255)
                                                                 tGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCa                                                                     >  2:52834/1‑150 (MQ=255)
                                                                      cacaCGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGc                                                                >  2:381065/1‑150 (MQ=255)
                                                                       acacGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGcc                                                               >  1:773612/1‑150 (MQ=255)
                                                                         acGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCgg                                                             <  1:413019/150‑1 (MQ=255)
                                                                         acGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCgg                                                             <  2:442308/150‑1 (MQ=255)
                                                                                            aaCTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTGGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCg                                          <  2:262758/150‑1 (MQ=255)
                                                                                                   aGACTTGCTGCCAGAGGGCCGCATCGCCACTGCCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTATTCACATTCGTGGGGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGa                                   <  2:156746/150‑1 (MQ=255)
                                                                                                           tgcCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATAtcgtcg                           >  1:528097/1‑150 (MQ=255)
                                                                                                             ccAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTa                         <  1:603968/150‑1 (MQ=255)
                                                                                                                gagGGCCGCATCGCCACTGTCCGGCATCCAGATgccgc                                                                                                                                      >  2:942166/1‑38 (MQ=255)
                                                                                                                gagGGCCGCATCGCCACTGTCCGGCATCCAGATgccgc                                                                                                                                      <  1:942166/38‑1 (MQ=255)
                                                                                                                   ggCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACGAGGGTCGTCAAATCCGTGGAGAAAACCGGAGCGTTCCGCGCCCGTCCAGAGTCGGCGCCGATATCGTCGTAGCGAcc                   >  1:1097447/1‑150 (MQ=255)
                                                                                                                        cATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAgt                                                                                                      >  1:742258/1‑62 (MQ=255)
                                                                                                                        cATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAgt                                                                                                      <  2:742258/62‑1 (MQ=255)
                                                                                                                              cACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTg                                                                              <  1:147872/80‑1 (MQ=255)
                                                                                                                              cACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTg                                                                              >  2:147872/1‑80 (MQ=255)
                                                                                                                                   tCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGt                                                                               >  1:602072/1‑74 (MQ=255)
                                                                                                                                   tCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGt                                                                               <  2:602072/74‑1 (MQ=255)
                                                                                                                                    gcGGCATCTGGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgat       <  2:601372/144‑1 (MQ=255)
                                                                                                                                    gcGGCATCTGGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGcgat       >  1:601372/2‑145 (MQ=255)
                                                                                                                                    ccGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGAtt  >  2:129146/1‑150 (MQ=255)
                                                                                                                                         aTCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAAt                                                                                   >  1:98775/1‑64 (MQ=255)
                                                                                                                                         aTCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAAt                                                                                   <  2:98775/64‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GCTGAATCAATTGCCGATCACAATCCGCCTCGAGAGCAGCACTCAACGGCTGGCCAGGCAGAGCCTGAACCACACGCTGGCCCTCGTTGCGCAACTGGCAGACTTGCTGCCAGAGGGCCGCATCGCCACTGTCCGGCATCCAGATGCCGCCAGTTGGCAATACGACCTCCGCTCGCCCCAGTGTGACCAGGGTCTTCAAATCCGTGGAGAAACCAGTGGCCGGGCGCGCCCGGCCAAAGTCGGCGCCGATATCGTCGTAGCGACCGCCCTGGGCGATCGATT  >  NC_002947/5559398‑5559679

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: ‑