Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,224,546 T→G 7.3% intergenic (‑614/+207) PP_3703 ← / ← PP_3704 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,224,5460TG 7.3% 157.3 / 4.7 55intergenic (‑614/+207)PP_3703/PP_3704hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/30);  new base G (2/2);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGG  >  NC_002947/4224401‑4224691
                                                                                                                                                 |                                                                                                                                                 
cccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCtgt                                                                                                                                                 <  1:1175279/148‑1 (MQ=255)
cccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCtgt                                                                                                                                                 >  2:1175279/1‑148 (MQ=255)
 ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCa                                                                                                                                              >  1:1250927/1‑150 (MQ=255)
 ccTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCa                                                                                                                                              >  2:263856/1‑150 (MQ=255)
  cTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCAGCTGTGCAt                                                                                                                                             >  1:601598/1‑150 (MQ=255)
   tATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCAtt                                                                                                                                            <  2:844527/150‑1 (MQ=255)
    aTCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTg                                                                                                                                           >  2:686427/1‑150 (MQ=255)
      cAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCt                                                                                                                                         <  1:686427/150‑1 (MQ=255)
               cGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCAc                                                                                                                                >  1:232010/1‑150 (MQ=255)
               cGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCAc                                                                                                                                >  1:231917/1‑150 (MQ=255)
                      gTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                         <  1:1113177/150‑1 (MQ=255)
                      gTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTaaa                                                                                                                         <  1:614356/150‑1 (MQ=255)
                         ccTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                      >  2:808587/1‑150 (MQ=255)
                             aGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                  <  2:231917/150‑1 (MQ=255)
                             aGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTAt                                                                                                                  <  2:232010/150‑1 (MQ=255)
                               ccTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAGCCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGa                                                                                                                <  2:1339748/150‑1 (MQ=255)
                                 tGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGc                                                                                                              <  1:263856/150‑1 (MQ=255)
                                  gggagcggGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCt                                                                                                             >  2:452363/5‑150 (MQ=255)
                                        ggggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGAc                                                                                                       >  1:23880/1‑150 (MQ=255)
                                         gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACa                                                                                                      <  1:224014/150‑1 (MQ=255)
                                         gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACa                                                                                                      <  1:911780/150‑1 (MQ=255)
                                         gggAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACa                                                                                                      <  2:673286/150‑1 (MQ=255)
                                                acCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTgg                                                                                               <  2:23880/150‑1 (MQ=255)
                                                 ccTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTggg                                                                                              >  2:757811/1‑150 (MQ=255)
                                                   tCGCTCCCGTAAGAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGGGTACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTa                                                                                            <  1:664269/150‑1 (MQ=255)
                                                     gCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATc                                                                                          >  2:1005928/1‑150 (MQ=255)
                                                       tCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATcgc                                                                                        >  1:210252/1‑150 (MQ=255)
                                                          cGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGcc                                                                                     <  2:439472/150‑1 (MQ=255)
                                                           gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGt                                                                                             <  1:716071/141‑1 (MQ=255)
                                                           gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGt                                                                                             >  2:716071/1‑141 (MQ=255)
                                                           gTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGccc                                                                                    >  1:868247/1‑150 (MQ=255)
                                                                 aaGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGt                                                                              >  1:1311242/1‑150 (MQ=255)
                                                                     cGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCg                                                                          <  2:255862/150‑1 (MQ=255)
                                                                       gCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGga                                                                        <  1:1005928/150‑1 (MQ=255)
                                                                       gCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGga                                                                        <  2:868247/150‑1 (MQ=255)
                                                                              aCGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGCTGATa                                                                 <  2:1091509/150‑1 (MQ=255)
                                                                               cGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATAc                                                                <  2:601598/150‑1 (MQ=255)
                                                                               cGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATAc                                                                <  2:210252/150‑1 (MQ=255)
                                                                                gTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACa                                                               >  2:1049983/1‑150 (MQ=255)
                                                                                    tAAGCGGCGCGACGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAAc                                                           <  2:969154/150‑1 (MQ=255)
                                                                                          gAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                      <  2:1198862/85‑1 (MQ=255)
                                                                                          gAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGc                                                                                                                      >  1:1198862/1‑85 (MQ=255)
                                                                                                  cGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCa                                             <  1:757811/150‑1 (MQ=255)
                                                                                                   gggAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                                            >  2:237992/1‑150 (MQ=255)
                                                                                                    ggAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATg                                           <  1:452363/150‑1 (MQ=255)
                                                                                                            ccGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTg                                   >  1:1131041/1‑150 (MQ=255)
                                                                                                                   aaaTAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc                            <  2:1311242/150‑1 (MQ=255)
                                                                                                                   aaaTAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgc                            <  1:237992/150‑1 (MQ=255)
                                                                                                                     aTAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAgcgc                          >  2:506605/1‑150 (MQ=255)
                                                                                                                      tAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGcc                         >  2:616123/1‑150 (MQ=255)
                                                                                                                             aGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaaggaag                  <  1:616123/150‑1 (MQ=255)
                                                                                                                             aGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCaaggaag                  <  1:506605/150‑1 (MQ=255)
                                                                                                                                  aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                                            >  1:578595/3‑119 (MQ=255)
                                                                                                                                  aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCAt                                            <  2:578595/117‑1 (MQ=255)
                                                                                                                                  aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTg                                   <  2:842570/126‑1 (MQ=255)
                                                                                                                                  aaccAGCAATGCACAGGTGCACTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTg                                   >  1:842570/3‑128 (MQ=255)
                                                                                                                                             cACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAggg  <  1:1049983/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCCTATCAATCCAAACGTAAATGTCCCTAAGCCTGGCGGCGGGGAACAACCTCGCTCCCGTAACAAAGACGGCCAGATACGTGATAAGCGGAGCGATGCGGGAAAACCCCGCGGAAAATAAACCTAGGTGGCCCAGCAGTGCACCTGTGCATTGCTGGTTGCCACGTGTAAAGGCTTATGAGCTAAAGACAGGACTGGGTATCGCGCCCATCTGTAGCGGAGATGATACAGAACCTTAGGATTGTCCATGTAAACCTGGAAAAGCGCCAAGGAAGGTTCTGTGGAGAAGGG  >  NC_002947/4224401‑4224691

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