Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,742,779 G→A 8.5% K12K (AAG→AAA)  PP_3307 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,742,7790GA 8.5% 127.1 / 3.4 47K12K (AAG→AAA) PP_3307hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base A (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

TGGCAGTGGCTCAGCCGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAACACGCAGGCTCGGGAAGCCGGACACCGTCGCA  >  NC_002947/3742630‑3742928
                                                                                                                                                     |                                                                                                                                                     
tggCAGTGGCTCAGCCGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTCTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCAGGGACAAAGACAAGTACACGGACAAg                                                                                                                                                       <  1:2240369/150‑1 (MQ=255)
  gCAGTGGCTCAGCCGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCa                                                                                                                                                     <  1:2323269/150‑1 (MQ=255)
       ggCTCAGCCGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAAc                                                                                                                                                >  2:200718/1‑150 (MQ=255)
          tCAGCCGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTa                                                                                                                                             >  2:2123187/1‑150 (MQ=255)
               cGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGc                                                                                                                                                      >  2:294572/1‑136 (MQ=255)
               cGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGc                                                                                                                                                      <  1:294572/136‑1 (MQ=255)
                 gAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGa                                                                                                                                      <  1:2123187/150‑1 (MQ=255)
                 gAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAACAGAAACGTAAGGCCGa                                                                                                                                      >  2:1957632/1‑150 (MQ=255)
                  aaGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAg                                                                                                                                     <  1:200718/150‑1 (MQ=255)
                     tGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGcac                                                                                                                                  >  1:1341495/1‑150 (MQ=255)
                           tGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAg                                                                                                                            <  1:1220188/150‑1 (MQ=255)
                               gcAGGTTCCGGCAGTTAGCTGAAAGGGCATCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACTTACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAAt                                                                                                                        >  1:509121/1‑150 (MQ=255)
                                              tAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAg                                                                                                         <  2:72422/150‑1 (MQ=255)
                                                  tGAAAGGGCCTCAGTTCTACCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGt                                                                                                     >  2:1365374/1‑150 (MQ=255)
                                                    aaaGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTgg                                                                                                   <  1:491141/150‑1 (MQ=255)
                                                    aaaGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTgg                                                                                                   <  2:1341495/150‑1 (MQ=255)
                                                          ccTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAACAGAAACGTAAGGCCGAGCACATCGAg                                                                                                                            <  2:1716171/119‑1 (MQ=255)
                                                          ccTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAACAGAAACGTAAGGCCGAGCACATCGAg                                                                                                                            >  1:1716171/1‑119 (MQ=255)
                                                            tCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAACAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGAt                                                                                           <  1:1957632/150‑1 (MQ=255)
                                                                   tGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCg                                                                                    >  1:1277697/1‑150 (MQ=255)
                                                                     cccTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGaa                                                                                  >  2:2275862/1‑150 (MQ=255)
                                                                     cccTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGaa                                                                                  >  2:2275521/1‑150 (MQ=255)
                                                                     cccTGGTTTTCAACCGTTCATGCAACCCATGCAGCAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGaa                                                                                  >  1:803121/1‑150 (MQ=255)
                                                                       cTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGCCAAGTACACGGCCAAGCAGAAACGTAAGACCGAACACAACGAGCAATCTGACGAGGACACCGGCCTGGGCACGGATGAAGCCGACGc                                                                                >  1:444511/1‑150 (MQ=255)
                                                                       cTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATg                                                                                          >  2:1247580/1‑140 (MQ=255)
                                                                       cTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATg                                                                                          <  1:1247580/140‑1 (MQ=255)
                                                                           ttttCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTcgc                                                                            >  1:1430703/1‑150 (MQ=255)
                                                                                        aTGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGACGAAGCTCGAGGCTGGGCCACCg                                                               >  2:2244414/1‑150 (MQ=255)
                                                                                         gatctcccATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGt                                                              <  2:1623625/145‑1 (MQ=255)
                                                                                               ccATGCAGGAGGAAATACCATGCCACGCGGGGAAAACCAGAAGTGCCGCCTGAAGCAGAAACGTAAGGCCTAGCACATCGAGCAATCTTACGACGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTcagatc                                                        >  1:1623625/1‑146 (MQ=255)
                                                                                                cATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCaacaaa                                                       >  2:2039245/1‑150 (MQ=255)
                                                                                                       gaggaAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAg                                                                                                                >  1:1744755/1‑86 (MQ=255)
                                                                                                       gaggaAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAg                                                                                                                <  2:1744755/86‑1 (MQ=255)
                                                                                                         ggaAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTgg                                                                      >  1:2044688/1‑126 (MQ=255)
                                                                                                         ggaAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTgg                                                                      <  2:2044688/126‑1 (MQ=255)
                                                                                                              tACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGaaa                                                                                                                                                 >  2:1810118/1‑46 (MQ=255)
                                                                                                              tACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGaaa                                                                                                                                                 <  1:1810118/46‑1 (MQ=255)
                                                                                                                       aCGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAaca                                >  2:493190/1‑150 (MQ=255)
                                                                                                                         gcgGGGAGAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGACCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAACACg                              <  2:509121/150‑1 (MQ=255)
                                                                                                                              gACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATc                                                                                                                               <  1:605808/48‑1 (MQ=255)
                                                                                                                              gACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATc                                                                                                                               >  2:605808/1‑48 (MQ=255)
                                                                                                                                 aaaGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGa                                                                                            >  2:1794466/1‑80 (MQ=255)
                                                                                                                                 aaaGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGa                                                                                            <  1:1794466/80‑1 (MQ=255)
                                                                                                                                 aaaGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGAGAGTATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAACACGCAGGCTCg                      <  1:1365374/150‑1 (MQ=255)
                                                                                                                                  aaGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACaaa                                                                                                          <  1:380733/65‑1 (MQ=255)
                                                                                                                                  aaGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACaaa                                                                                                          >  2:380733/1‑65 (MQ=255)
                                                                                                                                           acacGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAACACGCAGGCTCGGGAAGCCGGa            <  1:1127504/150‑1 (MQ=255)
                                                                                                                                                  cAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAACACGCAGGCTCGGGAAGCCGGACACCGTc     <  2:1430703/150‑1 (MQ=255)
                                                                                                                                                     gCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAACACGCAGGCTCGGGAAGCCGGACACCGTCGCa  >  2:350924/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCAGTGGCTCAGCCGGAAGTGGATCTGGCGCAGGTTCCGGCAGTTAGCTGAAAGGGCCTCAGTTCTGCCCTGGTTTTCAACCGTTCATGCAACCCATGCAGGAGGAAATACCATGCCACGCGGGGACAAAGACAAGTACACGGACAAGCAGAAACGTAAGGCCGAGCACATCGAGCAATCTTACGAGGACAAAGGCGTGGGCAAGGATGAAGCCGAAGCTCGCGCCTGGGCCACCGTCAACAAACAATCCGGTGGCGGTGAAAAACACGCAGGCTCGGGAAGCCGGACACCGTCGCA  >  NC_002947/3742630‑3742928

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