Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,401,128 A→G 7.3% N10N (AAT→AAC)  PP_3879 ← HK97 family portal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,401,1280AG 7.3% 152.8 / 2.8 55N10N (AAT→AAC) PP_3879HK97 family portal protein
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (2/2);  total (28/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.02e-01

ACGCTCTCCGCCAGCAGCGTGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCAGAAATGAGCGCAGCGGGCACATGGAT  >  NC_002947/4400982‑4401273
                                                                                                                                                   |                                                                                                                                                 
aCGCTCTCCGCCAGCAGCGTGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTCCCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCAtt                                                                                                                                                 >  6:302505/1‑149 (MQ=255)
 cGCTCTCCGCCAGCAGCGTGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCAt                                                                                                                                                  <  7:34905/147‑1 (MQ=255)
      tcCGCCAGCAGCGTGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTgccgc                                                                                                                                            >  4:173685/1‑148 (MQ=255)
             gcagcGTGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCaag                                                                                                                                        >  6:327153/1‑145 (MQ=255)
                gcGTGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTg                                                                                                                                 <  6:95979/149‑1 (MQ=255)
                   tGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCt                                                                                                                               >  3:206817/1‑148 (MQ=255)
                       acacGTCTGCAATACAGTAATTGCCAACGCTCGGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCAGCC‑AAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTgcc                                                                                                                                              <  8:70142/128‑1 (MQ=255)
                                cAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAAt                                                                                                                   <  3:173685/147‑1 (MQ=255)
                                   gACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGcc                                                                                                              >  1:188471/1‑149 (MQ=255)
                                    aaacttaTTTCCAACCCCCTTTCCGGGGTGACCCGTTA‑TCCCGCTGCTTTCCGGCTATTCCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGcct                                                                                                             <  7:52564/143‑1 (MQ=255)
                                    aCAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGcct                                                                                                             >  1:4518/1‑149 (MQ=255)
                                    aCAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGcct                                                                                                             >  5:272337/1‑149 (MQ=255)
                                      aGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGcctcc                                                                                                           >  5:356506/1‑149 (MQ=255)
                                        tAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTgaagaa                                                                                                                         <  4:89865/133‑1 (MQ=255)
                                        tAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTgaagaa                                                                                                                         >  3:89865/1‑133 (MQ=255)
                                           ttGCCAACGCGCTGTCTGGTGTGA‑CCGGTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTctg                                                                                                      <  6:272337/149‑1 (MQ=255)
                                                     gcTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGcctcc                                                                                                           >  4:95965/1‑134 (MQ=255)
                                                     gcTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGcctcc                                                                                                           <  3:95965/134‑1 (MQ=255)
                                                     gcTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGttt                                                                                            <  4:206817/149‑1 (MQ=255)
                                                     gcTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGttt                                                                                            <  1:229782/149‑1 (MQ=255)
                                                         tCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTgaga                                                                                        <  2:4518/149‑1 (MQ=255)
                                                         tCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTgaga                                                                                        <  5:327153/149‑1 (MQ=255)
                                                             gtgtGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTgaga                                                                                        >  6:135077/1‑145 (MQ=255)
                                                             gtgtGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTgaga                                                                                        <  5:135077/145‑1 (MQ=255)
                                                                  a‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCgg                                                                               >  8:169598/1‑149 (MQ=255)
                                                                      gTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCGTTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTct                                                                                                       <  3:241127/122‑1 (MQ=255)
                                                                      gTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCGTTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTct                                                                                                       <  3:245782/122‑1 (MQ=255)
                                                                      gTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCGTTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTct                                                                                                       >  4:241127/1‑122 (MQ=255)
                                                                      gTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCGTTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTct                                                                                                       >  4:245782/1‑122 (MQ=255)
                                                                       ttACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTgag                                                                                         <  5:72321/135‑1 (MQ=255)
                                                                       ttACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTgag                                                                                         >  6:72321/1‑135 (MQ=255)
                                                                              ggCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCttt                                                                    <  2:188471/149‑1 (MQ=255)
                                                                              ggCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCttt                                                                    <  6:356506/149‑1 (MQ=255)
                                                                                    tgAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTctgggctgg                                                                                                >  1:125278/1‑115 (MQ=255)
                                                                                    tgAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTctgggctgg                                                                                                <  2:125278/115‑1 (MQ=255)
                                                                                     gAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCttt                                                                    <  5:339391/142‑1 (MQ=255)
                                                                                     gAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCttt                                                                    >  6:339391/1‑142 (MQ=255)
                                                                                        cGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGgccagc                                                          >  2:113941/1‑149 (MQ=255)
                                                                                        cGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGgccagc                                                          >  2:160405/1‑149 (MQ=255)
                                                                                            tAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAg                                                      >  8:26461/1‑149 (MQ=255)
                                                                                                     ccccGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCttt                                                                    <  2:246498/126‑1 (MQ=255)
                                                                                                     ccccGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCttt                                                                    >  1:246498/1‑126 (MQ=255)
                                                                                                       tctcgcagcCAAATACTTCTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGAt                                            <  1:127280/145‑1 (MQ=255)
                                                                                                         gcgcCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGCCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAg                                         <  5:302505/149‑1 (MQ=255)
                                                                                                           gccgccAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGac                                            >  2:127280/1‑143 (MQ=255)
                                                                                                                      ctactGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCa                            <  7:169598/149‑1 (MQ=255)
                                                                                                                                  tCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCAGAAATGAgcg                  >  5:126643/1‑147 (MQ=255)
                                                                                                                                     gAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCAGAAATGAGCGCAGc              <  5:6648/148‑1 (MQ=255)
                                                                                                                                     gAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCAGAAATGAGCGCAGc              >  6:6648/1‑148 (MQ=255)
                                                                                                                                           aGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGcc                                      <  6:128369/118‑1 (MQ=255)
                                                                                                                                           aGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGcc                                      >  5:128369/1‑118 (MQ=255)
                                                                                                                                               cTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAg                                         >  6:171322/1‑111 (MQ=255)
                                                                                                                                               cTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAg                                         <  5:171322/111‑1 (MQ=255)
                                                                                                                                                ttCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCAGAAATGAGCGCAGCGGGCACATGGAt  >  5:236219/1‑149 (MQ=255)
                                                                                                                                                  cATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCAGAAACGAGCGCAGCGGGCACATGGAt  >  2:93532/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACGCTCTCCGCCAGCAGCGTGACACACGTCTGCAAGACAGTAATTGCCAACGCGCTGTCTGGTGTGA‑CCGTTACTCCGGCTGCTGAGCGGCTAGACCCTACCCCGCGCCGCCAAAAACTACTGCCGCCATCCGAAACCAGCCCTTCATTGCCGCCAAGAAGATTGCTGAAGAACATAATCAGCCTCCTTCTGGGCTGGGTTTGAGAGGTACTCGGGCTGCGGCTTTGTCGGCCAGCCAAGCCCAGCCGATCAGGCCCAATCCGGCAGAAATGAGCGCAGCGGGCACATGGAT  >  NC_002947/4400982‑4401273

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