Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,478,845 G→C 9.8% intergenic (+98/+347) PP_2172 → / ← PP_2173 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,478,8450GC 9.8% 147.0 / 5.2 51intergenic (+98/+347)PP_2172/PP_2173hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (18/26);  new base C (2/3);  total (20/31)
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.32e-01

CTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGAT  >  NC_002947/2478706‑2478992
                                                                                                                                           |                                                                                                                                                   
cTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGcgctcg                                                                                                                                             >  6:295081/1‑143 (MQ=255)
   aGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAgg                                                                                                                                         <  3:206754/149‑1 (MQ=255)
   aGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAgg                                                                                                                                         <  6:171476/149‑1 (MQ=255)
   aGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAg                                                                                                                                          <  3:118040/148‑1 (MQ=255)
     cAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgt                                                                                                                                       <  2:216044/149‑1 (MQ=255)
      aaCGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgtg                                                                                                                                      <  4:199494/149‑1 (MQ=255)
               gggTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTgttttgtttt                                                                                                                             >  3:124228/1‑149 (MQ=255)
               gggTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTgttttgtttt                                                                                                                             <  6:127166/149‑1 (MQ=255)
                         cAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCTCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGAcc                                                                                                                      <  1:45595/146‑1 (MQ=255)
                              aaTGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCAAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAc                                                                                                              <  2:228490/149‑1 (MQ=255)
                                   tCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagc                                                                                                          <  4:172820/148‑1 (MQ=255)
                                   tCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagc                                                                                                          >  1:228633/1‑148 (MQ=255)
                                   tCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagc                                                                                                          >  3:172820/1‑148 (MQ=255)
                                    cGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagc                                                                                                          >  8:68297/1‑147 (MQ=255)
                                       tgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTc                                                                                                     <  4:124228/149‑1 (MQ=255)
                                       tgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTc                                                                                                     <  3:250211/149‑1 (MQ=255)
                                                cGCATCCATGCATTACCGGCTGTAGAGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCGCCCACAACCTGCTAGCGCTTAGCTTTCGCAGTGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTgg                                                                                            <  3:293718/149‑1 (MQ=255)
                                                  cATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAg                                                                                                   >  5:92539/1‑140 (MQ=255)
                                                  cATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAg                                                                                                   <  6:92539/140‑1 (MQ=255)
                                                   aTCCATGAAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCACCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATAACGAGCGAGTCAGCACGTGGCgg                                                                                         >  7:68966/1‑149 (MQ=255)
                                                      cATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgat                                                                                      <  3:143609/149‑1 (MQ=255)
                                                               ccGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCg                                                                              <  1:220376/148‑1 (MQ=255)
                                                                cGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCg                                                                              <  2:250940/147‑1 (MQ=255)
                                                                cGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGt                                                                            <  3:241300/149‑1 (MQ=255)
                                                                  gCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTaa                                                                          >  3:290295/1‑149 (MQ=255)
                                                                    tGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcatc                                                                                                                   >  2:223658/1‑106 (MQ=255)
                                                                    tGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcatc                                                                                                                   <  1:223658/106‑1 (MQ=255)
                                                                      tAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACca                                                                       >  1:61433/1‑148 (MQ=255)
                                                                         cGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGt                                                                   >  4:115351/1‑149 (MQ=255)
                                                                                  gagCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCgg                                                          >  2:191428/1‑149 (MQ=255)
                                                                                     cTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCAt                                                       >  7:109378/1‑149 (MQ=255)
                                                                                        cgcgTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagc                                                                                                          <  1:263942/95‑1 (MQ=255)
                                                                                        cgcgTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAcgagc                                                                                                          >  2:263942/1‑95 (MQ=255)
                                                                                        cgcgTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATgg                                                     <  1:191428/148‑1 (MQ=255)
                                                                                              ggCCGTGCTAATTGAACGGGTTACCCCGCACACGCTGCGAGCGCTTAGCTTTCGCAGTGTGTTTTGTTTTGGAGACGATCATCACGAGAGAGTCAGCACGTGGCGGGATCGATCTCTGTAACCACAGTTACCTCCGGCATGTCCAAGGt                                              <  4:290295/149‑1 (MQ=255)
                                                                                                ccGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGTCATGGCCAAGGTAc                                            <  4:123833/149‑1 (MQ=255)
                                                                                                ccGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTAc                                            <  3:173624/149‑1 (MQ=255)
                                                                                                    gCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGCCCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTg                                        >  5:64591/1‑149 (MQ=255)
                                                                                                    gCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTg                                        >  4:116602/1‑149 (MQ=255)
                                                                                                      tAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACTATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTAc                                                                <  1:79461/123‑1 (MQ=255)
                                                                                                      tAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTAc                                                                >  2:79461/1‑123 (MQ=255)
                                                                                                            aaCGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGc                                <  7:68297/149‑1 (MQ=255)
                                                                                                                    aCCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGttt                         <  6:64591/148‑1 (MQ=255)
                                                                                                                    aCCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGttt                         <  3:116602/148‑1 (MQ=255)
                                                                                                                     ccccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTc                        <  8:109378/148‑1 (MQ=255)
                                                                                                                       ccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGttttgttttg                                                                                                                            >  2:144493/1‑46 (MQ=255)
                                                                                                                       ccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGttttgttttg                                                                                                                            <  1:144493/46‑1 (MQ=255)
                                                                                                                            cacCCTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGAGCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaagaag                >  2:157509/1‑149 (MQ=255)
                                                                                                                                   cGAGCGCTGATCTTTCTCAGGGTGTTTTGTTTTGGAGACCATCATCACGATCGAGTCAGCACTTGGCGTGATCGATCTCGGTAACCACAGTTACCTCCGGCATTGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGAc          <  8:68966/148‑1 (MQ=255)
                                                                                                                                   cGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGa           <  2:61433/147‑1 (MQ=255)
                                                                                                                                          tGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTAGCAAGAAGGTAGACGGCACGAt  >  7:240680/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGAT  >  NC_002947/2478706‑2478992

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: