Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,357,496 G→C 10.6% V289L (GTG→CTG) ‡ PP_3831 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,357,4960GC10.6% 128.6 / 7.5 47V289L (GTG→CTG) ‡PP_3831hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base C (3/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

ACTGCCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAG  >  NC_002947/4357349‑4357642
                                                                                                                                                   |                                                                                                                                                  
aCTGCCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAg                                                                                                                                                    <  6:224735/148‑1 (MQ=255)
    ccAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGc                                                                                                                                               >  5:63666/1‑149 (MQ=255)
    ccAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGc                                                                                                                                               <  8:189597/149‑1 (MQ=255)
     cAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGcc                                                                                                                                              <  2:216584/149‑1 (MQ=255)
      aaGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGACACGCGCCTGGGCAATACGCCGGCAGTGTGCCg                                                                                                                                             >  8:284213/1‑149 (MQ=255)
        ggACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgc                                                                                                                                           <  3:299580/149‑1 (MQ=255)
              tCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCgctgct                                                                                                                                     <  2:6178/149‑1 (MQ=255)
                                  gcCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTg                                                                                                                 >  3:269739/1‑149 (MQ=255)
                                    cTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTg                                                                                                                 >  2:304891/1‑147 (MQ=255)
                                         aTTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgcg                                                                                                                                          <  5:121715/117‑1 (MQ=255)
                                         aTTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgcg                                                                                                                                          >  6:121715/1‑117 (MQ=255)
                                            gAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTAt                                                                                                        <  2:415/148‑1 (MQ=255)
                                               ccTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCt                                                                                                    >  1:141767/1‑149 (MQ=255)
                                                  gCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTCGAGCACTATGCCTAgg                                                                                                 >  2:101681/1‑149 (MQ=255)
                                                    tAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTgggg                                                                                               >  2:107282/1‑149 (MQ=255)
                                                    tAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTgggg                                                                                               >  3:152162/1‑149 (MQ=255)
                                                       aGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGgc                                                                                              >  2:137882/1‑147 (MQ=255)
                                                       aGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGgc                                                                                              <  1:137882/147‑1 (MQ=255)
                                                           cTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTgg                                                                                        <  1:107282/149‑1 (MQ=255)
                                                           cTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTgg                                                                                        <  5:268629/149‑1 (MQ=255)
                                                            ttGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGc                                                                                       >  3:28596/1‑149 (MQ=255)
                                                                ccTGGGAGCCGGAGAGTGAAGCCAATCGCCTGGTCGCCGGAATAGTCCGACAGGTAGCCACGCGCCTGGGCAATACTCGGGGAGTGTGCCGGCGCTGCTAAATCCTCCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAt                                                                                   <  1:219116/149‑1 (MQ=255)
                                                                    ggAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAttt                                                                                 <  1:304891/147‑1 (MQ=255)
                                                                        ccGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGa                                                                            >  3:156430/1‑148 (MQ=255)
                                                                           gagaGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTCGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCCGCAGTTCAGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGaaaaa                                                                        >  3:226610/1‑149 (MQ=255)
                                                                                tGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCgt                                                                    <  6:63666/148‑1 (MQ=255)
                                                                                   aGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTgc                                                                >  4:163263/1‑149 (MQ=255)
                                                                                         aCGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAg                                                          <  4:152162/149‑1 (MQ=255)
                                                                                             cAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCt                                                      <  7:240250/149‑1 (MQ=255)
                                                                                             cAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCt                                                      >  7:107487/1‑149 (MQ=255)
                                                                                              aGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTg                                                     >  3:212258/1‑149 (MQ=255)
                                                                                                gTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGa                                                   >  6:93537/1‑149 (MQ=255)
                                                                                                 tCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCgtgt                                                                  >  8:171023/1‑133 (MQ=255)
                                                                                                 tCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCgtgt                                                                  <  7:171023/133‑1 (MQ=255)
                                                                                                      gcAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCgga                                             <  4:28596/149‑1 (MQ=255)
                                                                                                               cGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTgg                                     <  7:284213/148‑1 (MQ=255)
                                                                                                                     gTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCt                              <  2:141767/149‑1 (MQ=255)
                                                                                                                        gCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCtgtt                           <  8:107487/149‑1 (MQ=255)
                                                                                                                        gCCACGCGCCCGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGCGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCtgt                            <  4:269739/148‑1 (MQ=255)
                                                                                                                          cACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCg                                                                             >  2:284762/1‑97 (MQ=255)
                                                                                                                          cACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCg                                                                             <  1:284762/97‑1 (MQ=255)
                                                                                                                                ccTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTttcttc                   <  4:175327/149‑1 (MQ=255)
                                                                                                                                      cAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCATTTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCAc             <  3:163263/149‑1 (MQ=255)
                                                                                                                                      cAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCAc             >  2:273484/1‑149 (MQ=255)
                                                                                                                                         tACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTATTAAGGCACtt           >  7:162783/1‑148 (MQ=255)
                                                                                                                                         agcGCCGGCACCCTGCCGGCGGTGCTACCTCCACCCGGCAGTTCTGGAGCACTCTGGCTTGGGGCGCCTGGCGAATTTGCCGAATAACCGTGTGCGCAACGGCCTGGAGCCGGAGGCAGTGGCCTGGCGGGTGTTTCTTCCGGCACGTg          >  8:300261/3‑149 (MQ=255)
                                                                                                                                         agcGCCGGCACACTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTg          >  7:37093/3‑149 (MQ=255)
                                                                                                                                                 cAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACag  <  5:138132/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACTGCCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAG  >  NC_002947/4357349‑4357642

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